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Results computed in lattice QCD+QED are presented for the electromagnetic 
mass splittings of the low lying hadrons. These are used to determine the renormal- 
ized, non-degenerate, light quark masses. It is found that mff s = 2.24 (10) (34), 
m js = 4.65 (15) (32), and mf^ = 97.6 (2.9) (5.5) MeV at the renormaliza- 
tion scale 2 GeV, where the first error is statistical and the second systematic. 
We find the lowest order electromagnetic splitting {m n + — tti^qed = 3.38(23) 
MeV, the splittings including next-to-leading order, (m n + — m 7r o)QED = 4.50(23) 
MeV, (m K + — m K o)qET> = 1.87(10) MeV, and the m u / m& contribution to the 
kaon mass difference, {m K + — rn K o)^ niu _ md ) = —5.840(96) MeV. All errors are 
statistical only, and the next-to-leading order pion splitting is only approximate 
in that it does not contain all next-to-leading order contributions. We also com- 
puted the proton-neutron mass difference, including for the first time, QED inter- 
actions in a realistic 2+1 flavor calculation. We find (m p — TOjJqed = 0.383(68) 
MeV, (m p — m n )( mu _ md ) = —2.51(14) MeV (statistical errors only), and the total 
2.13(16)(70) MeV, where the first error is statistical, and the second, 
part of the systematic error. The calculations are carried out on QCD ensembles 
generated by the RBC and UKQCD collaborations, using domain wall fermions and 
the Iwasaki gauge action (gauge coupling (3 = 2.13 and lattice cutoff a -1 rj 1.78 
GeV). We use two lattice sizes, 16 3 and 24 3 ( (1.8 fm) 3 and (2.7 fm) 3 ), to address 
finite volume effects. Non-compact QED is treated in the quenched approximation. 
The valence pseudo-scalar meson masses in our study cover a range of about 250 
to 700 MeV, though we use only those up to about 400 MeV to quote final results. 
We present new results for the electromagnetic low energy constants in SU(3) and 
SU(2) partially-quenched chiral perturbation theory to the next-to-leading order, 
obtained from fits to our data. Detailed analysis of systematic errors in our results 
and methods for improving them are discussed. Finally, new analytic results for 
SU(2)l x SU(2)R-plus-kaon chiral perturbation theory, including the one-loop logs 
proportional to a em m, are given. 
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I. INTRODUCTION 

The mass splitting in the meson and baryon systems is an interesting topic in hadron 
spectroscopy. It is related to the quark masses which are fundamental parameters of the 
standard model. The mass splitting in the pseudo-scalar meson octet is a signature of 
the breaking of the strong isospin symmetry by the electromagnetic (EM) interaction and 
non-degenerate quark masses. 

The hadron spectra are rich in diversity due to two origins: the nonperturbative quan- 
tum dynamics of the strong interaction and the presence of flavor symmetry breaking. In 
the standard model, the latter originates from non-degenerate quark masses as well as the 
difference between up-type quarks and down-type quarks. These sources of flavor symmetry 
breaking affect significantly the hadron spectra less than (1 ~ 2) GeV. In the baryon octet, 
for instance, the mass difference of the proton and neutron is crucial to the phenomenological 
model of nuclei because it plays an important role in neutron /3-decay, which is related to 
the stability of nuclei. If the up and down quark masses were degenerate, the proton would 
be heavier due to the EM interaction, and our Universe would not exist! Even though the 
mass differences in the baryon octet spectrum have already been measured in experiments to 
good accuracy, it is important to confirm that we can predict these splittings in the standard 
model from first-principles computation. Parameterizing the calculated splittings in terms 
of low energy constants (LEC) is also useful for effective theories like chiral perturbation 
theory. 

The mass of a hadron is determined by both quantum chromodynamics (QCD) and 
quantum electrodynamics (QED), though the vast majority of the mass is due to QCD. 
For the QCD interaction, since the coupling constant is large in the low energy regime 
(E < 1 GeV), perturbation theory is not applicable, and we must turn to the techniques 
of lattice gauge theory to solve QCD for the hadron spectrum. On the other hand, the 
EM contributions to the masses, which break degeneracies due to flavor SU(3) (and isospin) 
symmetry since the up quark has charge +2/3 e and the down and strange quarks have 
charge —1/3 e, are expected to depend on the small QED coupling constant, a em = e 2 /47r 
1/137. However, since the hadrons are formed from bound-states of the quarks, there is 
no systematic way to treat the contributions in weak-coupling perturbation theory. Thus, 
calculations are done nonperturbatively in a combined lattice QCD+QED theory (indeed, 
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even if the strong coupling constant were small, one is forced to a nonperturbative solution for 
QCD because of confinement). The state-of-the-art in lattice calculations is such that sub- 
percent errors (statistical and systematic) on low-lying hadron masses and other observable 
quantities are becoming the norm(for a broad review, see pQ and references therein). As the 
precision of lattice calculations improve, the EM splittings become more and more relevant. 
Indeed, the splittings themselves can be computed with sub-percent precision, at least for 
the statistical errors[2]. 

It is well known that the lowest order EM effect, the so-called Dashen term[3], which 
enters at 0(a em ), is the dominant contribution to the charged- neutral pion mass difference. 
In the chiral limit where the quarks are all massless, it is also true for the kaons. This 
theorem, known as Dashen's theorem, is broken by terms of order 0{a em m) away from the 
chiral limit. Using an effective theory of QCD known as chiral perturbation theory, these 
corrections can be identified in the lattice calculation, and the non-degenerate quark masses 
determined by matching to the experimentally measured mass splittings [I], 

m n ± - m,o = 4.5936(5) MeV, (1) 
m K ± - m K o = -3.937(28) MeV, (2) 
m n -m p = 1.2933321(4) MeV. (3) 

In fact, any three hadron masses are enough to determine the three light quark masses, and 
we choose m n +,mK+ and m^o for reasons explained later. The determination of the up 
quark mass, m u , is particularly interesting since one can check the simplest solution to the 
strong CP problem, m u = 0. 

In the study presented here, we work with lattice domain wall fermions (DWF) [5l E] for 
the quarks and the Iwasaki gauge action for the gluons. We use an ensemble of gluon con- 
figurations with a single lattice spacing, generated by the RBC and UKQCD collaborations 
using (2+l)-flavors of dynamical quarks, a pair of degenerate quarks for the up and down, 
and a heavier strange quark [7|. The photons are simulated in non-compact, quenched QED, 
as was done in the pioneering quenched QCD+quenched QED calculations [SJ E] • 

There are several differences between this work and our previous one j2]. The most 
obvious is that the dynamical strange sea quark has been included here. In Ref. [2] the QED 
gauge potential was fixed to the Coulomb gauge, and here we work in Feynman gauge in QED 
on finite volume, as described in Ref. [10]. Next, we perform fits to full next-to-leading order 



6 



(NLO) partially-quenched chiral perturbation theory (PQ%PT), including photons, for both 
SU(3) L x SU(3) R chiral symmetry [IT] and SU(2) L x SU(2) R -plus-kaon 0II21CE3], where the 
latter treats only the pion-triplet as pseudo-Nambu-Goldstone bosons. The NLO PQ^PT 
for SU(2)l x SU(2)R-plus-kaon, including photons, is new, and is presented here for the first 
time. Calculations by the RBC and UKQCD collaborations [7J IX1HTT] . and more recently 
PACS-CS[18l [19], have shown that the physical strange quark mass is out of reach for NLO 
SU(3) chiral perturbation theory. Since we also wish to determine the strange quark mass in 
our calculation, we have developed the chiral perturbation theory for the SU(2) L x SU(2) R - 
plus-kaon system, including photons. In addition, since the photons are not confined, finite 
volume effects are expected to be large, so we work with two lattice sizes, 16 3 and 24 3 , with 
the same lattice spacing to investigate these effects. The leading EM finite volume effects 
have been computed in PQxPT [10], which we also use in our analysis. 

This paper is organized as follows. In Sec. [TTJ we summarize the chiral perturbation 



theories used to fit our lattice calculations (details are given in the Appendix Bl). In Sec. Ill 



the basic framework and details of the lattice simulations are given. Section IV contains 
results and discussion of the calculation, including the fitted LEC's and the quark masses. 



Section [V] discusses systematic errors, and in Sees. VI and VII we give final values for the 
quark masses and meson splittings, respectively. We examine the impact of (m u — ma) on 



the decay constant ratio, fx/fir in Sec. VIII The nucleon mass splitting is computed in 
Sec. |IX| Finally, this work is summarized in Sec. |X} 

We reported the preliminary results from SU(3) PQ^PT study in this work in Ref. [201 
I2T] . The MILC collaboration also presented their first results on the EM splittings, using 
improved staggered fermions and non-compact, quenched QED configurations in Ref. [22J. 



II. CHIRAL PERTURBATION THEORY 

We briefly review the framework and formulas of partially quenched chiral perturbation 
theory relevant for our 2+1 flavor calculation. The EM corrections in SU(2) chiral pertur- 
bation theory coupled to kaons are new. Details are given in Appendix [Bj 

Recently it has been shown that SU(3) chiral perturbation theory is poorly convergent for 
quark masses near the physical strange quark mass, and that a straightforward and effective 
solution is to treat the strange quark mass m s as large compared to the light quark masses 
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mi in an expansion in mi/m s [H [12j [TH [19]. We carry out fits to the data using both SU(3) 
and SU(2) chiral perturbation theory. We find the poor convergence extends to the EM 
sector as well, and use SU(2) chiral perturbation theory to quote our central results. 

Before proceeding, it is important to discuss the order in chiral perturbation theory to 
which we work. For the SU(3) case where the kaon is a Nambu-Goldstone boson, the leading 
order (LO) includes all terms that are 0(p 2 ) and 0(e 2 ), and the next-to-leading order (NLO) 
includes all terms that are 0(p 4 ), 0(e 2 p 2 ), and 0(e 4 ), where the conventional power counting 
is 0(e) ~ 0(p). This counting is the same for the square of the masses and the mass-squared 
splittings. 0(e 4 ) contributions have so far been ignored [21 |TT| since 0(e 4 ) <C 0(e 2 p 2 ) in 
practice, and we also follow this here. 

In the SU(2) theory coupled to kaons the power counting becomes a bit more complicated 
for the kaon (for the pion it is the same as in the SU(3) case). Since the kaon is no longer 
a Nambu-Goldstone boson, LO for the mass-squared is now 0(p°), and NLO is 0(p 2 ) and 
0(e 2 ). The mass-squared EM splitting, however, remains the same order of magnitude 
as in (partially-quenched) SU(3) chiral perturbation theory. That is, to obtain the NLO 
contributions to the mass-squared splittings, we must work to NNLO for the masses. Since 
the aim here is to include all effects up to and including 0(e 2 p 2 ) terms in the meson mass- 
squared splittings as well as — m u , we include all 0(e 2 p 2 ) contributions to the kaon 
mass. Because we compute m^ — m u from the neutral-charged kaon mass-squared difference, 
the pure QCD effects at 0(p A ), including the one-loop logarithms [12], cancel and are not 
included in our analysis. 

Finally, to avoid confusion we emphasize that in this paper we only calculate correlation 
functions for "charged", or "off-diagonal" mesons. However, since we are free to change 
the charges and masses of the valence quarks making up these mesons, the total charge of 
the (unphysical) meson may happen to be zero. Sometimes we refer to these as "neutral" 
mesons, but it must be kept in mind these never correspond to the 7r° meson which requires 
so-called disconnected quark diagrams in its correlation function as well as the full treatment 
of "diagonal" mesons in PQxPT. 
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A. SU(3) L x SU(3) R 

The partially quenched chiral perturbation theory has been worked out in Ref. [TT] . 
and we adopt their notation. For three non-degenerate sea quarks and two non-degenerate 
valence quarks, labeled by "1" and "3" , the meson mass-squared at NLO is 

M 2 3 = X i3 + -e^Qis 

A8L r 6 -2ALl _ 16L r 8 -8L r 5 2 

^2 Xl3Xl H ^2 Xl3 

-48e 2 ^L^ 1 2 3 x 1 - 16e 2 §- 4 Llqj 3Xl3 

-e 2 [l2K? r + YlKf- - \2K* r - YlKf r \ q 2 Xl3 

-e 2 [4K? + AKt] q 2 pXl3 

+e 2 [4K 9 Er + AK%] q 2 pXp 

+12e 2 K* r q 2 13Xl 

+8e 2 [Kg + Kg] g 2 3 Xi3 

-e 2 [8Kg + AK%] qi q3 X i3 

■1 ^ pm l Xm . 1 1 r>» 1 Xp 

^4 7«~~2 Xi4 log —2-914 + Xi5 log — o-9i5 + Xi6 log — rfte 913 

Fq 167T^ \ yLT yU^ ylT / 

,2^1/ i X34 i X35 , X36 \ 

+2e "^4 7fi~^ *34 l0 § ~<?34 + X35 log —935 + X36 log "936 9l3 
Xq 107T \ fjj pb fjj J 

+e 2 5 mrcs (9i 2 + 9 3 2 )- (4) 



+ 3 
-2e 2 



Indices 1 — 3 always refer to valence quarks, 4 — 6 to sea quarks. The coefficients -R™ 13 and 
Rq nv are the residue functions written in terms of quark masses and are defined in Ref. [TT] . 
The index p implies summation over valence indices 1 and 3, and if q is also present, then 
the sum is over pairs (1,3) and (3,1). The indices (m,n) signify a sum over pairs (tt, rj) and 
(rj, 7r). X ij = B (rrii + rrij) is the LO mass-squared for a meson made of quarks with masses 
rrii and rrij, qij = qi — qj where qi is the electric charge of the ith quark in units of the 
fundamental charge e. X i — 2-B q (to4 + m 5 + m 6 )/3 and q 2 = {q\ + q\ + 9§)/3. Xw and Xrj are 
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given by the solution of 

Xtt + X v = 2 Xu (5) 
4 

XwXv = 3 5o(m 4 m 5 + m 5 m 6 + m 4 m 6 ). (6) 

The LO Dashen term is proportional to the low energy constant (LEC) C and the fine 
structure constant a em . B and F [] are the LO QCD LEC's, the L's are the Gasser- 
Leutwyler LEC's at NLO, and the K's are the EM LEC's at 0(a era m). S mics is a pure 
lattice artifact LEC associated with the finite size of the extra dimension for DWF. 

We note from Eq. Q that masses of mesons ~ qq' made from degenerate valence quarks 
q, q' with equal charges do not have logarithmic corrections at NLO. This happens for the 
SU(2) case as well. 

Following Ref. [11], the EM LEC's can be written in terms of five independent linear 
combinations of the K's, which is all that can be determined from lattice calculations, 

Yi = K? r + K% r - K? r - Ki\ (7) 

Y 2 = K 9 Er + Jjf* (8) 

Y 3 = -K? - Kf + 2K% + 2K?{, (9) 
Y 4 = 1Kf r + 1Kf r + 2K* r + K* r , 

Y 5 = K§\ (10) 

The EM mass-squared splitting of the pseudo-scalar meson is defined as AM 2 = M 2 (e ^ 
0) — M 2 (e = 0). In terms of the Yj's, it becomes 
,o 2Cc o 



AM = -^q( 3 



F 2 



r ^0 



-12e 2 Yig 2 Xi3 + 4e 2 Y2g 2 x P + te 2 Y 3 q 2 l3X iz ~ ^e 2 Y iqi q 3X iz + l2e 2 Y 5 q 2 13X i 



-2e -— Xw log — gi4 + Xi5 log — 9i5 + Xie log — qie ) ?i3 

-Tq X07T \ jJL jJj jJj 

2 C 1 (. ,__%34_ , ,__X35_. , .. n __X36 



+ 2e "^4 7^ *34 log "^-<?34 + X35 log ^"935 + X36 log "^-<?36 <?13 

Fq 16ti 2 \ jj, 2 n 2 n 2 



(<?i3) 2 e 2 f /xis 



1«t> » 3 |3 1og^]-4 

,2r /„2 



+e 2 <W cs (g 2 + g 3 2 )- (11) 



1 Fq is normalized such that the physical decay constant is roughly 92 MeV. 
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Note that Y\ is proportional to the sea quark charges. Since we work with quenched QED, 
this LEC can not be obtained from our calculation. 



We carry out the fit in Section IV B 2 with the finite volume correction to the chiral 
logarithms taken into account. The finite volume correction to the leading-order chiral 
logarithms was computed in Ref. [10J, 



SM^L) = mUl) - MUoc) 



1 1 



^Cl3 Xl3 



M{Jx^L) t 1 1 



3 16tt 2 F 2 nU 



L 2 



+ 



R Lr, Xl3 



M{Jx P L) 



3 16n 2 F 2 qnv 



L 2 



-2e 2 



C 1 



<?13 



F 4 16tt 2 

x q u — h qis To ^ 9ie 



V- 



V- 



V- 



+2e 2 



C 1 



<?13 



F 4 167T 2 

( M{^XmL) M(^x^L) M(^XmL) 
x 934 T7, V ?35 tt, V qm 



L 2 



o (gi3) e 2 « 
4vr L 2 

(qn) 2 e 2 \K (VxTsL) 



±7T 



L 2 



13 



L 2 



L 2 



(12) 



M{x) is the function appearing in the finite volume correction to the chiral logarithm induced 
by the tadpole diagram, 



M(x) = 47T 



dX 
A 2 



eXp ("^ A ) T(A) 



T(A) = h? 3 0, i 



A 



(13) 



where $3(1*; r) is a Jacobi-theta function, 



$3 (v, r) = exp (nrin 2 + 2ixvin) 
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The other functions and a constant k are given by [10 



K 



f°° rl\ 

l ^ 5(A) = 2.837-.., (14) 



3 

3 



S{X) = -U»s[0,i\)) -1-Aij, (15) 

H(x) = vr ^ erf h^y^ 5(A) , (16) 

K{x) = 4x J o ^_(i- e -fc*)5(A), (17) 

where erf(x) is the error function, 



2 r 

erf (x) = —= / dse 
Jo 



B. SU(2) L x SU(2) R -plus-kaon 

Some time ago Roessl [32] worked out the low energy SU(2) Lagrangian of pions coupled 
to a kaon. Recently, the RBC and UKQCD collaborations showed that SU(3) chiral pertur- 
bation theory is poorly convergent for quark masses near the strange quark mass but that 
SU(2) chiral perturbation theory coupled to a kaon worked well for pions with masses less 
than about 400 MeV at NLO [7j. In Ref. [7], the unitary Lagrangian was extended to the 
partially quenched case. Here, we extend both works to include the EM interactions to the 
order a em m for the kaon mass, including the one-loop diagrams proportional to a em . For the 
pion mass, we begin with the partially-quenched SU(3) Lagrangian in Ref. [11] and expand 
in mi/m s . 



1. pions 



We derive the SU(2)l x SU(2)r result for the pion mass-squared splitting by expanding 
Eq. (11) in (mi, m.3, 777.4, rah)/™®, where rriQ is the strange sea quark mass, mi and 777,3 are 
taken as non-degenerate light valence quark masses, and m 4 and m 5 the light sea quark 
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masses, 
AM 2 = ^? 



^0 



-48e ^4^3 16e -^L 5 q 13 Xi3 

-12e 2 yig 2 xi3 + 4e 2 F 2 g 2 Xp + 4e 2 y 3 g 2 3 xi 3 - 4e 2 F 4 gig 3 Xi3 + 12e 2 F 5 g 2 3 ^^ 

2 3 Xl3 2 , 2 1 2 

- 6 16^ Xl3l0g ^ 3 + e 4^^ 3 



/i 

2 C 1 ( , Xl4 . , XlB i X34 , X35\ 

^?13 <?14Xl4 log —j- + gisXlS log — =- - g 3 4X34 log — ^ - 235X35 log — 3- 

+e 2 5 mrcs (g 2 + g 3 2 )- (18) 



In Eq. (18) all of the low energy constants now depend implicitly on the strange sea quark 
mass which is fixed (we rename them below to distinguish them from their SU(3) counter- 
parts). In addition the Dashen term has absorbed contributions from the NLO SU(3) LEC's 
and the logs which do not depend on the charges or masses of the up and down quarks, 

2Ce 2 2 , in 2 V 2 X& 2e 2 C 2 x& AQ 2 @ T r 2 X& / in \ 

-Ff** + 12e Y ^Y-T^F* q ^ hg ^- A8e F* L * q "T' (19) 

Including the contributions in pure QCD [TJ, the pion mass-squared to NLO becomes 
M 2 = Xra {l + ||(24 2) - L?)^ + ^(2L?> - z£>) Xu 

\ Tei^ (^ 3 X7r log | + ^ 3 Xl log I + ^ * 3 log I) } 



+ 2 



2C( 2 )e 2 2 
+ — ^2— ?13 

-10 2-ia(2) -2 , j 2xx( 2 ) 2 , a 2xr( 2 ) 2 . 2 V (2) . 10 2x^(2) 2 Xi + X5 

- 12e ; g xh + 4e y 2 v ; g p x P + 4e F 3 V ; g 13 Xis - 4e F 4 V >g 3 Xi3 + 12e K> ; g 13 — - — 

2 3 Xl3 2 , 2 1 2 

" e 16^ Xl3l0g ^ gi3 + e 4^ Xl3gi3 

2 C (2) 1 _ f_ _ ,__Xl4 , _ .. ,__Xl5 _ ., ,__X34 _ .. n __X35 



~ e ~EH ^<?13 <?14Xl4 log — + <?l 5 Xl5 log - ^34X34 log ~ <? 35 X35 log , 

F 4 87H \ ji 2 ji 2 fx 2 n 2 

+ e 2 6 mr Jq 2 1 +ql). (20) 

The LO LEC's F and B are the counterparts of Fq and i?o from the SU(3) theory, and 
the other SU(2) LEC's are denoted by an explicit superscript "(2)". R l - k is given in SU{2) 
partially quenched 

Ri = (Xi - Xa) (Xi - Xs) / 21 n 
Jfc ~ (Xi - X 3 ) (Xi ~Xk)' 
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The finite volume correction to Eq. (20) is given by 



5M 2 {L) = M 2 (L) - M 2 (oo) 



1 Xis ( R7T M (Vx^L) t d1 M (VxIL) t d3 M (VXiL) 



2 16vr 2 F 2 



13 



+ R 



■MfVxiIF) , M{Jx^L) 



" 2C 16^ gi3 t, 914 Z 2 " + 915 



-g 3 4 ^ <?35 



L 2 



L 2 
L 2 



o (gis) e 2 w 
4tt L 2 



(4tt)< 



L 2 



(22) 



The constant k and various functions are defined in Eqs. (13), (14), (16) and (17). 



2. kaons 

The kaon mass can be obtained from the tree-level Lagrangian, following Refs. [T2(, 123] . 
by constructing the kaon from one light and one "heavy" quark and writing down all op- 
erators with the desired symmetries in a non-relativistic theory where the power counting 
is straightforward. The needed relativistic Lagrangian is then constructed such that in the 
limit that the kaon is heavy, the non-relativistic theory is recovered. This has been done in 
the case of QCD to NNLO in Ref. [12] and to NLO in partially quenched QCD in Ref. [7]. 
Here we add the order e 2 p 2 terms induced by the EM interactions. Once the tree-level La- 
grangian is known, the one-loop corrections can be computed. The 0(e 2 ) Lagrangian and 
details of the one-loop calculation are given in the Appendix [Bj 

The 0(e 2 p 2 ) Lagrangian is quite complicated, with many operators appearing. While 
we have listed all possible operators in the Appendix [B] that contribute, we have not yet 
reduced them to a linearly independent set using relativistic invariance and the equations 
of motion. Still, this is enough to give the general quark mass and charge dependence. In 
the following, this is given by the generic LEC's £3 ~ 
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From Eqs. (B36), (B39) and (B42), the mass-squared of the kaon is 

M 2 K = M 2 - AB(A 3 m 1 + A 4 (m 4 + m 5 )) 

+e 2 (2 (a™ + 4 2 '") q ] + A'f'^l + 24*W) 

^ ' s=4,5 " 

+cV (xf + g 3 ) 2 + xf\ qi - q 3 ) 2 + xf \q\ - g 3 2 )) 
+e 2 ^ (4^1 + + x?\ qi - q ,f + ^(^ - <&) 

+e 2 ^m rcs (gi 2 + g 2 ), 



(23) 



where we have included the explicit chiral symmetry breaking LEC S mieB , the same as for 
the pion. Here the subscript "1" stands for a light valence quark, u or d, and "3" for the 
strange valence quark (charge). "4" and "5" refer to the u and d sea quarks, respectively. 
To avoid confusion we note that LEC's without superscripts denote the pure QCD LEC's of 
Refs. [12] and [7j while those with superscripts are EM LEC's defined in Appendix IB} The 



finite volume correction to Eq. (23) is given in Appendix B Eqs. (B40) and (B41) 



Notice that the LO "Dashen" term is different than for the pion: the latter is a single 
LEC proportional to q 2 3 while the former consists of three LEC's and depends on the u, d, 
and s charges separately. This is a consequence of the different chiral symmetries assumed 
in the two cases. 

We remind the reader that we do not keep terms of order p 4 and e 4 . 



III. LATTICE FRAMEWORK 



Following Ref. [8], the lattice calculation employs combined QCD+QED gauge configura- 
tions. A combined gluon-photon gauge link is simply the product of two independent links, 
a SU(3) color matrix for the gluons and a U(l) phase for the photons. 

U x „ = U$x(uV) Qt , (24) 

where Qi = eqi is the charge of the quark with flavor %. It is the combined link that appears 
in the lattice Dirac operator, in the usual gauge-invariant way. The gluon and photon links 
were generated independently in our calculation, so the sea quarks were not electrically 
charged. This quenched QED calculation suffers a systematic error that is expected to be 
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0(a em a s ) from a simple vacuum polarization argument. In chiral perturbation theory, the 
charged sea quarks first contribute at 0(a em m va i) for the valence quark mass m va i, as we have 
seen in Sec. |llj This drawback can be eliminated with the technique of re- weighting J23H2H], 
which is becoming common in large scale dynamical calculations [T5] - fTTl 127] , and is under 
active investigation by us [281 EH]. In a different context, combined dynamical simulations 
have also been performed for the first time [30], where the sea quarks are charged from the 
beginning. 

For the QCD configurations, we use the 2 + 1 flavor QCD configurations generated with 
DWF and the Iwasaki gauge action ((3 = 2.13) by the RBC and UKQCD collaborations [7J 
□U EI]. The lattice sizes are 16 3 x 32 and 24 3 x 64. The lattice spacingQis cT 1 = 1.784 (44) 
GeV, as determined from the Q baryon mass on the larger lattice, and which yields physical 
volumes (1.76 fm) 3 and (2.65 fm) 3 , respectively. The domain wall height M5 and the size of 
the extra dimension L s are 1.8 and 16, respectively. The residual quark mass in the chiral 
limit for pure QCD is found to be m res = 0.003148(46) and 0.003203(15), for the 16 3 and 24 3 
lattice sizes, respectively. The latter is slightly larger than the value 0.00315(2) determined 
in Ref. [7] on a smaller ensemble of configurations. 

The ensembles and number of measurements on each are summarized in Table HI The 
stopping criterion in the conjugate- gradient algorithm used to compute quark propagators 
was 10~ 8 , the same as in Ref. [7j. To increase our statistics on some of the ensembles, 
two or more different locations of the source are used on each configuration (see Table [T]). 
The QCD configurations are separated by 20 or 40 Monte Carlo time units to suppress the 
autocorrelations in them. Our calculation is for the pseudo-scalar meson at the unitary point 
(mi = W3 = m4 = 777.5) and the partially quenched point (arbitrary quark mass combination). 

The quenched QED configurations were generated on the non-compact manifold [21 [TO] . 
Here we employ the Feynman gauge instead of the Coulomb gauge which was used previously 
[2] in our two-flavor calculation. Since the mass is a gauge invariant quantity, the result 
should be consistent within the statistical error, up to the effects of zero-modes. Further, 
removal of the modes also results in the satisfaction of Gauss' Law on the torus [lOj . An 
advantage of the non-compact QED formalism is that the U(l) gauge potential can be 
chosen randomly with the correct Gaussian distribution in momentum space, then Fourier 
2 This result is slightly larger than the published one, 1.729 (28) GeV, in Ref. [7] because it was determined 
on a larger ensemble of lattices. It is also larger than the result of a combined fit, including new 32 3 , 
f3 = 2.25, ensembles [T5¥T7] . Later, we use the slight difference as a systematic error. 
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transformed to coordinate space, so there are no autocorrelations in the ensemble. Finally, 
yet another advantage is that there is no lattice-artifact photon self-interactions in the action. 
To couple to the fermions, the non-compact field is exponentiated to produce the photon 
link, U { x 11 — exp (ieA X)f j,), where e — a/ 47ra em « 0.30286. 

Since the QED interaction does not confine, it is possible that the finite volume may 
induce a significant systematic error. We thus do our simulation on both 16 3 and 24 3 lattice 
configurations with the same lattice spacing. This allows direct investigation of the finite 
volume effect in the mass spectrum. 



lat Tn sea 


myal 


Trajectories 


A -/V" meas 


^src 


16 3 0.01 


0.01, 0.02, 0.03 


500-4000 


20 


352 


4,20 


16 3 0.02 


0.01, 0.02, 0.03 


500-4000 


20 


352 


4,20 


16 3 0.02 


0.01, 0.02, 0.03 


500-4000 


20 


352 


4,20 


24 3 0.005 0.001, 0.005, 0.01, 0.02, 0.03 


900-8660 


40 


195 





24 3 0.01 


0.001, 0.01, 0.02, 0.03 


1460-5040 


20 


180 





24 3 0.02 


0.02 


1800-3580 


20 


360 


0,16,32,48 


24 3 0.03 


0.03 


1260-3040 


20 


360 


0,16,32,48 



TABLE I. Ensembles of QCD gauge field configurations generated by the RBC and UKQCD 
collaborations [71 1171 [3"T] for j3 = 2.13 with the Iwasaki gauge action that were used in this work. 
A is the separation between measurements in molecular dynamics time units. iV meas denotes the 
total number of measurements, and i src is the Euclidean time-slice location of the source. 
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IV. RESULTS 

In this section we present our results, focusing on the 24 3 ensembles, for the electromag- 
netic pseudo-scalar mass splittings (AM 2 ), EM LEC's in SU(3) and SU(2) chiral perturba- 
tion theory describing the pseudo-scalar masses, and the quark masses. Results from the 
16 3 ensemble are used for estimating finite volume effects which are discussed extensively in 



Sec. V Before turning to the results for AM 2 , we first describe lattice-artifact electromag- 
netic effects induced by the finite size of the fifth dimension of the DWF used to simulate 
the four dimensional u, d, and s quarks. 

In the following the notation uu (dd) denotes a meson whose two-point correlation func- 
tion is made from just the connected quark diagram using degenerate light quarks with equal 
charges, q = 2/3 (—1/3). Such a meson is neutral, but should not be confused with the 7r°, 
which requires disconnected quark diagrams. 

A. Electromagnetic effects in m res 



We first calculate the residual mass m res [32H31] from the pure QCD configurations. 
Then we consider the residual mass from the combined QCD+QED configurations so that 
the QED contribution to m rcs can be extracted. In the lattice DWF, m res is determined 
from the ratio 

(£w*K(o) 

R(t) = \* f-, (25) 

^j 5 a (*,tK(o) 

where t is the Euclidean time, Jg is a pseudo-scalar density evaluated at the mid-point of 
the extra dimension, 7r a denotes the usual 4d pseudo-scalar density, and the superscript a 



is a non-singlet flavor index. The correlation functions in Eq. (25) are computed from wall 
source, point sink, quark propagators. 

The residual mass is an ultra-violet, additive shift of the input, bare quark mass. Be- 
cause we are interested only in the EM meson mass-squared splittings, the leading order 
dependence of m res on the bare quark mass cancels, and we use a mass-independent residual 
mass in our later analysis that can be identified by extrapolating R(t) for the unitary quark 
masses to mj = with a suitable t-average. 
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16 3 24 3 


TYli 




chiral limit 


0.003148(46) 0.003203(15) 


0.005 


N/A 0.003222(16) 


0.01 


0.003177(31) 0.003230(15) 


0.02 


0.003262(29) 0.003261(16) 


0.03 


0.003267(28) 0.003297(15) 



TABLE II. The QCD residual mass for 16 3 and 24 3 lattice sizes. The data correspond to unitary 



mass points. The fit range for R{t) (defined in Eq. (25) is 9 < t < Nt/2 



Table [IT] shows the numerical result of the residual mass computed from the QCD con- 
figurations alone. In each case, R(t) was averaged over the range 9 < t < N t /2 for the size 
N t of the lattice in the time direction, after folding the correlation function about Nt/2. 
Figure [T] shows the chiral extrapolation of m res . The residual mass at the chiral limit is very 
close between the 16 3 and 24 3 lattices. m res is around 0.003, which is comparable to the 
lightest input sea quark mass, mi = 0.005, and larger than the smallest valence quark mass 
mf — 0.001, so the effect of the explicit violation of chiral symmetry from finite L s is not 
negligible in our calculation. Our measured values are roughly consistent with those found 
by the RBC and UKQCD collaborations [2 EE]. 

Next we consider the QED contribution to the residual quark mass. The QED contribu- 
tion from quark flavor % can be expressed as 

m reSii (QCD + QED) - m rcs (QCD) = e 2 C 2 ql (26) 

where m reSi j(QCD + QED) means the residual mass computed on the combined QCD+QED 
configurations and m res (QCD), the residual mass computed on the pure QCD configurations. 
Both are evaluated at m/ = and the former with physical quark charge qi. C2, which is of 
order 0(m res ), parametrizes the QED contribution to the additive shift of the quark mass. 
Although we compute this correction via the Ward-Takahashi identity for DWF [32], using 



a neutral meson made with degenerate, equally charged quarks, the form of Eq. (26) is 
completely consistent with a calculation in weak-coupling perturbation theory, say from the 
one-loop self-energy Feynman diagram for a quark with charge qi. In our chiral perturbation 
theory power counting, the QED contribution to the residual mass is 0(a em m res ) and must 
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0.0035 



0.0034 



0.0033 



0.0032 



0.0031 



0.003 




0.01 0.02 0.03 0.04 



m, 



0.0035 



0.0034 



0.0033 



0.0032 



0.0031 



0.003 




0.01 0.02 0.03 0.04 



m, 



FIG. 1. The QCD residual mass for 16 3 (upper) and 24 3 (lower) lattice sizes. The data correspond 
to unitary mass points. The linear chiral extrapolation to the m/ = limit is also shown on the 
plot. 
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L s uu dd 
16 3 lattice size 

16 2.597(23) 2.532(22) 

32 0.309(16) 0.301(16) 
24 3 lattice size 

16 2.585(7) 2.519(7) 



TABLE III. C2 (xlO ) in Eq. (26), representing the electromagnetic contribution to the residual 



mass, uu and dd denote the type of correlation function used to extract C2. 
therefore be included in our NLO analysis discussed in the next section. 



To compute the residual mass and extract the EM contribution via Eq. (26), we use uu 
or dd correlation functions in Eq. ( J25j ) |^] . The total contribution to the meson mass-squared 
due to explicit chiral symmetry breaking is, as in the case of pure DWF QCD, just the 
sum of contributions from each quark in the meson, modulo higher order than 0(a em m res )- 



corrections. Table III shows the results for C2 from uu and dd correlation functions. They 
agree well up to two digits, which implies that the 0(otl m m res ) contribution is quite small. 
These differences are higher order in chiral perturbation theory relative to the one we work 
to in this paper, so we neglect them. We note that C2 e 2 qf is the expected size, 0(a em m Tes ). 
The attained statistical precision on C2, which is impressive, of course stems from the fact 
that m reSj j(QCD + QED) and m res (QCD) are computed on exactly the same set of gluon 
configurations, so they are highly correlated, and the QCD fluctuations cancel between 



them. In addition, C2 appears to be insensitive to the volume (see Tab. Ill), presumably 
because the residual mass arises from the UV, short distance, regime. 



B. Meson mass splittings 

The electromagnetic mass splittings are determined from the pseudo-scalar masses com- 
puted with e / and e = 0, using the same gluon configurations. We use the additional 
trick of averaging correlation functions over ±e, configuration-by-configuration [2| 135]. 
3 In an earlier paper [2] we mistakenly included an independent contribution, proportional to qiqj, to the 
residual mass for the charged mesons made of quarks with charges qi and qj. This is clearly inconsistent 
with flavor conservation and the definition of a renormalized quark mass defined in perturbation theory. 
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In Fig. [2j the improvement due to the ±e averaging is demonstrated for the meson mass- 
squared splitting. The vertical axis shows the ratio of the statistical error without the trick to 
that with the trick, so that larger values indicate smaller statistical error for the ±e averaging 
trick. In most cases there is a large decrease (~ 1/10) in the error over the naive factor of 
\/2 that would result simply from doubling of the measurements (dashed line), while the few 
points with ratio exactly equal to one correspond to combinations that are trivially invariant 
under the change e — > — e. This procedure corresponds to including the QED configuration 
— A Xtlt for each A Xtfl in the path integral and can exactly cancel unphysical 0(e) noise with 
finite statistics which would have obscured the physical 0(e 2 ) signal of interest, only the 
latter of which survives in the infinite statistics limit. Together, the complete procedure 
yields mass-squared splittings with sub-percent statistical precision. 




12 3 4 

|( qi -q 3 )/(e/3)| 



FIG. 2. A comparison of the statistical errors for the meson mass-squared splitting with and 
without the ±e averaging trick [2 |3S]. The vertical axis shows the ratio of the error without the 
average to that with the average, so that larger values indicate smaller statistical error from the ±e 
averaging trick. In most cases there is a large decrease (~ 1/10) in the error over the naive factor 
of y/2 that would result simply from doubling of the measurements (dashed line) . The few points 
with ratio exactly equal to one correspond to combinations that are trivially invariant under the 
change e — > — e, i.e., vn\ = 771,3 and q± = —q^. 
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The pseudo-scalar meson masses are obtained from single state fits to wall source, point 
sink correlation functions with periodic boundary conditions in time with use of the fit 
function 

C m (t - t STC ) = A[e -M(t-t src+ N t )%N t + e -M(N t -t + t src )%N t ^ (2?) 

where M is the ground state meson mass, and t src is the time slice where the source is placed. 
To improve statistics in some cases, we average results from two sources (see Table [T]). The 
fitting procedure is done with the standard x 2 minimization (maximum likelihood), and the 
error on the mass is obtained by the standard jackknife method. Since the meson correlation 
function is symmetric about the midpoint (from the source) in the time direction, we fold the 
data about this point and fit with a time range smaller than N t /2. Based on the obtained 
effective masses (a representative example is shown in Fig. [3]), for all correlation functions 
we chose a fit range of 9 < t — t s < Nt/2. 

0.3 i 1 1 , , 1 1 — i 




0.1 1 1 1 1 1 1 1 — 1 

10 20 30 40 50 60 

t 

FIG. 3. Representative effective masses. Lattice size 24 3 . m sea =0.005, mi=m3=0.01, q\ = 1/3 
and t/3 = (upper points) and m sea =mi=m3=0.005, q± = 1/3 and #3 = —1/3 (lower points). The 
horizontal lines represent the fit result. 



The pseudo-scalar meson masses are tabulated in Tabs. XIII - XIV We have extracted 
the masses in two ways, one being from the fits to the correlation functions using the full 
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covariance matrix and the other being uncorrelated fits following [7J [T7]. The values of 
X 2 /dof for the covariant fits are roughly one for the 24 3 ensembles, but somewhat higher 
in some cases for the 16 3 ones and for the heavier quark masses on both ensembles. Such 
behavior for the 16 3 ensembles was seen in the earlier, pure QCD, analysis using these 
configurations, and was attributed to an inferior gauge field evolution algorithm 1311 . An 



improved algorithm was used to generate the 24 3 ensemble. From Tabs. XIII - XIV the 
masses and errors determined with either fit method agree quite well. Our final analysis is 
based on the masses from the uncorrelated fits in order to be consistent with the analysis 
in Refs. [H [T7| from which we take the pure QCD LEC's. The typical statistical error on 
the mass is at the half of a percent level and smaller. 

The meson mass-squared splittings are given by AM 2 = M 2 (e 7^ 0) — M 2 (e = 0), and 
the errors are again computed using a jackknife procedure. As an example, in Figs. [4] 
and M AM 2 for the dd meson is shown. Only the unitary points appear in the figure. A 



full summary of the mass-squared splittings is given in Tabs. |XV| and |XVI| The promised 
statistical precision is observed. Even though the errors on the masses themselves are of 
the same order as the mass difference, the splitting is statistically well resolved under the 
jackknife analysis thanks to the strong statistical correlation between e = and e ^ 0. 

We pause to compare the observed explicit chiral symmetry breaking effects to those 
expected from the discussion of the residual mass in the previous section. In the chiral limit, 
ui j = — m res (QCD), and in the absence of EM induced explicit chiral symmetry breaking 
(L s — > 00), the neutral meson mass-squared should vanish (up to a 2 m corrections which we 
ignore), and so too should the splittings. But it is clear from Figs. [4] and [5] that the dd mass- 
squared splitting does not (the same is true for the uu meson). Following the discussion 



in Sec. IV A and from the result of the pseudo-scalar mass-squared at the lowest-order in 
chiral perturbation theory, the shift in the splitting in the chiral limit should be 2i?oC*2e 2 gJ 
or 25C*2e 2 gJ, depending on whether we choose SU(3) or SU(2) chiral perturbation theory. 
A simple linear fit, also shown in Fig. |4| suggests this is true. Note that at NLO there are no 
logs in the splitting of neutral mesons made from only connected quark propagators, that is, 
a "charged" meson whose net charge happens to be zero. Further, by making L s larger, this 
lattice artifact should be (exponentially) reduced, which is also clear from the Fig. [4] where 
for L s = 32 the shift has been reduced by roughly a factor of ten, and the splitting nearly 
vanishes. Similar results hold for the uu mesons. The result based on the Ward-Takahashi 
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< 




FIG. 4. AM 2 for the dd meson, with L s = 16 (upper set of line and plots) and 32, 16 3 lattice 
size from the SU(3) fit. The extrapolated values (box) are e 2 5 mTes {li + qf). For comparison, we 
also show the values of BoC2e 2 (q 2 + q 2 ) obtained from the Ward-Takahashi Identity (values are 
slightly shifted horizontally to the left for clarity). 5 mrcs is obtained from the fit range of 0.01-0.02 
for L s = 16 and 0.01-0.03 for L s = 32. The error on BoC2e 2 (q 2 + g|) comes mostly from the error 
on Bq. 



Identity depends also on the value of B or B, depending on whether we choose SU(3) or 
SU(2) chiral perturbation theory, which introduces some uncertainty. On the other hand, 
Figs. [4] and [5] clearly show this effect is due to finite L s chiral symmetry breaking, and that it 
can be precisely subtracted from the physical splitting by introducing a new lattice-artifact 
LEC to the fit, e 2 5 m[CS (q 2 + q 2 ). We conclude that the explicit chiral symmetry breaking 
artifacts induced by finite L s and QED interactions are precisely quantifiable at NLO in 
chiral perturbation theory and that higher order terms can be safely neglected, so these 
artifacts can be robustly eliminated, just as in the case for pure (DWF) QCD. 
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0.00025 



0.0002 



0.00015 



0.0001 




0.01 0.02 0.03 

m f +m res 



FIG. 5. Same as Fig. [4] but for lattice size 24 3 and L s = 16. 



1. infinite volume fits 



In Figs. [6] and [7[ the meson mass-squared splittings are shown for the unitary quark 
mass points, for both 16 3 and 24 3 ensembles, respectively. For now, we concentrate on 
the 24 3 ensemble, and fit the mass-squared splittings to the infinite volume, NLO, chiral 
perturbation theory formulas described in Sec. [XTJ The formulas require the values for the 
pure QCD LEC's, some of which we have not computed. The pure QCD LEC's, including 
Fq, F, B , B, and the Gasser-Leutwyler L's, have been calculated already by the RBC and 
UKQCD collaborations from a larger ensemble of which the present one is a subset. We use 
these values in our fits, in a combined super-jackknife analysis so that the statistical errors 
on the QCD parameters are fed into our analysis. 

Figure [7] shows the fit to the full SU(3)l x SU(3)r NLO formula, which is summarized in 
Tables IV and[VJ The quark mass range in the fit is m 1) m 3 < 0.01, and the x 2 /dof for these 
uncorrelated fits is about two. x 2 degrades significantly if larger quark mass points are used 
in the fits. Only unitary points are shown in the figure for clarity while all of the (allowed) 
quark mass and charge combinations for the mesons have been used in the analysis. For 
the 24 3 ensemble, this amounts to 52 data points for mi, m 3 < 0.01. The charged meson 
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splittings should not vanish in the chiral limit, rrif = — m res ; this is just the LO Dashen 
term proportional to a cm and the lattice-artifact chiral symmetry breaking. The neutral 
meson splittings do not vanish either due to the latter. The chiral logarithms reduce the LO 
Dashen term relative to the value given by a simple linear ansatz. Recall that the splittings 
of "neutral" mesons made from connected quark diagrams only do not contain logs at NLO, 
so their chiral behavior is particularly simple. 

Figures [7] and [8] show similar fits for SU(2) L x SU(2) R -plus-kaon chiral perturbation the- 
ory for the pions and kaons, respectively. Here we use the same range for the light quark 
masses, and for the kaons the valence strange quark is fixed to either 0.02 or 0.03. x 2 /dof is 
similar to the SU(3) case for the pion and also significantly degrades when the quark mass 
range is extended upwards. For the kaon fits y 2 is small. The total number of data points 
in the fits are 52 for the pions and 36 for the kaons. The SU(2) LEC's are also summarized 
in Tables ED -ED 

Before proceeding, we address a subtly in the kaon fits that was not recognized until after 
the correlation functions had already been computed. Our original plan was to use an SU(3) 
chiral perturbation theory analysis only, for quark masses in the range 0.005 — 0.03, and 
non-degenerate meson correlation functions were computed for these masses in all possible 
combinations. However, learning first from the pure QCD analysis [7] [15], [16], and later 
from our own, it became clear that 0.02 and 0.03 were too heavy, and that SU(2) chiral 
perturbation theory would be needed to access the physical strange quark mass regime. We 
decided to include a lighter valence mass point, 0.001, to augment our fits, but since this was 
a new, separate calculation, only the mass-degenerate mesons could be computed. Thus, in 
our kaon fits, we have only two valence and two sea quark mass combinations available for 
the region m u ^ < 0.01. Now the subtly: it turns out these combinations of quark masses 



and charges are not enough to constrain all 10 LEC's appearing in Eq. (23). There is one 



direction in the multi-dimensional parameter space that is not linearly-independent from the 

if or 4 s 



rest. Fixing any one of the LEC's to zero, except A^' 1 ' 1 ' or A^' 2 \ results in a stable fit with 



the same x 2 , but with different values of the LEC's. While these fits all agree exactly when 
evaluated at the data points used in the fits, they differ elsewhere. There are two ways to fix 
this problem of an accidental flat direction in the \ 2 function at our disposal. First, keeping 
the same quark mass range, use the technique of singular-value-decomposition (36] (SVD) 
to determine all 10 LEC's. Second, increase the number of sea or valence quark mass points 
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in the fit, so the parameter directions are all linearly independent. While treating the (next 
available) mass 0.02 quark simultaneously as light and strange contradicts our assumption 
that mi/m s <C 1, nevertheless it allows the LEC's to be linearly independent, and only 
slightly increases x 2 which is still small. In practice, we only added the 0.02 valence quark 
mass to the kaon fit, keeping the light sea quark mass < 0.01. As it happens, the quark 
masses determined from these two methods agree well, giving confidence that the SVD fit 
procedure, which we use for our central values, is reliable. Further, in the case where 0.02 
data points were used, setting each of the LEC's to zero in turn resulted in much bigger \ 2 
values except for Xg — Xg , the ones related to sea quark masses which are not constrained 
as well, = gave the smallest % 2 . In each of these cases the quark masses agreed 

within statistical errors to the full SVD fit. We use the difference in the central values of the 
quark masses from the two procedures as an estimate of one of the systematic errors due to 
fitting. 



From Table |VJ we can see a large effect on C going from SU(3), where it is almost 
zero, to SU(2) where it is almost ten times larger. Recall that in the SU(2) theory the 
contributions of the strange quark terms in the SU(3) theory are absorbed into (see 



Eq. (19)). This situation is reminiscent of the pion decay constant in pure QCD computed 
on these lattices and for the same range of quark masses; the logs in that case also tend to 
significantly reduce the LO contribution over a simple analytic function, and the physical 
value [7J [15J [16]. Here, especially in the SU(3) case, the effect is even more dramatic. The 
other pion electromagnetic LEC's are roughly the same in both theories. In the SU(2) 
case, the size of the NLO EM correction turns out to be smaller than the LO one, showing 
compatibility with the chiral expansion. Finally, in Table [V] we show LEC's corresponding 
to the phenomenological parameter set presented in Ref. [UJ. The fact that the SU(3) 
NLO LEC's computed here (left-most column) do not agree is not surprising since the LO 
LEC, C, is clearly underestimated by a large degree. Note that to compare values of C, a 
factor of a -4 needs to be introduced, as well. We discuss the Dashen term further in Sec. fV 
after presenting the finite volume fits. 
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SU(3) inf.v SU(2) inf. a 



W°B 


2.15(11) 


2.348(44) 


10 2 F 


3.43(19) 


4.55(10) 


10 6 (2L 6 - L 4 ) 


-2.6(29.6) 


2.9(45.3) 


10 4 (2L 8 - L s ) 


5.42(29) 


4.36(31) 


10 5 L 4 


1.7(5.5) 


2.48(89) 


10 4 L 5 


2.02(63) 


5.49(47) 


10 3 m res 


3.131(27) 


3.131(27) 


a-^GeV) 


1.784(44) 


1.784(44) 



TABLE IV. The QCD LEC's from RBC/UKQCD collaboration's infinite volume fits on 24 3 lattices 
with SU(3) and SU(2) PQxPT |17j . They were computed from a larger ensemble of lattices than 
used in [7j. All of the QCD LEC's are defined at the chiral scale A x = 1 GeV. The labels in the 
first column correspond to SU(3) definitions; the analogous LEC for SU(2) is given in the third 
column. 

2. finite volume fits 

Next we include in our fits the finite volume corrections to the chiral logarithms using 



Eq. (12) for the SU(3) fit, and Eq. (22) for the pion and the results in Appendix B2 for the 
kaon in the SU(2) fit. We continue to use the pure QCD, infinite volume, LEC's from |17j . 
Since the finite volume effects in QCD are very small compared to the QED ones, we ignore 
the former. Figure [9] shows the modified fits for the pions and kaons on the 24 3 lattice. The 
LEC's are given in Table [V] and |VI[ The largest change by far is in C, the LO Dashen term, 
which roughly doubles in the SU(2) case and increases by a factor of four in the SU(3) case. 
Note, it is still much larger for the SU(2) fit. This is consistent with the observed large effect 
in the charged meson splitting compared to the neutral. Fortunately, this huge change does 
not greatly affect the values of the quark masses, as we shall see. 



C. quark masses 



Having determined the LEC's to NLO describing the pseudo-scalar masses in chiral per- 
turbation theory, we now turn to fixing the physical quark masses at the (arbitrary) low 
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SU(3) SU(2) 





inf.v. 


f.v. 


inf.v. 


f.v. 


S T I ( 3 1 -1- n h p n o m 


10 7 C 


2.2(2.0) 


9.3(2.4) 


18 3(1 8) 


32 9(2 3) 


41 


1 D 2 Vn 


± .UOl 1U 1 






i.ou 1 \ j 


n 1 q 


LKJ 1 3 


— 1 1 85C741 

± ± .UU Nil 


—5 17(70) 


— 1 61 Cfi2l 


—4 00(62") 


1.25 


io 3 y 4 


13.4(1.7) 


9.7(1.7) 


10.6(1.1) 


8.3(1.1) 
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Dashen's term(MeV) 


0.40(37) 


1.68(39) 


1.88(18) 


3.38(23) 


3.7 


X 2 /dof 


2.11(73) 


2.12(78) 


2.27(78) 


2.19(82) 





TABLE V. The SU(3) PQxPT and SU(2) pion PQxPT QED LEC's from fits of the mass-squared 
splittings measured on the 24 3 lattices. All of the LEC's are defined at chiral scale A x = 1 GeV 
and are given in lattice units. The quark mass range in the fits is mj^ <= 0.01. "inf.v." and 
"f.v." means infinite and finite volume fits, respectively. "SU(3)+phenom." refers to a parameter 
set presented based on phenomenology and using SU(3) xPT [II]- Labels in the first column 
correspond to SU(3) definitions. "Dashen's term" is the LO result for the mass splitting in the 
chiral limit. 

energy scale of 2 GeV. First, the bare quark masses are determined by solving Eq. Q or 



Eqs. (20) and (23) evaluated at the physical meson masses [I] (in MeV) 



M w ± = 139.57018 ± 0.00035 (28) 
M K o = 497.614 ± 0.024 (29) 
M K ± = 493.667 ± 0.016, (30) 

where only the central values are used in our analysis since the errors are negligible compared 
to the lattice results. Using a -1 = 1.784(44) GeV and the pure QCD nonperturbative 
renormalization constant Z m = 1.546 (2) (43) [T7| computed by the RBC and UKQCD 



collaborations, MS light quark masses are given in Table VII for infinite volume, finite 
volume, SU(3), and SU(2) fits. We have not included the 0(a em ) renormalization of the 
quark mass from QED interactions. These are similar to those found in our earlier two flavor 
work [2] , also using DWF, but which used a more crude chiral perturbation theory analysis 
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10 2 M 2 4.804(88) 6.89(10) 7.37(36) 
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- 
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0.376(42) 
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-8.73(86) 
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X 2 /dof 0.4578(52) 0.2869(40) - 0.4578(52) 0.2869(40) 



TABLE VI. Kaon QCD and QED LEC's extracted from 24 3 lattice size data. LEC's are in lattice 
units. The kaon is composed by one light- (mi) and one strange- (m.3) quark. We choose mi < 0.01 
and m,3=0.02 or 0.03. The light sea quark is chosen as m sea < 0.01. The mass of the strange sea 
quark is fixed at 0.04. The kaon QCD LEC's are quoted from RBC/UKQCD's work [17J. x 2 /dof 
refers to the fit using the SVD method |36j . 

that did not include logarithms. The strange quark mass is somewhat lower here, which may 
be a real flavor-dependent effect [TJ [2]. We also note that in the combined continuum limit 
analysis mentioned earlier, the RBC and UKQCD collaborations find that the strange quark 
mass is even smaller [T5HT7"] . The average light quark mass is close to the value determined 
in pure QCD [H [15l HE] . 
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FIG. 6. Meson mass-squared splittings. 16 3 lattice size. Infinite volume linear fit (upper panel) 
and infinite volume SU(3) chiral log fit (lower panel). The fit range of the linear fit is 0.01-0.03. Fit 
ranges of chiral log fits, 0.01-0.03 (solid line) and 0.01-0.02 (dashed line). Data points correspond to 
ud, uu and dd mesons, respectively, from top to bottom. Only unitary points are shown, although 
all of the partially quenched points were used in the fit. 
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FIG. 7. Meson mass-squared splittings. 24 3 lattice size. Infinite volume linear fit (upper panel), 
and infinite volume SU(3) and SU(2) chiral log fits (lower panel). The fit range of the linear fit 
is 0.005-0.03. Fit range of chiral log fits is 0.005-0.01. The solid (dashed) line in the lower panel 
represents the SU(3) (SU(2)) fit. Data points correspond to ud,uu and dd mesons, respectively, 
from top to bottom. Only unitary points are shown, although all of the partially quenched points 
were used in the fit. 
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FIG. 8. Kaon mass-squared splitting and infinite volume SU(2) kaon fit. The mass of the strange 
valence quark is fixed at 0.03, and m sca = 0.005. Different lines in the plot correspond to different 
charge combinations of the valence quarks. 



SU(3) SU(2) 





inf.v 


f.v 


inf.v. 


f.v. 


m u 


2.606(89) 


2.318(91) 


2.54(10) 


2.24(10) 


m d 


4.50(16) 


4.60(16) 


4.53(15) 


4.65(15) 


m s 


89.1(3.6) 


89.1(3.6) 


97.7(2.9) 


97.6(2.9) 


m d - m u 


1.900(99) 


2.28(11) 


1.993(67) 


2.411(65) 


m ud 


3.55(12) 


3.46(12) 


3.54(12) 


3.44(12) 


m u /m d 


0.578(11) 


0.503(12) 0.5608(87) 0.4818(96) 


m s /m ud 25.07(36) 


25.73(36) 


27.58(27) 


28.31(29) 



TABLE VII. The u, d and s quark masses determined from QCD+QED interaction on 24 3 lattices. 
The values are given in MeV and the MS scheme at renormalization scale fi = 2 GeV. SU(3) or 
SU(2) mean quark masses from SU(3) PQxPT or SU(2) PQxPT + kaon theory. 



34 



0.0016 r — — i 1— 

0.0014 - 

0.0012 - ^ ^^^^^ 
0.001 - 
2 0.0008 - ,3/""^ 
0.0006 - / 
0.0004 / 
0.0002 i 

o I 1 1 1 

0.01 0.02 0.03 0.04 

m f +m res 



0.0011 




0.0007 1 1 1 1 1 

0.005 0.01 0.015 0.02 

m 1+ m res 

FIG. 9. 24 3 SU(2) chiral log infinite volume and finite volume fits for pion (upper) and kaon 
(lower) mass-squared splittings. Lines correspond to fit results. The fit range is 0.005-0.01. The 
solid (dashed) line represents the infinite (finite) volume fit. In the upper panel, the fit curves are 
evaluated for degenerate unitary light quarks. For the lower panel, the curves are evaluated for 
msea = 0.005 and m% = 0.03. Data points in the plot correspond to qi = 2/3e and qs = — l/3e, 
but all partially quenched points allowed by the fit range were used in the fit. 
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V. SYSTEMATIC ERRORS 

In this section we examine the important systematic errors in our calculation: the chiral 
extrapolations, finite volume, non-zero lattice spacing, and QED quenching. In each case 
we estimate the size of the effect on the values of the quark masses and investigate the effect 
on the LO electromagnetic LEC's. Similar systematic uncertainties have been given for the 
pure QCD sector [3 HS1 US]- 

To estimate the systematic errors, the change in a quantity is computed under the in- 
fluence of a change in how that quantity is computed, for example, by using a different fit 
formula. Since the data is the same, or there is significant overlap, in each case, we compute 
the change under the (super-)jackknife procedure in order to assess its significance. Central 
values of all quantities are quoted for the finite volume, SU(2) chiral perturbation theory 
fits which we believe give the most accurate results. The systematic errors computed in the 
following come from comparison to these central values. 

A. chiral extrapolations 

Previous studies have used the difference in analytic and chiral perturbation theory fits 
to estimate the chiral extrapolation error that stems from using unphysical heavy quarks [7J 
[III EDS, [T7J [37J EE] • One can also estimate the error in chiral perturbation theory alone by 
comparing the relative sizes of LO, NLO, or even NNLO corrections to a given quantity. For 
the latter to work, the estimates of the higher order contributions must be accurate. 

It is perhaps not surprising to find that the meson mass-squared splittings show little 
trace of the chiral logarithms. For the mass range of pions in this study, it is well known 
that low energy observables like the meson mass-squared or decay constant exhibit more 
or less linear dependence on the quark mass. In Fig. [7j the charged pseudo-scalar splitting 
appears linear over the range of unitary points shown in the figure. Nevertheless, the fits 
to our data do show that NLO chiral perturbation theory (chiral logs) is consistent with 
the data. A similar conclusion was reached in the pure QCD case [74 IT5T4T7] . To NLO in 
chiral perturbation theory, there are no logs for the neutral mesons made from connected 
quark propagators like those studied here, and indeed the neutral splittings, too, appear to 
be quite linear. 
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We do point out one aspect of the EM logs that leads one to expect a noticeable affect. 
They behave like am log m, not m 2 logm as the pure QCD logs do. A factor of a has 
replaced a factor of the quark mass. In fact they are like the quenched logs in pure QCD in 
this respect. 

The first step in estimating the systematic error is to determine the fit range, or range 
of quark masses included in the fit. The available ranges are summarized in Table |TJ In 
Refs. [7J rHHTT] it was shown that for the same ensembles used in this work, SU(3) and 
SU(2) chiral perturbation theories give sensible fit results for pion masses less than about 400 
MeV, or for bare quark masses satisfying m/ < 0.01. It is possible that the range is different, 
perhaps larger, for the EM splittings. After all, most of the pure QCD contributions at LO 
and NLO completely cancel in the EM splittings (some of the pure QCD LEC's survive at 
0(a em m)). We work with uncorrelated fits, though our data are highly correlated, because 
there are too many mass and charge combinations to accurately determine the correlations 
on this finite statistical ensemble. The uncorrelated fits have been shown to agree with 
correlated ones when the covariance matrix is well determined, and when it is not, the 
correlated fits break down [TTl 139] . As already mentioned, when the quark mass range is 
extended upwards, for both SU(3) and SU(2) (pion) fits, x 2 /dof increases noticeably, by 
more than a factor of two. Since we use uncorrelated fits, this \ 2 is n °t an absolute test of 
goodness of fit, though we expect changes do indicate relative goodness of fit. Thus we stick 
with the range m-y, m% < 0.01 for the light quarks to quote central values and to estimate 
systematic errors. One of the systematic errors is the difference in the central values for the 
quark masses determined from this restricted range and those values computed from the 
range mi, m 3 < 0.02 for the light quarks, mi, m 3 < 0.01 corresponds to valence pions in the 
range 250-420 MeV. 

For the mass range m\,m^ < 0.01, the most important change, which is anticipated 
in Fig. [7j is that the Dashen term increases significantly when the logs are omitted from 
the SU(3) fit. C increases by about a factor of five, although it is still small compared to 
the value one would obtain from the physical splitting. For SU(2), the situation is much 
different; C changes very little, about two percent. Presumably, the large logs containing 
the strange quark mass contribute to the LO term in this case, and the remaining effect of 
the light logs is not as important. The higher order terms change more, but without logs to 
affect their running, there is not much sense in comparing the changes. 
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In Fig. 10 the LO and NLO contributions in finite- volume SU(2)l X SU(2)r chiral pertur- 
bation theory for the charged pion mass splittings are shown. At the values of quark masses 
used in our calculation, after accounting for the S mrcs contribution, the NLO contributions 
are about 50-100 % of the LO contribution. It is interesting to see how — m u is affected 
at the various orders in chiral perturbation theory. Using the LEC's determined in the full 
SU(2)l x SU(2)R-plus-kaon fits, we find that the NLO contributions increase — m u by a 
bit less than 2%. 

Taking all the above uncertainties due to fitting into account, we estimate systematic 
errors of about four and zero percent for the up and down, and strange quark masses, 



respectively. These are collected in Tab. VIII 
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FIG. 10. The LO and NLO in finite volume SU(2) chiral perturbation theory contributions to the 
EM meson mass splitting. The dashed line corresponds to the finite volume fit result. The data 
points shown are for charged mesons with q\ = 2/3 and q% = —1/3. The lower horizontal line 
gives the contribution of the lattice artifact e 2 S rnics (qf + g|) while the upper horizontal line gives 
the sum of this contribution and Dashen's Term (in other words, their difference is just the LO 
contribution). The solid line corresponds to the total LO+NLO+e 2 ^^^ 2 + g|) contributions 
based on the fitted, finite volume LEC's, but evaluated with the infinite volume logarithms. The 
fit curves are evaluated for degenerate unitary light quarks. 
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B. finite volume 



The effect of finite voiume on the measured charged-meson splittings is large, as we have 



seen. In Fig. 11 the difference between the measured 16 3 and 24 3 EM splittings is about 
15-20%. The LO LEC C changes dramatically, by about a factor of two, when the finite 
volume corrections at NLO are included in the chiral perturbation theory fits (see Table [V] 



and Fig. 10 ). Besides the usual special functions that replace the infinite volume logs, a large, 
negative constant appears in the finite volume formula, —3Kq 2 3 /4:7rL 2 [10] with k ~ 2.837, 
which cancels against an enhanced value of C. 

To estimate how reliable the NLO finite volume corrections are, one can use the LEC's 
from the 24 3 fits to predict the finite volume shift in the 16 3 splitting. The fit and prediction 



are shown in Fig. 11. First, the SU(2) fit agrees well with the 24 3 results for mj < 0.01 
which is the quark mass range used in the fit. For larger masses the fit deviates significantly 
from the data and suggests that NLO chiral perturbation theory is not reliable for these 
masses. Even for rrif = 0.01, where we may trust NLO chiral perturbation theory, the 
theory over-predicts the shift on the 16 3 lattice by about a factor of two. The NLO LEC's 
Y 3 and Y 5 also have large finite volume shifts. 



From Table fVT] the shifts in the kaon mass-squared LEC's are much smaller, especially 
the ones representing the LO Dashen term {A^ ,l \ A^\ A^ l,l \ and A^' 2 ^). 

Even though some of the LEC's show large finite volume shifts, the ultimate shifts in 
the quark masses are smaller. The largest, about 14%, occurs for the up quark mass. The 
down quark mass is affected much less, about 3%, and the shift in the strange quark mass 
is negligible. 

From the pure QCD calculations, we know the finite volume effects in the 24 3 meson 
masses are at about the one percent level pT5HT7] . and therefore the QED finite volume 
corrections dominate. 



The finite volume errors on the quark masses are summarized in Table VIII 



C. non-zero lattice spacing 

Since our calculation has only been done at a single lattice spacing, we can not estimate 
the non-zero lattice spacing errors directly. However, by now there is much evidence that 
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FIG. 11. Finite volume effect in the measured EM splittings. All of the data points have q\ = 2/3 
and qs = —1/3. Circles and squares correspond to 24 3 and 16 3 lattice sizes, respectively. The solid 
line is from the finite volume fit on 24 3 ensembles. The dashed line is the theoretical prediction for 
16 3 lattices based on the LEC's extracted from 24 3 finite volume fit. The fit curves are evaluated 
for degenerate unitary light quarks. 

these 0(a 2 + m res a) discretization errors are small in pure DWF QCD, and they should 
largely cancel in the splittings. Even assuming they do not cancel, there is no reason to 
expect they are enhanced over the pure QCD case. In the first QCD calculation using the 
24 3 ensemble, it was estimated that scaling errors were at about the four percent level for 
low energy quantities like the pion decay constant and the kaon [7j. Since then, a new 
calculation at the same physical volume but smaller lattice spacing has shown this estimate 
was about right, or perhaps a bit conservative [15HI7]. Of course, here we are interested 
only in the mass splittings. The pion and kaon masses are fixed to their continuum values, 
so they have no scaling errors. Instead, the lattice spacing errors enter in the LEC's and 
the physical quark masses. Therefore we assign a robust four percent scaling error to the 
quark masses, which will be eliminated in up-coming calculations on the finer lattice spacing 
ensemble [ToT - TP?] . This error also encompasses the uncertainty in setting the lattice scale 
itself, which as mentioned earlier differs by about 2 ~ 3 percent from the scale given in 
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Ref. jZj. 

The non-zero lattice spacing errors on the quark masses are summarized in Table VIII 



D. QED quenching 

As mentioned our calculation is done in quenched QED where the sea quarks are neutral. 
In chiral perturbation theory, we have neglected terms of order O(aem^sea), including logs. 
From a weak-coupling perturbation theory perspective in QCD+QED, we have neglected 
vacuum polarization effects at order 0(a em a s ). For the pions, the consequence is that the 
single (linear combination) LEC Y\ can not be determined. For the kaons there are several 



LEC's that can not be determined (see Eq. (23)). However, we do note that sea quark charge 



effects from the logs can be included a posteriori in our determination of the quark masses. 

Since the LEC's absorb changes of scale in the logs, one way to estimate the effect of the 
missing LEC's, or counter-terms, is to mark the change in the quark masses when these logs 
are included, or not. This leads to a negligible change in the quark masses. From Table |Vj 
the other EM LEC's have magnitudes roughly in the range 0.01 to 0.001. Setting Y± at 
the high end, Y\ = ±0.01, the quark masses again change very little. Of course, the LEC's 
calculated with g soa = will differ from those with g sea ^ 0, by 0(az em ). This is higher 
order for all the LEC's determined here except C for the pions and A^\ A^\ A^' 1 ' , 

(s 2) 

and A K ' for the kaons. Taking all of the above into account, we quote a conservative two 
percent systematic error in our quark mass determination, stemming from the quenched 
approximation to QED. 

Of course the above is only a rough estimate, so presently we are investigating the use 
of so-called re- weighting techniques to eliminate the quenching effects [24T[2"6| l28| |29| HO] . 
Re- weighting is simply the use of ratio (s) of fermion determinants in observable averages in 
order to include the desired dynamical-quark effects. The calculation of a determinant which 
is non-local in the fields is quite expensive, so stochastic estimators must be used to make the 
calculation tractable. Re-weighting in the strange quark mass to the a posteriori determined 
physical value has proved quite useful and efficient in recent 2+1 flavor simulations [15- 
131271. 
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VI. QUARK MASSES 

We use finite volume SU(2) chiral perturbation theory and the light quark mass range 
m < 0.01 to obtain our central values of the physical quark masses. 

The physical strange quark mass is determined from the kaon mass-squared which is an 
implicit function of the bare sea and valence strange quark masses, through its LEC's which 
are calculated for fixed valence strange quark masses 0.02 and 0.03 and fixed sea strange 
quark mass 0.04. Assuming that the m s dependence is modest in this region, the physical 
kaon mass-squared is determined from a linear extrapolation in the valence strange quark 
mass. A similar procedure was carried out in |7] where three data points in the range 0.02- 
0.04 showed the kaon mass-squared is well approximated by a linear function (it turns out 
that the physical strange quark corresponds to about 0.035). Because we have only carried 
out calculations at a single strange sea quark mass value of 0.04, the kaon mass-squared can 
not be evaluated at the physical strange sea quark mass. This partial quenching leads to a 
small systematic error that was conservatively estimated in [7J to be 2% for m s which we 
adopt here. It is added in quadrature to the total systematic error for m s which appears 
below. The systematic error on the light quark masses is about 0.7% which is negligible 
compared to the other systematic errors, so we ignore it. 

The statistical errors come from fits underneath a standard jackknife analysis. The QCD 
LEC's come from an analysis of the extended RBC/UKQCD 24 3 ensembles; the results are 
consistent with those in Ref. [TJ. All of the fits and corresponding LEC's are analyzed under 
a super-jackknife analysis so that statistical errors on all quantities, from all ensembles, are 
included. The systematic errors assigned have been discussed in this section. The mass- 
independent quark mass renormalization factor is 

Z^(/! = 2GeV) = 1.546(2)(43), (31) 

which is obtained via nonperturbative renormalization using the RI/SM0M 7m scheme [TTl 
|4"1TI4"5] . The second error is systematic, including 0((fia) 2 ), which will be removed when 
we take the continuum limit in future work (the 0(a era ) QED correction to Z m is omitted). 
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The final values are 

(32) 
(33) 
(34) 
(35) 
(36) 
(37) 
(38) 

where the first error is statistical, and the second is a total systematic error, derived by 



m u 


= 2.24 ±0.10 ±0.34 MeV 


™d 


= 4.65 ±0.15 ±0.32 MeV 


m s 


— cm.u in z.y in o.o ivie v 


m d - m u 


= 2.411 ± 0.065 ± 0.476 MeV 


m ud 


= 3.44 ±0.12 ±0.22 MeV 


m u /m d 


= 0.4818 ± 0.0096 ± 0.0860 


m s /m ud 


= 28.31 ±0.29 ± 1.77, 



adding the individual errors summarized in Table VIII in quadrature. We remind the reader 
that these central values are obtained from our SU(2), finite volume fits on the 24 3 ensembles. 

We note that the up quark mass obtained here is different from zero by more than 
six standard deviations, which seems to rule out the m u = solution to the strong CP 
problem. However, there is an extensive literature concerning this scenario to which we 
refer the interested reader. For a discussion of extracting the up quark mass by using 
chiral perturbation theory, and its consequences, see [4"6T[4"8"] . The possibility of instanton 
effects additively shifting the up quark mass is discussed in many places [49— 52j . In [53J, 



renormalization scheme dependence of the renormalized quark mass was discussed in the 



context the isospin breaking. The effect vanishes in the (perturbative) MS scheme. At 
this point, there seem to be no common consensus on if there is any non-perturbative 
contribution, which is related to the aforementioned instanton effects, and how large it 
might be. Our results are potentially susceptible to this uncertainty, as are all other quark 
masses renormalized in a perturbative scheme. 
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TABLE VIII. Summary of quark mass systematic errors. Central values quoted from the finite vol- 
ume, SU(2), chiral perturbation theory fit. Masses given in MeV. The quark mass renormalization 
error comes from the nonperturbative QCD result [TTj plus a one percent error from QED, added 
in quadrature. Systematic errors are given as a percent (%). The algebraic sign of each change 
comes from the difference (quantity under change) — (central value). 





value (stat. error) 


fit fv lat. 


spacing 


QED quenching 


; m s quenching 


; renorm 


m u 


2.24(10) 


+4.02 +13.50 


4 


2 




2.8 


m d 


4.65(15) 


+3.55 -2.48 


4 


2 




2.8 


m s 


97.6(2.9) 


+0.23 +0.07 


4 


2 


2 


2.8 


m d - m u 


2.411(65) 


+7.77 -17.35 


4 


2 




2.8 


m ud 


3.44(12) 


+2.75 +2.71 


4 


2 




2.8 


m u /m d 


0.4818(96) 


+5.45 +16.40 


4 








m s /m ud 


28.31(29) 


+2.91 -2.56 


4 


2 


2 
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VII. MESON MASS SPLITTINGS 

In Tab. |V|we give the contribution to the charged pion mass splitting in the chiral limit, 
or Dashen's term. The physical splitting, given in Eq. Q, is 4.5936(5) MeV. The SU(3) fit 
gives a very small value, about half an MeV. The finite volume fit dramatically increases 
the value, but it is less than half the physical value. The SU(2) fit gives a bigger value 
still, and after including finite volume corrections, it gives the LO EM correction to the 
pion mass difference (m n ± — win ) qed = 3.38(23) MeV. Coincidentally, this is about the 
same value obtained from the linear fit, 3.22(25)MeV. The value of m 2 ± in the chiral limit, 
which comes from the LO EM correction is 929(64) MeV 2 and is similar to the values, 
using a sum rule and lattice-computed vector and axial-vector correlation functions in pure 
QCD, reported in [5U |33]. Our value for (m„-± — m 7r o)q ED is roughly consistent with, but 
two statistical standard deviations smaller than, the value from phenomenology and SU(3) 
chiral perturbation theory reported in [TT], 3.7 MeV. 

The above suggests that NLO contributions at the physical quark masses may be as large 
as 25% of the total pion mass difference, m n + — m n o. Away from the chiral limit, there 
are corrections to m^o that we have not computed in the lattice calculation (disconnected 
diagrams), nor in chiral perturbation theory (logs). However, we can estimate some of them 



by evaluating Eq. (20) for m u = m d = m ud , q± = q 3 = q u and averaging it with the case for 

Qi = Q3 = qd 



2 X 

M (q 1 , q 3 ; mi) = - {M 2 (m 1 , q%, m 1 , q x ) + M 2 (m l3 g 3 , m ls q 3 )} ■ (39) 



This form can be inferred for the 7r° made with degenerate light valence quarks (m 3 = m x ) 
in our current study in which only the connected valence quark diagram is computed and 
QED is quenched. We focus on the one-particle irreducible two-point function S^o^ 2 ) of 
7T°, and pick out the part depending on the valence EM charges induced from the connected 
diagram. E^o(p 2 ) can be divided into a pure QCD part S^ CD (p 2 ) and a QED correction 



45 



E^ EI V) at order e 2 , 

= V) + £?o E V) ' (40) 
^ CD (P 2 ) = (^) 2 tr ((r 3 ) 2 ) A QCD (p 2 ) = A QCD (p 2 ) , (41) 

£2fV) = (^) 2 tr (r 3 Qr 3 g) ^(p 2 ) + 2 x (^'tr ((r 3 ) 2 Q 2 ) £> 2 (p 2 ) 

= l(qI + al)D 1 (p 2 ) + ( q 2 1 + q l)D 2 ( P *). (42) 



where r a (a = 1, 2, 3) denote the Pauli matrices and Q = diag(gi, g 3 ). In Eq. (42), the 
first term originates from the Feynman diagram in which a virtual photon is exchanged 
between two valence quark lines, while a photon propagates on the same valence quark lines 
and induces the second term. Because the functions D 12 (p 2 ) are given by QCD dynamics 
weighted by the photon propagator, the self-energy S 7r +(p 2 ) of the charged pion is also 
expressed in terms of these functions 

^+( P 2 ) = £Q c V) + £?fV), (43) 
£?fV) =tr (r + Qr_Q) J D 1 (p 2 )+tr((r + r_ + r_r + )g 2 ) D 2 {p 2 ) 

= qiq3D 1 (p 2 ) + (qi + ql)D 2 (p 2 ). (44) 



From Eqs. (40)-(44), the charge dependence of m 2 and m 2 + = M 2 (m 1; g 1; m 1; q 3 ), to the 



order relevant to us, is found as 

M 2 (m 1; Ql , m u q 3 ) = K + q^Ft + (qj + q 2 ) F 2 , 

mlo = K + l - (ql + q 2 3 ) F 1 + {q\ + q 2 3 ) F 2 

= M 2 (q u q 3 ; mi ), (45) 

where K denotes the QCD part to NLO of chiral perturbation, and F 12 the 0(e 2 ) and 
0(e 2 mi) part. The chiral symmetry as well as QED gauge invariance should give ^i| mi=0 = 
—2 F 2 \ m =0 to reproduce the EM charge dependence (q± — q 3 ) 2 of the LO EM correction 



to rn^ + . Using Eq. (39) for m 2 , we find the LO + NLO EM pion mass difference at the 
physical point to be m n + — m^o = 4.50(23) MeV. Phenomenology predicts that a small part 
of the total NLO correction is due to m u -m d ^ 0, 0.17(3) MeV [56] and 0.32(20) MeV [57]. 

For the kaons, the pure EM mass difference is rrix+ — m^o = 1.87(10) MeV, while the 
contribution from m u — m d ^ is —5.840(96) MeV. Here, of course, the result includes all 
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NLO corrections, and LEC's from the finite volume fit are used. These values are obtained 
by taking the SU(2) formula for the kaon mass-squared Mj c (m 1 ,qi,m 3 ,q 3 ), Eq. (23) 



M 2 K (m u , 2/3, m„ -1/3) - M 2 K (m u , -1/3, m„ -1/3) = A^M 2 K + A^~ m ^M 2 K 

+0{e 2 {m u -m d )) (46) 

where the contributions to the mass-squared splitting are defined as 

A (EM) M 2_ = M 2 K ( mudj 2/3, m ud , -1/3) - M 2 K {m ud , -1/3, m ud , -1/3), (47) 
A (m u - md ) M 2 K = M 2,( mu; q, m s , 0) - M 2 K (m d , 0, m s , 0). (48) 

A (EM)/K-m d ) M ^( M ^ + Mr± ^ are quotec i a bove. So out of a physical mass-squared 
splitting (M K o) 2 - (M K ±) 2 = 3902.7 MeV 2 , about -47(2)% is A^M 2 K and +148(2)% is 

The breaking of Dashen's theorem can also be parametrized by AE [58J, 

AE = Mjcjmi, q u m 3} q 3 ) - M 2 K {m x , q 3 , m 3 , q 3 ) _ ^ 
M 2 (mi, q h m h q 3 ) - M 2 (m 1 , q 3 , m h q 3 ) 

where mi is the light quark mass and m 3 is the strange. M 2 (m 1; q 3 , mi, q 3 ) is used here to 

represent m 2 ; no significant change of AE is observed in our numerical study even when 



the average (39) is adopted for m 2 in place of M 2 (mi, q 3 , mi, q 3 ). In the SU(3) chiral limit 
AE = since the LO Dashen terms are the same in the numerator and denominator. If the 
strange quark mass is fixed to its physical value, then it does not vanish, and can be much 
larger than zero, even in the light quark chiral limit. Notice that AE vanishes trivially in 
both SU(2) and SU(3) theories when mi — > m 3 . 



We show AE for our data in Fig. 12 where the artifact S mies (q 2 + q 2 ) has been subtracted 
for each value of the meson mass-squared. In the upper panel, fit results are shown for SU(3). 
The fit evaluated at the simulated mass points does a reasonable job of reproducing the data, 
though as m 3 increases differences emerge. This is not surprising since only mi,m 3 < 0.01 
points were used in the fit, and including larger values yielded significantly poorer fits. More 
troublesome is the light quark extrapolation which yields a large value of AE at the physical 
point, which can be understood from two primary causes. First, the numerator is quite large 
since m 3 is evaluated at the physical strange quark mass, leading to a large 0(ctm) correction 
to the charged kaon mass-squared. Second, the denominator becomes quite small because 
the LO Dashen term is quite small in the SU(3) fit (compared to NLO terms). Both facts, 
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of course, signal a breakdown in SU(3) chiral perturbation theory which renders the SU(3) 
AE unreliable. As noted in [TT], the sea quark charge LEC's drop out of AE, and only 



known logarithms remain. Adding these to the (cyan) physical curve in Fig. 12 changes it 
only slightly. 



In the lower panel of Fig. 12 we show analogous results for the SU(2) fits. While the 
SU(2) fits are more reliable since the LO contribution is larger compared to NLO, the 



latter corrections are still large (recall Fig. 10). As expected, the fit agrees better with the 
data points for larger values of m 1; but the extrapolated value at the physical point and 
infinite volume is still much larger than the data points. We find in quenched QED that 
AE = 0.628(59) where the error is statistical only. This is much larger than the value 
reported in our previous two flavor paper [2] and not much smaller than phenomenology 
and SU(3) chiral perturbation theory [11]. The main difference is that here we use full NLO 
chiral perturbation theory with finite volume corrections while in [2] only simple analytic 
fits were used. To properly address these large corrections, one needs to simulate with larger 
volumes and smaller quark masses, a project that is now underway. 

Perhaps more useful for other pure QCD simulations are the "physical" values of and 
itlk in pure QCD deduced from our SU(2) fits with m u = m d = m ud : 

m (QCD) = 13 4. 98 (23) MeV, (50) 
m (QCD) = 494 521(58) MeV. (51) 

The small statistical errors result because the physical pion and kaon meson masses were 
used to determine the physical quark masses from our fit. 



Finally, based on the quark masses in Eqs. (32)- (34) and (36), we examine the ratio 
introduced in Ref. [59] . 

_ m d - m u 2m ud 

"^quark mass — . j V / 

m s - m ud m s + m ud 

which is equal to 

k = {M&, ~ M^ ± ) QCD Ml 

M 2 (m dl 0, m a , 0) - M 2 K (m u , 0, m 8 , 0) M 2 (m ud , 0, m ud , 0) 
M^(m s ,0,m ud ,0) - M 2 (m ud 0,m udj 0) M|(m s ,0,m^0) ' 
up to NNLO in SU(3) ChPT [56]. For SU(3) we obtain 

«quarkmass = 0.00201(3), (55) 
«meson = 0.00201(3), (56) 
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FIG. 12. Breaking of Dashen's theorem for the quenched QED case. The unphysical contribution 
Sm rcs {<li + <zi) nas been subtracted from the data. Data for two values of the strange quark, 0.02 
and 0.03, are shown. The curves correspond to the SU(3) fits (upper panel) and SU(2) fits (lower 
panel). The cyan bands denote the infinite volume extrapolations with one standard deviation 
statistical errors, using the LEC's extracted from the finite volume fits; the sea and strange quark 
masses are fixed at their physical values. 
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while for SU(2), 

K q uarkmass = 0.00176(4), (57) 
^meson = 0.00191(3), (58) 

where the errors are statistical only. For SU(3) the values are quite consistent with each 
other, while for SU(2) there is a small difference. In [HH], k extracted from rj — > 7r°7r + 7r~ 
decays is 0.0019(3) while the 0(p 6 ) analysis in [57] gives k = 0.00260 at m s /m U( i = 24. 
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VIII. ISOSPIN BREAKING EFFECTS ON THE KAON DECAY CONSTANT 

In our results the up quark mass is about 35% smaller than average of the up and down 
quark masses, m u d- In principle, this isospin breaking effect may cause visible effects on 
phenomenologically important quantities when they are measured with sufficient accuracy. 
As we saw in the previous section, a major part of the Kaon mass splitting comes from the 
quark mass difference, m u — m^. 

Here we examine isospin breaking effects on the Kaon decay constant, fx- By combining 
the experimental decay widths, T(K — > ^(7)) and T(7r — > z//i(7)), and f n and fx, one can 
extract the corresponding ratio of CKM matrix elements \§U\ ■ In the latest global analysis by 
the FlaviaNet Working group on Kaon Decays ED, iF 77^ is obtained from experimental 

fir Vud 

results with an accuracy of 0.2%. The ratio of the decay constants used are from their world 
average of lattice QCD simulations, and is 

?f = 1.193(5) [0.4%]. (59) 

Jit 

We address a question: how far does the value of fx shift when the light quark mass in 
the Kaon is changed from from m u d to m n ? Some lattice determinations of fx use m u d as 
the light quark mass while the experiments measure decays of the charged Kaon to obtain 
fx±, which is made of an up (and strange) quark. So it is relevant to know if the shift 
fx(n^ud) — fK( m u) is comparable in size to the total error on the ratio, 0.4%. We note 
the analyses of V us /V u d in JJITl EI] (see also [H2]) correct for the QED effects of the decay 
constants, and we only consider the decay constant for e = but m u 7^ rrid in this section. 



In Fig, f 3 , /^(mj obtained by the RBC/UKQCD collaboration [7] is plotted as a function 
of valence light quark mass m x . The sea and valence strange quark masses are fixed. The 
square points are from light sea quark mass m; = 0.01 (~ 40MeV) while the circle data are for 
mi = 0.005(~ 22 MeV). The curves are from the partially quenched SU(2) ChPT fits. The 
upper two curves denote fxi^x) at fixed degenerate sea quark masses mi = 0.01 (upper) 
and 0.005 (lower), while the doted curve is evaluated for unitary quark mass, m x = mi. The 
lower three, almost degenerate, curves are for m L = 0.7 m ud , m ud , and 1.3 m u d- 

The inset magnifies the region close to the physical point. The filled square is fx for equal 
up and down quark masses, mi = m x = m u d- When the valence quark mass is decreased to 
a 30% smaller value, 0.7 x m ud , fx decreases by about 1%, if we simultaneously decrease the 
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light sea quark masses to 0.7 x m ud (empty square). This setting of quark masses (empty 
square) underestimates the value of fx in Nature, since the down sea quark mass is also 
decreased to 0.7 x m ud for the empty square^} 

The more accurate estimation of fx for non-degenerate valence up and down quark masses 
is the empty circle, where the degenerate sea quark mass is fixed to mi = m u d, and only the 
valence quark mass is set to the lighter mass, m x = 0.7m u d. We note that the non-degenerate 
quark mass effect in the sea sector is suppressed by {m u — m d ) 2 , and setting degenerate sea 
quark mass to mi = m ud is a good approximation to estimate the fx shift due to the isospin 
breaking in the up and down quark masses. Because of the (accidental) decrease in the slope 
of fx{m x ) around the physical sea quark mass mi = m u d, the difference between /x(m u d) 
and fx (0.7 x m u d) is only about —0.2%, which is nevertheless sizable compared to the total 
error of 0.4% in the current world average of fx/ A- 

A similar analysis was done in [58], where fx+ was properly estimated at m x = m u . An 
indirect error on fx+ induced from (their) EM uncertainty in m u /md (~ 19%) was estimated 
to be ~ 0.07%. So their shift of fx due to the quark mass difference between m ud and m u 
would be roughly 0.07/0.19 x {m u d/m u — 1) ~ 0.25 % from their value of m u /m u d ~ 0.6. 
This shift is slightly larger than our estimation, 0.2%, in part because m u /m u d in [58] is 
smaller than ours by about 10%. 



4 We thank Chris Sachrajda for pointing this out. 




FIG. 13. Kaon decay constant in pure QCD |7J computed from the same ensembles as used in 
this work. The valence strange quark mass is fixed to 0.03 and the sea strange quark mass is 0.04. 
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IX. NUCLEON MASS SPLITTINGS 

Isospin breaking also occurs in the nucleon system. The proton is slightly lighter than 
the neutron, which makes the proton a stable particle. In conjunction with baryon PQ%PT, 
the lattice simulation helps us understand the relation between the baryon masses and their 
quark content [55] . 

In Nature, m p —m n =— 1.293321(4) MeV as determined by experiments, and it is explained 
by two mechanisms. One is the EM interaction. The proton is a charged particle, but the 
neutron is neutral, so the QED interaction makes the proton heavier. The other is due to 
non-degenerate u, d quark masses. The valence quark content in the proton and neutron is 
uud and udd, respectively. So the proton is lighter than the neutron due to the fact that 
the d quark is heavier than the u. Combining these two effects in our lattice calculations, 
we can compute the p-n mass splitting. 

For non-degenerate quark masses, we study the splitting using the pure QCD configura- 
tions. The nucleon mass in two flavor QCD is given by baryon PQxPT, to NLO [63] . 

m p = M + -(5a + 2f3)m u + -(a + Af3)m d + -a(mj + mi) (60) 
m n = M + -(a + 4/3)m u + ^(5a + 2f3)m d + ^a(mj + m t ) (61) 

where m u , m d are the masses of the valence quarks and rrij, mi are the masses of the sea 
quarks. The mass difference between the proton and the neutron is 

(m p - m n ) {md _ mu) = --(4a - 2(3){m d - m u ). (62) 



We note only the sum of sea quark masses, mj + mi, appears in Eqs. (60) and (61) and the 
difference mj — mi appears first at NNLO in any observable due to the symmetry under 
switching sea up quark to sea down quark. So our degenerate sea up and down quark mass 
is enough to extract the isospin breaking to NLO (We will ignore possible contributions of 
0(e 2 (m u -m d ))). 

Next, we test the EM induced mass splitting on QCD+QED configurations with unitary 
(and therefore degenerate) mass points. The lowest order mass difference is parametrized 
as: 

(m p - m n ) QED = a em (A + A x m ud ) (63) 
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where m u d = (m u + md)/2, and the dependence on a em is made explicit to remind ourselves 
that the splitting vanishes in the absence of QED. 

All of the above LEC's here can be extracted from fits to lattice data. 

We first extract the nucleon masses from the two-point correlation function. The corre- 
lation function measured on the lattice with anti-periodic boundary condition in time has 
the form [64] : 

G{t) = (1 + 7 4)A B+ e- M ^ - (1 - lA )A B+ e- M B+ 

+(1 + l4 )A B -e- M B- W-*) - (1 - l4 )A B -e- M B-\ (64) 

where B + represents the nucleon state which has positive parity and B~ represents the 
excited state of the nucleon which has negative parity. N t is the time-size of the lattice. 
Since the mass of the excited state is much heavier than the ground state, we neglect its 
contribution. The nucleon and anti-nucleon terms left in the correlation function are picked 
up by multiplying G(t) by the projection operator 1 ± 74 and taking the trace. Then we 
average these two terms by taking t — )■ N t — t for the anti-nucleon to improve the statistics of 
our measurements. The ±e trick is also used when QED configurations are included. Finally 
the nucleon masses are extracted from single state fits to point-sink correlation functions as 

G(t) = Ae~ M \ (65) 

where M is the ground state nucleon mass, and A measures the overlap between the nucleon 
state and the nucleon interpolation operator. 

Initially, nucleon correlation functions were computed from the same wall source propa- 
gators used for the meson splitting analysis. However, on the 24 3 ensembles these exhibited 
poor plateaus and had poor signals for the EM neutron-proton mass difference. We then 
switched to box sources (of size 16 3 ), which gave much better plateaus and signals, but only 
on the unitary points because of the additional computational cost. Thus, for the 16 3 and 
24 3 QCD configurations the masses come from wall source correlation functions while for 
the 24 3 QCD+QED configurations, the masses are from box source correlation functions. 



The configuration information of the additional measurements is listed in Tab. |IX| Figure [14 
shows representative plateaus for the sea quark mass 0.005 ensemble. 



The nucleon masses are listed in Tabs. XVII and XVIII They come from a standard x 2 



minimization with correlated fit, and the error on the mass is from the standard jackknife 
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lat ?n soa Tn VSu \ 


Trajectories 


A N meas 


tsrc 


24 3 0.005 0.005 


900-8000 


20 


355 





24 3 0.01 0.01 


1460-8540 


40 


534 


0,16,32 


24 3 0.02 0.02 


1800-3560 


20 


534 0,i 


5,16,24,32,48 


24 3 0.03 0.03 


1260-3020 


20 


534 0,i 


5,16,24,32,48 



TABLE IX. Summary of additional configurations used for the box source nucleon calculation on 
the 24 3 lattices. QCD gauge configurations generated by the RBC and UKQCD collaborations 
[T71l31j. A is the separation between measurements in molecular dynamics time units. The Iwasaki 
gauge coupling is /3 = 2.13. 



method. The results in Tab. XVII for the unitary masses and non-zero a em on the 24 3 
ensembles are consistent with the pure QCD results obtained on the same ensembles reported 
in Ref. [17] , except for the mi = 0.005 case. The proton and neutron masses are about 
three standard deviations smaller than in the pure QCD case, or roughly three percent. It is 
of interest to further investigate how large the EM effect is on the nucleon masses themselves, 
as well as on the mass difference. Of course, in Nature there is no way to measure the nucleon 
mass due to QCD alone. 

Figure [15] shows the mass difference between the proton and neutron due to the QED 
interaction for the unitary points. If there is no QED interaction and m u = ma, then 
m n = m p , which is the result of isospin symmetry. When the QED interaction is included, 
the proton is heavier than the neutron, and the mass difference decreases with quark mass 



as observed in Fig. 15 The 24 3 result is larger than the 16 3 result, once again signaling 
finite volume corrections. This simulation is on the unitary points, but m u ^ in nature. 
When we extrapolate (m p — m n )QED to the physical point, we use the average light quark 
mass m u d, as determined in the previous section. Finally we find that (m p — m n )QED is 



about 0.4 MeV (see Tab. |Xj). From Fig. 15 there is a visible flattening of the splitting at the 
lightest quark mass for the 24 3 lattice size. Using only the lightest two quark masses in the 
extrapolation, we obtain (m p — m n )QED = 0.63(23) MeV. The difference between the two 
results is used to estimate the systematic error in the chiral extrapolation. 

Since the photon is not confined, the EM proton-neutron mass difference could suffer 
from a large finite volume effect. In order to estimate this artifact, we use the Cottingham 
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0.85 




0.85 




FIG. 14. Proton effective masses, 24 3 lattice size, m; = 0.005. The upper panel is for the unitary- 
point and box source. The lower panel is for a non-degenerate case and wall source. 
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formula [9j |65 



5m e i e = 2nam^ ^ 



G E (q) 
L 3 ^ \q\ 



2 



q 2 + 4m 2 2m 2 \y q 2 



Am 2 , 
1 + _ - 1 (66) 



na 1 v-^ 2 / / 4m 2 lit 

^ ma9 --^ T3 ^\ q \G M ( q ) .^i + _-i_- i + ^ /4m2 ), (67) 

where <5m e / e {5m mag ) is the electric (magnetic) contribution to the nucleon mass m. We 
evaluate the above formulae at the physical point, using the dipole form for the nucleon elec- 
tromagnetic form factors, G P E (Q 2 ) = G p M {Q 2 )/ji p = G r ^ 1 {Q 2 ) / \i n = Gd(Q 2 ), where /i p (/i n ), 
are proton (neutron) magnetic moment, and Gn(Q 2 ) = 1 / (1 + Q 2 / A 2 ) 2 with A 2 = 0.71GeV 2 . 
For G E (Q 2 ), we use the Galster parametrization of G n E (Q 2 ) = AQ 2 /(4m 2 + BQ 2 ) ■ G D (Q 2 ) 
with A = 1.70, B = 3.30 [66]. We obtain (m p - m^g^ = 0.04 MeV for 16 3 volume, and 
(m p — m„)o^ ) t '' ) = 0.16 MeV for 24 3 volume. Since the formula yields (m p — m n )Q^ = 0.77 
MeV for the infinite volume limit, the finite volume artifact corresponds to an underestimate 
of 0.73 MeV and 0.61 MeV for 16 3 and 24 3 , respectively. The tendency for the larger volume 
to correspond to larger (m p — m n )qED is qualitatively consistent with the lattice results 
presented here. 

Next we compute the mass splitting due to non-degenerate u and d quark masses, which 
is expected to switch the sign of the mass difference, in accord with Nature. Figure [16] shows 
the fit of the proton and neutron mass difference due to non-degenerate u, d quark masses 
computed on the QCD configurations. The LEC's and values of the splitting at the physical 



point are summarized in Tab. XI Figure 16 confirms that (m p — m n )( md - m \ is proportional 



to m d — m u , which is predicted by baryon PQ^PT (Eq. (62)). The slope is extracted and 
the physical (m p — m n ) ( md -m„) is estimated by setting m d — m u to its physical value, again 
as determined in the previous section. Our result is in good agreement with the one in 
Ref. EES- 



The quark mass dependence of m p —m n is simple in baryon chiral perturbation theory 



to NLO in pure QCD, as seen in Eq. (60). The leading quark mass dependence for the 
EM splitting is unknown, so we assume that it is linear, and at this stage the measured 
values likely can not be used to discern a more complicated form anyway. In contrast, 
chiral perturbation theory for the nucleon mass itself predicts several non-analytic terms at 
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NLO, and the careful extrapolation to the physical point is an important topic of current 
calculations. Because we have few data points, and our quark masses are relatively heavy, 
we do not attempt such an extrapolation here. 

Combining the contributions from the EM interaction and non-degenerate u, d quark 
masses, we give the physical p-n mass splitting. We find m p —m n = —1.93(12) and —2.13(16) 
MeV, for 16 3 and 24 3 lattice sizes, respectively, which is larger than the experimental result 
(—1.293321 (4) MeV), but remarkable given that compared to the mass itself, the splitting 
is a 0.1% effect in Nature. The errors above are statistical only, and their small size is due 
to the facts that the difference is calculated on exactly the same configurations and with the 
±e averaging trick. 

To estimate the systematic error on the EM splitting from the chiral extrapolation we 
take the difference between the extrapolation using all of the data points (on the 24 3 lattice) 
and the lightest two mass points, or roughly 0.3 MeV. The finite volume effects, while quite 
noticeable at the simulated quark masses, are smaller in the quark mass extrapolated result. 
To roughly estimate the finite volume effect, we consider the difference in the 16 3 and 24 3 
results which is about 0.05 MeV when all of the data are used in the fits, and roughly 0.3 
MeV if only the lightest points on the 24 3 lattice are used. In light of the much larger artifact 
predicted by the Cottingham formula, we take the more conservative estimate of 0.3 MeV. 
The finite volume error on the pure QCD splitting appears to be under better control, and 
we simply take the difference of the two as an additional finite volume effect, or ~ 0.25 MeV. 
The QCD splitting depends somewhat strongly on the value of m u —m d , and given the ~ 20% 
uncertainty in this quantity, we estimate the systematic error due to the extrapolation by 
varying m u — rrid over this range. This yields roughly a 0.5 MeV uncertainty. Adding all of 
these errors in quadrature, we find m p — m n = — 2.13(16)(70) MeV. The result and errors are 



summarized in Tab. XII Clearly further calculations are needed, at smaller quark masses 
to improve the extrapolation, with a different lattice spacing to take the continuum limit, 
and on a larger volume to improve the infinite volume extrapolation. 
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lattice size lO 2 ^ M. X 2 /dof (m p - ?ti?i)qed (MeV) 
16 3 2.42(95) 1.26(38) 0.002(96) 0.33(11) 
24 3 2.72(55) 1.80(22) 0.7(1.2) 0.383(68) 

TABLE X. Proton and neutron mass difference due to the QED interaction. The LEC's are 
extracted from the nucleon data on the unitary points. (m p — Tn n ) qed is given at the physical 
quark mass m u d determined in this work. 



lattice size -|(4a - 2/3) x 2 /dof {m p - m n )( md _ mu) (MeV) 



16 3 


-1.452(45) 


1.1(1.2) 


-2.265(70) 


24 3 


-1.612(92) 


0.06(24) 


-2.51(14) 



TABLE XL Proton-neutron mass difference due to non-degenerate u, d quark masses, computed 
on QCD configurations only. (m p — m n )( md _ mu ) is calculated at the physical value of (md — m u ) 
determined in this work. 




FIG. 15. The proton- neutron mass difference due to the QED interaction computed for unitary 
points. 24 3 (circle) and 16 3 (square) lattice sizes. The solid (dashed) line corresponds to a linear 
fit to the 24 3 (16 3 ) data points. 
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10 20 30 40 50 60 
(m d -m u )/Zm (MeV) 



FIG. 16. The proton-neutron mass difference for m u — / and e = 0. The solid (dashed) line 
corresponds to a linear fit to the 24 3 (16 3 ) data points, shown by circles (squares). 



lattice size m p — m n (MeV) fit error (MeV) finite vol. error (MeV) 
16 3 -1.93(12) 

24 3 -2.13(16) 0.58 0.39 

TABLE XII. Estimated result of the proton- neutron mass difference in Nature (systematic errors 
as described in the text). 
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X. CONCLUSION 

In this paper we have investigated the EM mass splittings of the low-lying hadrons from 
first principles in the framework of lattice QCD+QED. Our simulations were based on the 
2+1 flavor DWF QCD configurations generated by the RBC and UKQCD collaborations 
and quenched, non-compact, QED configurations generated by us. The mass splittings could 
be determined with very high statistical accuracy since the QCD part of the fluctuations 
in the hadron masses largely cancels in the splittings. The precision is further enhanced by 
applying our ±e trick [2J [35] to cancel 0(e)-noise on each configuration, before averaging 
over the QCD ensemble. The statistical errors on the pseudo-scalar splittings are at an 
impressive sub-one-percent level, as are the errors on the masses themselves. 

The explicit chiral symmetry breaking induced by the finite extra 5th-dimension of DWF 
was studied in detail and shown to be under good control. This is important because the 
leading 0(a em m ms ) effect is comparable in size to the physical effects under investigation. 

We fit the pseudo-scalar meson mass-squared splittings to the theoretical predictions of 
partially-quenched chiral perturbation theory, including photons, to extract the EM low 
energy constants of the effective theory, up to NLO. We presented new analytic results 
for the kaon mass-squared in Sec. [TT] and in the Appendix |Bj The fits were done to both 
SU(3)l x SU(3)r and SU(2)l x SU(2)R-plus-kaon theories, the latter being necessary to 
determine the strange quark mass [7t lTH4T9] . When using the finite volume PQxPT formulas, 
we found that the NLO corrections relative to LO are about 25% for the physical pion masses, 
neglecting 0(ag m )-terms in the ir° mass that come from the axial anomaly (disconnected 
graphs) and are expected to be small [TT]. Simple linear fits also work as well as the 
complicated NLO chiral perturbation theory ones, as has been seen in the case of pure 
QCD [3 [ToTTTTj . Indeed, our data do not show significant curvature, so while they do not 
seem to require the presence of chiral logs from a theoretical point of view, they are consistent 
with them. The EM splittings and LEC's are significantly affected by the finite volume of 
the lattice, as expected since the long range interactions of the photons are not confined. For 
our final values, we used the finite volume formulas for the chiral logs computed in Ref . [10] . 
The lattice-extracted, SU(3)l x SU(3)r LEC's were found to be somewhat inconsistent with 
the result of the phenomenological analysis in Ref. [TT] , although the latter were fit using an 
ad hoc set of choices for the LEC's. This may also be due to a lack of convergence of SU(3) 
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chiral perturbation theory in the range of quark masses used here, or finite volume effects, 
or both. 

The masses of the light quarks were also determined from our calculation. This is the 
first time EM interactions have been included directly in the quark masses determined from 
2+1 flavor calculations. We employed the physical masses of the n ± , K° and K mesons 
as input to fix the quark masses in PQxPT. The SU(2) L x SU(2) R -plus-kaon theory was 
used to quote our final values since the physical strange quark mass is outside the range of 



convergence of SU(3)l X SU(3)r chiral perturbation theory. They are given in Eqs. (32)- 



(36), along with statistical and systematic errors. The down- up mass difference and quark 



mass ratios are given in Eqs. (35)-(38). These quark masses, up to EM effects, are consistent 
with the pure QCD values given in Ref. [7J, which is not surprising since the pure QCD LEC's 
were taken from an identical analysis of extended ensembles of configurations [17] used there. 
Concerning the solution of the strong CP problem, it is of interest that our value for the up 
quark mass is different from zero by many (~ 6 — 7) standard deviations. 

The Dashen term, or LO EM contribution to the pion mass difference is (m„± — m w o )qed = 
3.38(23) MeV in our calculation, coming from the SU(2) chiral perturbation theory, finite- 
volume-corrected fit, which is our most reliable one. The error is statistical only. However, 
the value from the linear chiral fit agrees within errors. It is also consistent with the values 
of m^± in the chiral limit recently reported in [54"t [55] , but somewhat smaller than the value 
from phenomenology and SU(3) chiral perturbation theory [11] and the value we reported 
for two flavor QCD in Ref. [2] . This suggests that NLO contributions at the physical quark 
masses may be as large as 25% of the total pion mass difference, and approximating the 7r° 
mass from the LEC's computed here, we find the LO+NLO EM contribution at the physical 
point is m n + — m n o = 4.50(23) MeV. Phenomenology predicts that a small part of the NLO 
correction is due to m u — rrid ^ 0, 0.17(3) MeV [56] and 0.32(20) MeV [57]. Similarly, we find 
for the kaons that the pure EM mass difference is (m#-± — ttir-c^qed — 1-87(10) MeV, while 
the contribution from m u — rrid ^ is —5.840(96) MeV. While these values are interesting, 
there is still systematic uncertainty in them which can only be removed by calculations with 
lighter quark masses and larger volumes. 

Finally, we also computed the proton-neutron mass difference, again for the first time 
in 2+1 flavor QCD+QED. Our result is somewhat bigger than the experimental one, but 
encouraging. We found m p — m n = 0.383(68) MeV for the EM mass splitting, and —2.51(14) 
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MeV from m u 7^ m^, both on the larger lattice (errors are statistical). Part of the systematic 
error, stemming mainly from finite volume and chiral extrapolations of the splittings, was 
estimated. The total splitting was found to be m p — m n = — 2.13(16)(70) MeV, where 
the first error is statistical and the second, part of the systematic error. The central value 
is from the 24 3 lattice; we have not attempted either continuum limit or infinite volume 
extrapolations. The sign and relative size of the EM effect compared to the — m u mass 
difference effect is as expected. 

In this work, quenched QED configurations were used to account for the EM interactions 
of the valence quarks, i.e., the sea quarks were neutral in our calculation. The systematic 
error due to this approximation can be removed by the re- weighing method j2H |26]. We 
are now undertaking such a study. In similar spirit to the most recent RBC/UKQCD pure 
QCD calculation [TBTfTT] on a finer lattice ensemble, a w 0.086 fm, the analysis presented 
here is being replicated on those ensembles in order to take the continuum limit. Similarly, 
calculations on a third set of ensembles being generated by the RBC and UKQCD collabo- 
rations are on-going, with a new modified Iwasaki gauge action [67], to better explore the 
chiral regime. 
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Appendix A: Tables of hadron masses and splittings 

TABLE XIII: Summary of pseudo-scalar meson masses ob- 
tained from fits to the pseudo-scalar two-point correlation 
functions on the 16 3 QCD+QED lattice configurations. Fit 

range is 9 < t < Nt/2 in each case. m sea is the mass of the 
light quark in the sea sector, mi and 777,3 are the masses of the 
valence quarks. The mass of the strange sea quark is fixed at 
0.04. "cov" and "uncov" refer to covariant and uncovariant 
fits, respectively. 

777g ea mi m 3 gi(^) q 3 (±) M(cov) x 2 /dof(cov) M(uncov) \ 2 7 'do j '(uncov) 

() 2);!()! , 1.21(89) 0.2420(16) 0.02(4) 

0.2837(13) 1.39(96 0.2828 15) 0.03(5 

0.3196(12) 1.50(99 0.3188(14) 0.04(6) 

0.2427 14) 1.21(89) 0.2417 16) 0.02(4 

0.3193(12) 1.55(1.01) 0.3185(13) 0.04(6) 

0.3518 11) 1.62(1.03 0.3510 12) 0.04(6 

0.3817(11) 1.65(1.05) 0.3809(11) 0.04(6) 

0.2834 13) 1.39(96) 0.2825 15) 0.03(5 

0.2836 13) 1.39(96 0.2827 15) 0.03(5 

0.3196(12) 1.50(99) 0.3187(14) 0.04(6) 

0.3518 11) 1.62(1.03) 0.3509 12) 0.04(6 

0.2430(14) 1.22(89) 0.2420(16) 0.02(4) 

0.2838 13) 1.39(96 0.2828 15) 0.03(5 

0.3197(12) 1.50(99) 0.3188(14) 0.03(6) 

0.3194 12) 1.55(1.01) 0.3186 13) 0.04(6 

0.3519(11) 1.62(1.04) 0.3510(12) 0.04(6) 

0.3194 12) 1.50(99) 0.3185 14) 0.04(6 

0.3818(10) 1.66(1.05) 0.3810(11) 0.04(6) 

0.2434 14) 1.21(89) 0.2424 16) 0.02(4 

0.2841 13) 1.39(96 0.2832 15) 0.03(5 

0.3200(12) 1.50(99) 0.3191(14) 0.04(6) 

0.3198 12) 1.55(1.01) 0.3189 13) 0.04(6 

0.3522(11) 1.62(1.04) 0.3514(12) 0.04(6) 

0.3821 11) 1.66(1.05 0.3813 11) 0.04(6 

0.3191(12) 1.55(1.01) 0.3182(13) 0.04(6) 

0.3516 11) 1.62 1.03 0.3507 12) 0.04(6 

0.3814(11) 1.65(1.04) 0.3806(11) 0.04(6) 

0.2842 13) 1.40(96) 0.2833 15) 0.03(5 

0.2844(13) 1.40(96) 0.2835(15) 0.03(5) 

0.3201 12) 1.51(1.00) 0.3192 14) 0.03(6 

0.3204 12) 1.51 1.00 0.3195 14) 0.03(6 

0.3200(12) 1.56(1.02) 0.3191(13) 0.04(6) 

0.3524 11) 1.63(1.04 0.3515 12) 0.04(6 

0.3525(11) 1.63(1.04) 0.3517(12) 0.04(6) 

0.3823 10) 1.66(1.05) 0.3815 11) 0.04(6 

0.2437(14) 1.22(89) 0.2427(16) 0.02(4) 

0.2843 13) 1.40(96) 0.2834 15) 0.03(5 

0.3201(12) 1.50(1.00) 0.3193(14) 0.04(6) 

0.3201 12) 1.56(1.01 0.3192 13) 0.04(6 

0.3524 11) 1.63 1.04 0.3516 12) 0.04(6 
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TABLE Xlll - continued from previous page 



X* / do f (coy 

1.66(1.05) 
1.22(90) 
1.21(89 
1.39(96 
1.39(96) 
1.50(1.00) 
1.50(99) 
1.56(1.01) 
1.63(1.04 
1.63(1.04 
1.67(1.05) 
1.39(96) 
1.50(1.00) 
1.63(1.04) 
1.23(90) 
1.41(96) 
1.52(1.00) 
1.57(1.02 
1.64(1.04 
1.67(1.05) 
1.21(89' 
1.39(96 
1.50(1.00 
1.56(1.02 
1.63(1.04 
1.67(1.05 
1.07(84 
1.02(82 
0.94(79 
1.07(84 
1.01(82 
0.97(80 
0.97(80 
1.02(82 
1.03(82 
0.94(79 
0.97(80 
1.07(84 
1.03(83 
0.95(79 
1.02(82 
0.97(80 
0.94(79 
0.97(80 
1.07(84 
1.03(82 
0.95(79 
1.02(82 
0.97(80 
0.97(80 
1.01(82 
0.97(80 
0.96(80 
1.04(83 
1.03(83 
0.95(80 
0.95(79 
1.02(82 



m sea 

~omo 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
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0.020 
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0.020 
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0.020 
0.020 
0.020 
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0.020 
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mi 
O030 
0.010 
0.010 
0.010 
0.020 
0.010 
0.030 
0.020 
0.020 
0.030 
0.030 
0.020 
0.030 
0.030 
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0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.020 
0.010 
0.010 
0.010 
0.020 
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0.010 
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0.010 
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-1 
-1 
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-1 
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-1 
1 

-1 



M(cov 
0.3824(10 
0.2436(14 
0.2443 14 
0.2849(13 
0.2851 13 
0.3208(12 
0.3210 12 
0.3206(12 
0.3530 11 
0.3532 11 
0.3830(11 
0.2845 13 
0.3205(12 
0.3526 11 
0.2440(14 
0.2847 13 
0.3206(12 
0.3204 12 
0.3528(11 
0.3827 10 
0.2454 14 
0.2860(13 
0.3219 12 
0.3217(12 
0.3542 11 
0.3841(11 
0.2521(14 
0.2918 
0.3271 
0.2519 
0.3267 
0.3587 
0.3881 
0.2916 
0.2918 
0.3270(12 
0.3586(10 
0.2522 14 
0.2919(12 
0.3271 12 
0.3268(11 
0.3587 10 
0.3268 12 
0.3882(10 
0.2526 14 
0.2922(12 
0.3275 12 
0.3271(11 
0.3591 10 
0.3885(10 
0.3265 11 
0.3584(10 
0.3879(10) 
0.2924 12 
0.2925(12) 
0.3275 12) 
0.3278(12) 
0.3273(11) 



12 

12* 

14* 

11 

10) 

10 

12) 

12* 



12) 
11) 
15) 
14) 
12) 
16) 



M(uncov) 
0.3816(11 
0.2426(16 
0.2433 16 
0.2840(15 
0.2841 15 
0.3199(14 
0.3201 14 
0.3198(13 
0.3522(12 
0.3523 
0.3822 
0.2836 
0.3196 
0.3518 
0.2430 
0.2838(14 
0.3197(14 
0.3195 13 
0.3520(12 
0.3819 11 
0.2444 16 
0.2851(15 
0.3211 14 
0.3208(13 
0.3533 12 
0.3833(11 
0.2528 15 
0.2923(13 
0.3274 12 
0.2525(15 
0.3270 11 
0.3588 11 
0.3882(10 
0.2920 13 
0.2922(13 
0.3273 12 
0.3588(11 
0.2528 15 
0.2923(13 
0.3274(12 
0.3270 
0.3589 
0.3271 
0.3883 
0.2532 
0.2927 
0.3277 
0.3274 
0.3592 
0.3887 
0.3267(11 
0.3586(11 
0.3880 10 
0.2928 13 
0.2929(13 
0.3278 12 
0.3281(12 
0.3276(11 



X 



' /dof (uncov) 

0.04(6) 

0.02(4) 

0.02(4 

0.03(5) 

0.03(5 

0.04(6) 

0.04(6 

0.04(6 

0.04(6 

0.04(6 

0.04(6 

0.03(5 

0.03(6) 

0.04(6 

0.02(4 

0.03(5 

0.03(6) 

0.04(6 

0.04(6 

0.04(6 

0.02(4 

0.03(5) 

0.04(6 

0.04(6) 

0.04(6 

0.04(6) 
0.10(12) 
0.11(13) 
0.10(12) 
0.10(12) 
0.11(12) 
0.09(11) 

0.07(9) 
0.11(13) 
0.11(13) 
0.10(12) 
0.09(11) 
0.10 12) 
0.11(13) 
0.10(12) 
0.11(12) 
0.09(11) 
0.10(12) 

0.07(9) 
0.10(12) 
0.11(13) 
0.10(12) 
0.11(12) 
0.09(11) 

0.07(9) 
0.11(12 
0.09(11 

0.07(9) 
0.12(13 
0.12(13 
0.10(12 
0.10(12 
0.11(12 
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TABLE Xlll - continued from previous page 



X /do f (coy 

0.97(80 
0.97(80 
0.97(80 
1.08(84 
1.03(83 
0.95(79 
1.02(82 
0.97(80 
0.97(80 
1.08(85 
1.08(84 
1.03(83 
1.03(83 
0.95(79 
0.95(79 
1.02(82 
0.97(80 
0.97(80 
0.97(80 
1.03(83 
0.95(79 
0.97(80 
1.09(85 
1.04(83 
0.96(80 
1.03(83 
0.98(81 
0.97(80 
1.08(85 
1.04(83 
0.96(80) 
1.03(82 
0.98(81) 
0.98(81) 
2.00(1.16) 
2.00(1.16 
2.02(1.16 
2.00(1.16 
2.05(1.17 
2.12(1.19 
2.20(1.21 
2.00(1.16 
2.00 1.16 



m sea 

^020 
0.020 
0.020 
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0.020 
0.020 
0.020 
0.020 
0.020 
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0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 



mi 

omo 

0.030 
0.030 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.020 
0.010 
0.030 
0.020 
0.020 
0.030 
0.030 
0.020 
0.030 
0.030 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.020 
0.010 
0.010 
0.010 
0.020 
0.020 
0.010 
0.030 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 



0.030 
0.020 
0.030 
0.010 
0.020 
0.030 
0.020 
0.030 
0.030 
0.010 
0.010 
0.020 
0.010 
0.030 
0.010 
0.020 
0.030 
0.020 
0.030 
0.010 
0.010 
0.020 
0.010 
0.020 
0.030 
0.020 
0.030 
0.030 
0.010 
0.020 
0.030 
0.020 
0.030 
0.030 
0.010 
0.020 
0.030 
0.010 
0.020 
0.030 
0.030 
0.020 
0.020 
0.030 
0.030 
0.010 
0.020 
0.030 
0.020 
0.030 
0.030 
0.030 
0.010 
0.020 
0.030 
0.020 
0.030 
0.030 



-2 
2 
-2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 








1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 



-1 
1 

-1 








1 

-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 




-2 
-2 
-2 
-2 
-2 
-2 
2 
2 
2 
2 
2 
2 
1 
1 
1 

1 
1 
1 





1 
1 
1 
1 
1 



1 

-1 
-1 
-1 
-1 
-1 
-1 



M(cov 
0.3592(10 
0.3593(10 
0.3887 10 
0.2529(14 
0.2925 12 
0.3276(12 
0.3275 11 
0.3593(10 
0.3889 10 
0.2528 14 
0.2535(14 
0.2931 
0.2933 
0.3283 
0.3285 
0.3281(11 
0.3599(10 
0.3600 10 
0.3895(10 
0.2927 12 
0.3279 12 
0.3595(10 
0.2531(14 
0.2928(12 
0.3280 11 
0.3277(11 
0.3596(10 
0.3891 
0.2547 
0.2943 
0.3295 
0.3292 
0.3611 
0.3906 
0.2505 
0.2901(13 
0.3256(12 
0.2503 15 
0.3249(12 
0.3570 11 
0.3865(10 
0.2899 13 
0.2901 13 
0.3255(12 
0.3569 
0.2506 
0.2902 
0.3256 
0.3249(12 
0.3570(11 
0.3253 12 
0.3866(10 
0.2510 15) 
0.2905 13) 
0.3259(12) 
0.3253 12) 
0.3574(11) 
0.3870(10) 



10 

14) 

12 ! 

12 ; 
11 ! 
10) 
10) 
15 



2.02( 
2.12( 
2.00( 
2.00( 
2.02( 



1.16) 

1.19 

1.16) 

1.16 

1.16 



2.06(1.17 
2.12(1.19 
2.02(1.16 



2.20 
2.00 
2.00 
2.02 



1.21 
1.16 
1.16 
1.16 



2.06(1.17 
2.12(1.19 
2.20(1.21) 



M(uncov) 
0.3594(11 
0.3595(11 
0.3889 10 
0.2536(15 
0.2929 13 
0.3279(12 
0.3277(11 
0.3595 
0.3890 
0.2534 
0.2542 
0.2935 
0.2937 
0.3285 
0.3288 
0.3283(11 
0.3601(11 
0.3602 11 
0.3896(10 
0.2931 13 
0.3282 12 
0.3596(11 
0.2538 15 
0.2932(13 
0.3283(12 
0.3280 
0.3598 
0.3892 
0.2554 
0.2947 
0.3297 
0.3294 
0.3613 
0.3907 
0.2507 
0.2901 
0.3254 
0.2504(16 
0.3247(12 
0.3567(12 
0.3862 
0.2898 
0.2901 
0.3253 
0.3567(12 
0.2507(16 
0.2901 14 
0.3254(13 
0.3248 12 
0.3568(12 
0.3251 13 
0.3863(11 
0.2511 16 
0.2905 14) 
0.3258(13) 
0.3251 12) 
0.3571(12) 
0.3866(11) 



X 



'/dof (uncov) 

0.09(11) 
0.09(11) 
0.07(9) 
0.10(12) 
0.12(13) 
0.10(12) 
0.11(12) 
0.09(11) 
0.07(9) 
0.10(12) 
0.10(12) 
0.12(13) 
0.12(13) 
0.10(12) 
0.10(12) 
0.11(12) 
0.09(11) 
0.09(11) 
0.07(9) 
0.12(13) 
0.10(12) 
0.09(11) 
0.10(12) 
0.12(13) 
0.10(12) 
0.11(12) 
0.09(11) 
0.07(9) 
0.11(13) 
0.12(13) 
0.10(12) 
0.11(12) 
0.09(11) 
0.07(9) 
0.16(12) 
0.12(10) 
0.10(9) 
0.16(12) 
0.10(9) 
0.09(8 
0.09(7) 
0.12(10) 
0.12(10) 
0.10(9) 
0.09(8) 
0.16(12) 
0.12(10) 
0.10(9) 
0.10(9 
0.09(8) 
0.10(9 
0.09(7) 
0.16(12) 
0.12(10) 
0.10(9) 
0.10(9 
0.09(8) 
0.09(7) 
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TABLE Xlll - continued from previous page 
m sea mi m 3 gi(^) ggTjO M(cov) x 2 /do/(cov) M(uncov) x 2 /dof(uncov) 



0.030 0.020 


0.020 








0.32461 


12 


1 2.051 


'1. 


.17 


1 0.32451 


12 


) 0.10 




0.030 0.020 


0.030 








0.3567( 


11 


1 2.12( 


>1. 


.19 


1 0.3564( 


12 


) 0.09 




0.030 0.030 


0.030 








0.3863( 


10 


1 2.201 


1. 


.21 


1 0.38601 


11 


1 0.09(7) 


0.030 0.010 


0.020 


-2 


-1 


0.29071 


13 


1 2.0K 


>1. 


.16 


1 0.29061 


14 


) 0.12(10) 


0.030 0.020 


0.010 


2 


1 


0.2908( 


13 


1 2.0K 


1. 


.16 


) 0.29081 


; i4 ; 


) 0.12(10) 


0.030 0.010 


0.030 


-2 


-1 


0.3260( 


12 


1 2.041 


*1. 


,17! 


1 0.32591 


13 


1 0.10(9) 


0.030 0.030 


0.010 


2 


1 


0.3263( 


12 


1 2.031 


1. 


.17 


1 0.3261( 


13 


) 0.10 




0.030 0.020 


0.020 


-2 


-1 


0.3255( 


12 


1 2.06( 


'l. 


■17! 


1 0.32531 


12 


) 0.10 


f\ 


0.030 0.020 


0.030 


-2 


-1 


0.35751 


11 


\ 2.13( 


> 1. 


.19 


1 0.35731 


ir 


) 0.09 


8 


0.030 0.030 


0.020 


2 


1 


0.35761 


> 11 < 


\ 2.13( 


>1. 


.19 


1 0.35741 


11' 


) 0.09 


8 


0.030 0.030 


0.030 


-2 


-1 


0.38721 


"10, 


1 2.2K 


>1. 


.21 


) 0.38681 


hi' 


) 0.09 




0.030 0.010 


0.010 


2 





0.25131 


15 


\ 2.001 


>1. 


.16 


1 0.25141 


16* 


) 0.16( 




0.030 0.010 


0.020 


2 





0.2908( 


13 


1 2.0K 


;i. 


.16 


1 0.2907( 


14' 


) 0.12( 


11] 


0.030 0.010 


0.030 


2 





0.3261( 


'12; 


1 2.031 




.17; 


1 0.32591 


13* 


1 0.10 




0.030 0.020 


0.020 


2 





0.32561 


> 12 < 


1 2.06( 


1. 


.17 


1 0.32551 


12' 


) 0.10 




0.030 0.020 


0.030 


2 





0.35761 


11 


1 2.131 


> i. 


.19 


) 0.35741 


> 12 ; 


) 0.09 


\i 


0.030 0.030 


0.030 


2 





0.3873( 


10 


1 2.211 


'i. 


•2l! 


) 0.38701 


11 


1 0.09 




0.030 0.010 


0.010 


2 


1 


0.25121 




1 2.0K 


;i- 


,ie; 


1 0.2513( 


; ie ; 


) 0.16( 


12) 


0.030 0.010 


0.010 


2 


-1 


0.25191 


15 


1 2.001 


i. 


.16 


) 0.25211 


16 


) 0.16( 


12 


0.030 0.010 


0.020 


2 


-1 


0.2914( 


13; 


1 2.00( 


1. 


.16 


1 0.29141 


14 


) 0.12 


10) 


0.030 0.020 


0.010 


2 


-1 


0.2916( 


13 


1 2.001 


'i. 


.16 


1 0.2915( 


14' 


) 0.12 


10 


0.030 0.010 


0.030 


2 


-1 


0.3267( 


IT 


1 2.031 


'i. 


.16 


) 0.32661 


^3' 


) 0.10(9) 


0.030 0.030 


0.010 


2 


-1 


0.32701 


12 


1 2.02( 


: i. 


,ie; 


1 0.32681 


> 13 ! 


) 0.10 




0.030 0.020 


0.020 


2 


-1 


0.32621 


12 


1 2.061 


'i. 


.17: 


1 0.3261( 


12 


) 0.10 




0.030 0.020 


0.030 


2 


-1 


0.3582( 


11 


. 2.131 


*i. 


.19 


1 0.35801 


12 


1 0.09 


a 


0.030 0.030 


0.020 


2 


-1 


0.35841 


11 


1 2.131 


'i. 


.19 


1 0.35811 


12 


) 0.09 




0.030 0.030 


0.030 


2 


-1 


0.3879( 


10; 


1 2.2K 


> i. 


.21 


) 0.38761 


ir 


) 0.09(7) 


0.030 0.020 


0.010 


-2 





0.2910( 


13, 


1 2.001 


'i. 


.16 


1 0.29091 


14' 


) 0.12(10) 


0.030 0.030 


0.010 


-2 





0.32641 


'12' 


1 2.02( 


1. 


,ie; 


1 0.3263( 


13* 


) 0.10 




0.030 0.030 


0.020 


-2 





0.35781 


> 11 < 


1 2.13( 




.19 


1 0.35751 


12' 


) 0.09 


(3 


0.030 0.010 


0.010 


-2 


-2 


0.25151 


15 


1 2.02( 


> i. 


.16 


1 0.25171 


16* 


) 0.16(12) 


0.030 0.010 


0.020 


-2 


-2 


0.29111 


13 


1 2.02( 


'i. 


.16 


1 0.2911( 


14' 


) 0.12( 


10) 


0.030 0.010 


0.030 


-2 


-2 


0.3265( 


12 


1 2.04( 


'i. 


.17 


1 0.32641 


13' 


) 0.10 


(9) 


0.030 0.020 


0.020 


-2 


-2 


0.32591 


11 


1 2.07( 


1. 


.18 


1 0.32571 


'12' 


) 0.10 


9 


0.030 0.020 


0.030 


-2 


-2 


0.35801 


11 


1 2.14( 




.19 


) 0.3577( 




) 0.10 


8 


0.030 0.030 


0.030 


-2 


-2 


0.3876( 




1 2.211 




,21 


i 0.38721 




1 0.09(7) 


0.030 0.010 


0.010 


-2 


2 


0.25311 


$ 


1 2.00( 




.16 


1 0.25331 




) 0.16(12) 


0.030 0.010 


0.020 


-2 


2 


0.29261 


3 


1 2.001 




.16 


1 0.29261 




) 0.12(10) 


0.030 0.010 


0.030 


-2 


2 


0.3280( 




1 2.031 




,16 


1 0.32781 


} M 


1 0.10(9) 


0.030 0.020 


0.020 


-2 


2 


0.3273( 




1 2.061 


} i. 


.17 


1 0.32721 


'\T 


) 0.10 


(9) 


0.030 0.020 


0.030 


-2 


2 


0.35941 


\r 


1 2.14( 


} i. 


.19 


) 0.35911 




) 0.09 


8) 


0.030 0.030 


0.030 


-2 


2 


0.38901 


10 


1 2.221 


1. 


,22 N 


1 0.38871 


If 


) 0.09 


7 
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TABLE XIV: Same as for Tab. |XIII[ except for lattice size 
24 3 . 



11 1 sea " t l 


'"3 




) <m 3 j 


M(cov) 


X 2 / do f {cow) M(uncov 


) X 2 /^o/(uncov) 


n nn^ n nn^ 


n non 
u.uzu 


-i 


n 
u 


0.2605 


>i< 


0.92 




0.2605 




0.07 




n nn^ n nn^ 
u.uuo u.uuo 


n n^n 

U.UoU 


- 1 


A 

u 


0.2984 


A 


0.99 


42 


0.2985 


,5 


0.05 




n nn^ n nm 
u.uuo U.U1U 


U.U1U 


- 1 


A 

u 


0.2393 


>i< 


0.97 


>£ 


0.2393 




0.08 


i 


n nn^ n ni n 
u.uuo U.U1U 


n non 
u.uzu 


- 1 


A 

u 


0.2798 


A 


0.94 


41 


0.2799 




0.06 




n nn^ n ni n 
u.uuo U.U1U 


n mn 

U.UOU 


- 1 


A 
U 


0.3154 


>i< 


0.98 


>$ 


0.3156 




0.04 




n nn^ n non 
u.uuo u.uzu 


n mn 
u.uzu 


- 1 


A 

u 


0.3153 


3, 


0.89 


40 


0.3154 


,5 


0.04 




n nn^ n non 
u.uuo u.uzu 


n n"?n 

U.UoU 


- 1 


A 

u 


0.3475 


3 


0.91 


40 


0.3477 


4 


0.04 


3) 


n nn^ n m n 
u.uuo U.U1U 


n nfn 

U.UOU 


n 
u 


n 
u 


0.3152 


'4 


0.98 


42 


0.3154 


5 


0.04 


4) 


n nn^ n n^n 

U.UUO U.UOU 


n n"?n 

U.UoU 


- 1 


A 

u 


0.3773 


*3 


0.91 


*40 


0.3776 


■i) 


0.04 


3) 


n nn^ n nm 
u.uuo U.UUl 


n nm 

U.UUl 


- 1 


A 

u 


0.1399 


& 


1.05 


>J?< 


0.1396 


>l< 


0.14(10) 


n nn^ n nn^ 
u.uuo u.uuo 


n nn^ 
u.uuo 


- 1 


1 


0.1912 


A 


0.96 


41, 


0.1912 


,6, 


0.13(10 


n nn^ n nn^ 
u.uuo u.uuo 


n ni n 

U.U1U 


i 

- 1 


1 

1 


0.2169 


>i< 


0.95 


>fo< 


0.2169 




0.10 




n nn^ n nn^ 
u.uuo u.uuo 


n non 
u.uzu 


- 1 


1 

1 


0.2610 


A 


0.91 


40, 


0.2611 


5 


0.07 




n nn^ n nn^ 
u.uuo u.uuo 


n n^n 

U.UOU 


- 1 


1 

1 


0.2989 


4< 


0.99 




0.2990 


>^ 


0.05 


i 


n nn^ n nm 
u.uuo U.U1U 


n m n 

U.U1U 


- 1 


1 

1 


0.2398 


A 


0.97 


41 


0.2398 


,5 


0.08 




n nn^ n ni n 
u.uuo U.U1U 


n non 
u.uzu 


- 1 


1 

1 


0.2803 


>1< 


0.93 




0.2804 




0.06 


i 


n nn^ n ni n 
u.uuo U.U1U 


n n"?n 

U.UoU 


- 1 


1 

1 


0.3160 




0.97 


42 


0.3161 


,5, 


0.04 




n nn^ n non 
u.uuo u.uzu 


n non 
u.uzu 


- 1 


1 

1 


0.3157 


3 


0.88 


40 


0.3159 


5 


0.04 


3 


n nn^ n non 
u.uuo u.uzu 


n mn 

U.UoU 


- 1 


1 

1 


0.3480 


'3 


0.90 


40 


0.3482 


4 


0.04 


3 


n nn^ n non 
u.uuo u.uzu 


n non 
u.uzu 


n 

U 


n 
u 


0.3150 


3 


0.89 


40 


0.3151 


5 


0.04 


3 


n nn^ n n^n 

U.UUO U.UOU 


n n"?n 

U.UoU 


- 1 


i 


0.3778 


; ! 


0.91 


>S< 


0.3780 


IS, 


0.04(3) 


n nn^ n nm 

U.UUO U.UUl 


n nm 

U.UUl 


- 1 


i 


0.1405 




1.05 


,43 


0.1402 


X 


0.15(10) 


n nn^ n nn^ 
u.uuo u.uuo 


n nn^ 
u.uuo 


- 1 


-i 


0.1907 


>^< 


0.97 




0.1907 




0.13(10) 


n nn^ n nn^ 
u.uuo u.uuo 


n mn 

U.U1U 


- 1 


-i 


0.2165 


A 


0.96 


41, 


0.2165 


5 


0.10(9) 


n nn^ n nn^ 
u.uuo u.uuo 


n non 
u.uzu 


- 1 


-i 


0.2606 


>i< 


0.92 




0.2606 




0.07 




n nn^ n nn^ 
u.uuo u.uuo 


n nfn 

U.UOU 


i 

- 1 


-i 


0.2985 


A 


0.99 


42 


0.2986 


,5, 


0.05 




n nn^ n m n 
u.uuo U.U1U 


U.U1U 


- 1 


-i 


0.2393 


>i< 


0.97 


>S< 


0.2394 




0.08 


i 


n nn^ n m n 
u.uuo U.U1U 


n non 
u.uzu 


- 1 


-i 


0.2798 




0.94 


41, 


0.2799 


,5, 


0.06 




n nn^ n ni n 
u.uuo U.U1U 


n n^n 

U.UOU 


- 1 


-i 


0.3155 


4 


0.98 


42' 


0.3157 


5 


0.04 


4 


n nn^ n non 
u.uuo u.uzu 


n non 
u.uzu 


- 1 


-i 


0.3153 


'3' 


0.89 


i 40 


0.3154 


5 


0.04 


3 


n nn^ n non 
u.uuo u.uzu 


n n"?n 

U.UoU 


- 1 


-i 


0.3476 


3 


0.91 


40 


0.3478 


4 


0.04 


3) 


n nn^ n non 
u.uuo u.uzu 


n n"?n 

U.UoU 


U 


A 

u 


0.3473 


>l 


0.91 


>X 


0.3475 


>! 


0.04 


i 


n nn^ n n^n 

U.UUO U.UoU 


n mn 

U.UOU 


- 1 


1 
-1 


0.3774 


A 


0.92 


,40 


0.3776 


,4; 


0.04 




n nn^ n nm 

U.UUO U.UUl 


n nm 

U.UUl 


-i 


1 
-1 


0.1399 


>% 


1.05 


>$ 


0.1396 


>J< 


0.15( 


10) 


n nn^ n nn^ 
u.uuo u.uuo 


n nn^ 
u.uuo 


U 


A 

u 


0.1904 


A 


0.96 


41 


0.1903 




0.13(10) 


n nn^ n n^n 

U.UUO U.UoU 


n n"?n 

U.UoU 


U 


A 

u 


0.3771 


& 


0.91 


>^< 


0.3773 


'£ 


0.04(3) 


n nn^ n nm 

U.UUO U.UUl 


n nm 

U.UUl 


n 

U 


A 

u 


0.1395 


A 


1.05 


,43 


0.1392 


X 


0.14(10) 


n nn^ n nn^ 
u.uuo u.uuo 


n nn^ 
u.uuo 


n 

U 


1 


0.1907 


>5< 


0.96 




0.1906 


& 


0.13(10) 


n nn^ n nn^ 
u.uuo u.uuo 


n mn 

U.U1U 


U 


1 

1 


0.2164 


A 


0.95 


41 


0.2164 


5 


0.10(9) 


n nn^ n nn^ 
u.uuo u.uuo 


n non 
u.uzu 


U 


1 

1 


0.2606 


4 


0.92 


40 


0.2606 


5 


0.07 


7 


n nn^ n nn^ 
u.uuo u.uuo 


n mn 

U.UoU 


n 

U 


1 

1 


0.2985 




0.99 


42 


0.2986 


5 


0.05 


5 


n nn^ n m n 
u.uuo U.U1U 


n non 
u.uzu 


n 

U 


1 

1 


0.2798 


> 4 < 


0.94 


> 41 < 


0.2799 


5 


0.06 




n nn^ n ni n 
u.uuo U.U1U 


n n"?n 

U.UoU 


U 


1 

1 


0.3155 




0.98 


! 42 < 


0.3157 


>l< 


0.04 


i 


n nn^ n non 
u.uuo u.uzu 


n nfn 

U.UOU 


n 
u 


1 
1 


0.3476 


A 


0.91 


,40, 


0.3477 


,4, 


0.04 




n 005 n 005 


01 


n 

U 


n 


0.2161 




0.95 




0.2161 




0.10 




0.005 0.005 


0.020 








0.2603 




0.92 




0.2603 


i 


0.07 




0.005 0.005 


0.010 


1 





0.2164 




0.95 




0.2164 




0.10 


i 


0.005 0.005 


0.030 








0.2982 


i 


0.99 




0.2983 




0.05 




0.005 0.010 


0.010 








0.2390 




0.97 




0.2391 




0.08 




0.005 0.010 


0.020 








0.2795 


i 


0.94 


ft 


0.2796 




0.06 


IS 


0.005 0.005 


0.010 


-2 





0.2172 


4* 


0.95 


AT 


0.2172 


5 


0.10 


9) 


0.005 0.005 


0.020 


-2 





0.2611 


> 4 < 


0.92 




0.2612 




0.07 


7 


0.005 0.005 


0.030 


-2 





0.2990 


> 4 < 


0.99 


42 < 


0.2991 




0.05 


5 


0.005 0.010 


0.010 


-2 





0.2401 


i 


0.97 




0.2402 




0.08(8) 


0.005 0.001 


0.001 


-2 





0.1410 




1.04 




0.1407 




0.15(11) 
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TABLE XIV - continued from previous page 
7773 Qi(\) Qs{\) M(cov) x 2 Jdofjcov) M(uncov) x' 2 /dof(uncov) 



m. 



7771 



0.06(5 
0.04(4 
0.04(3 
0.04(3 
0.04(3) 
0.13(10) 
0.10(9 
0.10(9 
0.07(7 
0.07(7 
0.05(5 



0.005 0.010 


0.020 


-2 





0.28051 


a: 


1 0.931 


l 4 K 


1 0.28061 


;s 


0.005 0.010 


0.030 


-2 





0.316H 


4 


! 0.971 


M 


1 0.31621 


5 


0.005 0.020 


0.020 


-2 





0.31601 


3 


1 0.89( 


;4o; 


1 0.3162) 


5 


0.005 0.020 


0.030 


-2 





0.34821 


3 


i 0.901 


40 


1 0.3484) 


4 


0.005 0.030 


0.030 


-2 





0.378K 


3 


i 0.9H 


;4o; 


1 0.3784) 


4, 


0.005 0.005 


0.005 


-2 


1 


0.19241 


5 


1 0.961 


41 


1 0.1923) 


6 


0.005 0.005 


0.010 


-2 


1 


0.21801 


4 


1 0.951 


41 


) 0.2180) 


5 


0.005 0.010 


0.005 


-2 


1 


0.218K 


4 


i 0.951 


41 


1 0.2181) 


5 


0.005 0.005 


0.020 


-2 


1 


0.2619( 


4 


1 0.9H 


40 


1 0.2619) 


'5' 


0.005 0.020 


0.005 


-2 


1 


0.262H 


% 


1 0.911 


'40 


1 0.2622) 


5' 


0.005 0.005 


0.030 


-2 


1 


0.29971 


% 


i 0.981 


} 42 


1 0.2998) 


5' 


0.005 0.030 


0.005 


-2 


1 


0.30011 


>4 


1 0.981 


*42 N 


1 0.3002) 


'5' 


0.005 0.010 


0.010 


-2 


1 


0.24081 


% 


1 0.961 




1 0.2409) 


5' 


0.005 0.001 


0.001 


-2 


1 


0.14191 


>6 / 


1 1.041 


'42 N 


) 0.1416) 


7 


0.005 0.010 


0.020 


-2 


1 


0.2812( 


% 


1 0.921 


'40, 


1 0.2813) 


5' 


0.005 0.020 


0.010 


-2 


1 


0.2814( 


l 4 < 


i 0.931 


;4o; 


1 0.2814) 


5; 


0.005 0.010 


0.030 


-2 


1 


0.31681 


4 


1 0.961 


41 


) 0.3169) 


5 


0.005 0.030 


0.010 


-2 


1 


0.31711 


4 


1 0.971 


42 


1 0.3172) 


5 


0.005 0.020 


0.020 


-2 


1 


0.31671 


3 


1 0.881 


[39] 


i 0.3169) 


5 


0.005 0.020 


0.030 


-2 


1 


0.34891 


'3; 


i 0.901 


40 


1 0.3491) 


'A 


0.005 0.030 


0.020 


-2 


1 


0.349H 


3 


1 0.901 


40 


) 0.3492) 




0.005 0.030 


0.030 


-2 


1 


0.37881 


3 


1 0.911 


'40 


1 0.3790) 




0.005 0.005 


0.005 


-2 


-1 


0.1914( 


5 


1 0.971 




1 0.1913) 




0.005 0.005 


0.010 


-2 


-1 


0.21701 


4 


! 0.961 


41 


1 0.2170) 


5 


0.005 0.010 


0.005 


2 


1 


0.217K 


4 


1 0.961 


41 


1 0.2171) 


5 


0.005 0.005 


0.020 


-2 


-1 


0.26101 


4 


1 0.92) 


40 


1 0.2611) 


5 


0.005 0.020 


0.005 


2 


1 


0.26131 


4 


1 0.921 


;4o; 


) 0.2613) 


'5' 


0.005 0.005 


0.030 


-2 


-1 


0.29891 


% 


1 0.991 


; 42 : 


1 0.2990) 


} 5' 


0.005 0.030 


0.005 


2 


1 


0.29921 


>4 / 


1 l.OOl 


; 4 2; 


l 0.2993) 


5' 


0.005 0.010 


0.010 


-2 


-1 


0.23991 


% 


i 0.971 


41 


1 0.2400) 


5' 


0.005 0.001 


0.001 


-2 


-1 


0.1408( 


>6 / 


i 1.04( 


'42 N 


1 0.1404) 


7 


0.005 0.010 


0.020 


-2 


-1 


0.28031 


% 


i 0.941 


> 41 < 


) 0.2804) 


5' 


0.005 0.020 


0.010 


2 


1 


0.28051 


% 


i 0.941 




) 0.2806) 


5' 


0.005 0.010 


0.030 


-2 


-1 


0.3160) 


>4 / 


1 0.98( 


; 4 2; 


1 0.3161) 


5' 


0.005 0.030 


0.010 


2 


1 


0.31621 


% 


1 0.981 


M 


l 0.3164) 


5' 


0.005 0.020 


0.020 


-2 


-1 


0.31591 


> 


1 0.89( 


'40; 


) 0.3160) 


5' 


0.005 0.020 


0.030 


-2 


-1 


0.3481) 


3 


1 0.9H 


40 


1 0.3483) 


4 


0.005 0.030 


0.020 


2 


1 


0.34821 


3 


1 0.911 


;4o; 


) 0.3484) 


4, 


0.005 0.030 


0.030 


-2 


-1 


0.37801 


3 


1 0.921 


40 


1 0.3782) 


'A 


0.005 0.005 


0.005 


2 





0.19161 


5 


1 0.961 


41 


! 0.1915) 


I 6 


0.005 0.010 


0.005 


-2 





0.21731 


4 


i 0.951 


41 


1 0.2173) 


5 


0.005 0.020 


0.005 


-2 





0.26151 


4 


1 0.921 


40 


1 0.2615) 


'5' 


0.005 0.030 


0.005 


-2 





0.29951 


>4 


1 0.991 


42 < 


1 0.2996) 


5' 


0.005 0.020 


0.010 


-2 





0.28071 


% 


1 0.931 


>41 


i 0.2808) 


5 ! 


0.005 0.030 


0.010 


-2 





0.31651 


'4; 


1 0.981 


*42 ( 


1 0.3166) 


5' 


0.005 0.030 


0.020 


-2 





0.34841 


>3 


1 0.901 


'40, 


) 0.3486) 


} 4 


0.005 0.005 


0.005 


2 


-2 


0.19381 


5 


1 0.95( 




1 0.1937) 


6' 


0.005 0.005 


0.010 


2 


-2 


0.21931 


% 


1 0.94( 


41 


1 0.2193) 


5' 


0.005 0.005 


0.020 


2 


-2 


0.26321 




i 0.901 


$ 


1 0.2633) 




0.005 0.005 


0.030 


2 


-2 


0.30111 




1 0.97( 


42, 


1 0.3012) 




0.005 0.010 


0.010 


2 


-2 


0.242H 




1 0.951 




1 0.2421) 




0.005 0.010 


0.020 


2 


-2 


0.28251 


& 


i 0.921 




1 0.2826) 


% 


0.005 0.010 


0.030 


2 


-2 


0.31811 


'4; 


1 0.961 


41 < 


1 0.3183) 


'5' 


0.005 0.020 


0.020 


2 


-2 


0.31791 




1 0.871 


'39 


1 0.3180) 


5' 


0.005 0.020 


0.030 


2 


-2 


0.35021 


3 


1 0.891 


'40, 


1 0.3504) 


4 ! 


0.005 0.030 


0.030 


2 


-2 


0.38001 


'3; 


i 0.901 


40 


1 0.3803) 


''4 


0.005 0.001 


0.001 


2 


-2 


0.14361 


6 


i 1.031 


42 


1 0.1433) 


7 


0.005 0.005 


0.005 


2 


2 


0.19181 


5 


1 0.971 


;4i; 


1 0.1917) 


6' 



0.08(8 
0.15(11) 
0.06(5) 
0.06(5 
0.04(4 
0.04(4 
0.04(3 
0.04(3 
0.04(3) 
0.04(3) 
0.13(10) 
0.10(9) 
0.10(9 
0.07(7 
0.07(7 
0.05(5 
0.05(5 
0.08(8 
0.15(11) 



0.06 
0.06 
0.04 
0.04 
0.04 
0.04 
0.04(3 
0.04(3 
0.13(10) 
0.10(9) 
0.07(7) 
0.05(5) 
0.06(5) 
0.04(4 
0.04(3 
0.13(10) 
0.10(9) 
0.07(7 
0.05(5 
0.08(8 
0.05(5 
0.04(4) 
0.04(3) 
0.04(3) 
0.04(4) 
0.15(11 
0.13(10 
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TABLE XIV - continued from previous page 
7773 Qi(\) Q3{\) M(cov) x 2 Jdofjcov) M(uncov) x 1 / do f {\mcav) 



m. 



7771 



0.005 0.005 


0.010 


2 


2 


0.21751 


;4; 


1 0.961 


l 4 K 


l 0.21751 


fi 


0.005 0.005 


0.020 


2 


2 


0.26151 


4 


! 0.931 


40 


l 0.26161 


5 


0.005 0.005 


0.030 


2 


2 


0.29941 


4 


1 1.001 


42 


1 0.29951 


5 


0.005 0.010 


0.010 


2 


2 


0.24031 


4 


1 0.981 


41 


1 0.24031 


5 


0.005 0.010 


0.020 


2 


2 


0.28081 


4 


i 0.94i 




l 0.28091 


5, 


0.005 0.010 


0.030 


2 


2 


0.31651 


4 


1 0.99( 


M 


1 0.31661 


5 


0.005 0.020 


0.020 


2 


2 


0.31621 


3 


1 0.901 


40 


l 0.31641 


5 


0.005 0.020 


0.030 


2 


2 


0.34851 


3 


1 0.921 


40 


l 0.34871 


4 


0.005 0.030 


0.030 


2 


2 


0.37841 


'3; 


i 0.931 


'4r 


l 0.37861 


'4' 


0.005 0.001 


0.001 


2 


2 


0.14131 


>6 


1 1.031 


> 42 < 


l 0.14091 


7' 


0.010 0.010 


0.010 





1 


0.24091 


3 


1 2.981 


'73' 


l 0.24111 


5' 


0.010 0.001 


0.001 


-1 





0.14141 


5 


1 4.031 


'85) 


1 0.14191 


8' 


0.010 0.010 


0.020 





1 


0.281H 


3 


1 2.771 


> ?1 < 


! 0.28151 


5' 


0.010 0.001 


0.001 


-1 


1 


0.142K 


5 


1 4.051 


'85) 


) 0.14261 


8' 


0.010 0.010 


0.030 





1 


0.31671 


3 


i 2.54( 


68 


1 0.31731 


5' 


0.010 0.001 


0.001 


-1 


-1 


0.14141 


5 


1 4.04i 


;ss; 


1 0.14191 


> 8 ' 


0.010 0.010 


0.010 








0.24061 


3 


1 2.981 




1 0.24081 


5 


0.010 0.020 


0.020 





1 


0.31651 


3 


1 2.411 


66 


l 0.31691 


5 


0.010 0.001 


0.001 








0.14101 


5 


i 4.021 


85 


l 0.14151 


8 


0.010 0.020 


0.030 





1 


0.3488( 


'3; 


i 2.181 


'63 N 


1 0.34931 


'4' 


0.010 0.030 


0.030 





1 


0.37861 


3 


1 1.991 


'60 


l 0.37911 


>4' 


0.010 0.010 


0.020 








0.28081 


3 


1 2.781 


> n < 


) 0.28131 


} 5 : 


0.010 0.010 


0.020 


1 





0.281K 


3 


i 2.77( 




l 0.28151 


5' 


0.010 0.010 


0.030 


1 





0.31661 


3 


1 2.541 


y 


l 0.31721 


5' 


0.010 0.020 


0.030 


1 





0.34871 


3 


1 2.181 


63 


1 0.34921 


4 


0.010 0.010 


0.010 


1 


1 


0.24101 


3 


i 2.97( 


73 


l 0.24111 


5 


0.010 0.010 


0.020 


1 


1 


0.281K 


'3; 


1 2.761 




l 0.28161 


'5' 


0.010 0.010 


0.030 


1 


1 


0.3167) 


3 


1 2.541 


[68] 


1 0.31731 


5' 


0.010 0.020 


0.020 


1 


1 


0.31651 


'3; 


i 2.411 




l 0.31701 


5' 


0.010 0.020 


0.030 


1 


1 


0.34881 


3 


i 2.181 


63 


) 0.34931 


4' 


0.010 0.010 


0.030 








0.31641 


'3; 


1 2.551 


68 


l 0.31701 


5' 


0.010 0.030 


0.030 


1 


1 


0.37861 


3 


1 1.991 


60 


l 0.37911 


} 4 


0.010 0.010 


0.010 


1 


-1 


0.2414( 


3 


1 2.981 


'73' 


1 0.24161 




0.010 0.010 


0.020 


1 


-1 


0.28161 


'3; 


i 2.76( 


71 


1 0.28201 


} 5 


0.010 0.010 


0.030 


1 


-1 


0.31711 


3 


1 2.531 


M 


) 0.31771 


5' 


0.010 0.020 


0.020 


1 


-1 


0.31691 


3 


1 2.401 


;ee; 


1 0.31741 


5' 


0.010 0.020 


0.030 


1 


-1 


0.34921 


3 


i 2.17( 


63 


1 0.34971 


4 


0.010 0.030 


0.030 


1 


-1 


0.379K 


3 


1 1.991 


;eo; 


1 0.37951 


4, 


0.010 0.020 


0.020 








0.31621 


3 


1 2.421 


w 


1 0.31671 


5 


0.010 0.020 


0.030 








0.34851 


3 


1 2.191 




l 0.34901 


4, 


0.010 0.030 


0.030 








0.37831 


3 


i 2.001 


; 6 o: 


l 0.37881 


4 


0.010 0.010 


0.020 


-2 


-1 


0.28171 


'3; 


i 2.741 


71; 


1 0.28211 


'5' 


0.010 0.020 


0.010 


2 


1 


0.28181 


3 


1 2.751 


71 


l 0.28221 


5' 


0.010 0.010 


0.030 


-2 


-1 


0.31721 


3 


1 2.521 


68 


1 0.31771 


5' 


0.010 0.030 


0.010 


2 


1 


0.31741 


'3; 


1 2.521 


;es; 


l 0.31801 


5' 


0.010 0.020 


0.020 


-2 


-1 


0.31711 


3 


1 2.401 




) 0.31761 


5 


0.010 0.020 


0.030 


-2 


-1 


0.34931 


'3; 


1 2.171 


63 


1 0.34981 


4 


0.010 0.030 


0.020 


2 


1 


0.3495( 


3 


i 2.17( 


63 


) 0.34991 


4 


0.010 0.001 


0.001 


-2 





0.14251 


4 


1 4.071 


; 86 ; 


1 0.14301 


2- 


0.010 0.030 


0.030 


-2 


-1 


r\ . i r-7(-\ r\ , 

0.37921 




1 1.98l 


,60, 


1 0.37971 




0.010 0.010 


0.010 


2 





0.24181 




1 2.961 




l 0.24191 




0.010 0.010 


0.020 


2 





0.28181 




1 2.751 




l 0.28221 


J 


0.010 0.010 


0.030 


2 





0.31731 


'3; 


i 2.52( 


'68 N 


l 0.31791 


'5' 


0.010 0.020 


0.020 


2 





0.31731 


3 


1 2.401 


66 


l 0.31771 


5' 


0.010 0.020 


0.030 


2 





0.34951 


3 


1 2.161 


63 


) 0.35001 


4' 


0.010 0.030 


0.030 


2 





0.37941 


'3; 


1 1.981 


60 


l 0.37991 


''4 


0.010 0.010 


0.010 


2 


1 


0.24161 


3 


i 2.951 


73 { 


! 0.24171 


5' 


0.010 0.001 


0.001 


-2 


1 


0.14351 


5 


1 4.091 


'86^ 


1 0.14401 


8' 



0.101 




0.071 


J 


0.051 




0.08( 




0.05( 




0.041 




0.041 




0.041 





15) 
16) 



0.04(3) 
0.15(11 
0.27(18 
0.56(31 
0.23(16 
0.56(32 
0.20(15 
0.56(31 
0.27(18 
0.20(15 
0.55(31 
0.18(15 
0.17 
0.23 
0.23(16 
0.20(15 
0.18(15 
0.27(18 
0.23 
0.20 
0.20 
0.18 
0.20 
0.17 
0.27 
0.23 
0.20 
0.20(15 
0.18(15 
0.17(15 
0.20(15 
0.18(15 
0.17(15 
0.22 
0.22 
0.20 
0.20 
0.20 
0.18 
0.18 
0.57(32 
0.17(15 
0.27(18 
0.22(16 
0.20 
0.20 
0.18 
0.17 
0.27 



16) 
15) 
15) 
15) 
15) 
15) 
18) 
16 
15 



16) 
16) 
15) 
15) 
15) 
15) 
15) 



15) 
15) 
15) 
15) 
18) 



0.58(32) 
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m sea mi 


m 3 




i M(cov) 


X /doj 


f (cov) M(uncov 


) X /dof{ 


uncov) 


0.010 0.010 


0.010 


2 


-1 


0.24251 


■5; 


1 2.961 


73 


1 0.24261 




1 0.27( 


18 


) 


0.010 0.010 


0.020 


2 


-1 


0.28251 


3 


1 2.74( 


71 


1 0.28291 


5 


1 0.22( 


'16 




0.010 0.020 


0.010 


2 


-1 


0.28271 


3 


1 2.74( 




) 0.28311 


5 


1 0.22( 


16 


) 


0.010 0.010 


0.030 


2 


-1 


0.31801 


3 


1 2.50( 


} M 


1 0.31861 


5 


1 0.20( 


'15 




0.010 0.030 


0.010 


2 


-1 


0.31831 


3 


1 2.50i 




1 0.31891 


; s : 


1 0.20( 


15 


) 


0.010 0.020 


0.020 


2 


-1 


0.31801 


3 


1 2.38( 


ee 


) 0.31841 


5 


1 0.20( 


15 




0.010 0.020 


0.030 


2 


-1 


0.35021 


3 


1 2.151 


63 


) 0.35071 


4 


1 0.18( 


15 


) 


0.010 0.030 


0.020 


2 


-1 


0.35031 


3 


1 2.15( 




1 0.35081 


4 


1 0.18( 


15 




0.010 0.030 


0.030 


2 


-1 


0.38011 


'3; 


i 1.97( 


60 


1 0.38061 


'4 N 


1 0.17( 


15 




0.010 0.001 


0.001 


-2 


-1 


0.14221 


5 


1 4.07l 


86 


1 0.14281 


> 8 < 


1 0.56( 


'32 




0.010 0.020 


0.010 


-2 





0.28201 


3 


1 2.75i 


71 


! 0.28241 


5 A 


1 0.22( 


16 




0.010 0.030 


0.010 


-2 





0.31761 


3 


1 2.52( 


W 


1 0.31821 


'5" 


1 0.20( 


15 




0.010 0.030 


0.020 


-2 





0.34961 


3 


1 2.16l 


63 


1 0.35011 


>4 


1 0.18( 


15 




0.010 0.001 


0.001 


2 


-2 


0.14531 


5 


i 4.14( 


'se; 


1 0.14571 


8 


1 0.59( 


'33 




0.010 0.001 


0.001 


2 


2 


0.14271 


5 


1 4.091 


[86] 


l 0.14331 


8 X 


1 0.56( 


'32 




0.010 0.010 


0.010 


-2 


-2 


0.24201 


3 


1 2.92i 


73 


1 0.24211 


5 


1 0.26( 


'18 


) 


0.010 0.010 


0.020 


-2 


-2 


0.282H 


3 


i 2.73( 


70 


1 0.28251 


5 


1 0.22( 


16 




0.010 0.010 


0.030 


-2 


-2 


0.31771 


3 


1 2.511 


68 


1 0.31821 


5 


1 0.20( 


15 


) 


0.010 0.020 


0.020 


-2 


-2 


0.31751 


3 


1 2.391 


66 


i 0.31791 


4 


1 0.20( 


15 




0.010 0.020 


0.030 


-2 


-2 


0.3498( 


'3; 


1 2.161 


'63; 


1 0.35031 


'4 N 


1 0.18( 


15 




0.010 0.030 


0.030 


-2 


-2 


0.37961 


3 


1 1.98( 


60 


) 0.38011 




1 0.17( 


14 




0.010 0.010 


0.010 


-2 


2 


0.24381 


3 


1 2.95( 


>73 


1 0.24391 


> 5 < 


1 0.27( 


'18' 




0.010 0.010 


0.020 


-2 


2 


0.28391 


3 


1 2.71i 


70 


1 0.28431 


5 


1 0.22( 


16 


) 


0.010 0.010 


0.030 


-2 


2 


0.31941 


3 


1 2.46l 


67 


1 0.31991 


5 


1 0.19( 


15 




0.010 0.020 


0.020 


-2 


2 


0.3192) 


3 


1 2.35( 


66 


) 0.31961 


5 


1 0.20( 


15 


) 


0.010 0.020 


0.030 


-2 


2 


0.35151 


3 


i 2.12( 


[62] 


1 0.35201 


4 


1 0.18( 


15 




0.010 0.030 


0.030 


-2 


2 


0.3814( 


'3; 


1 1.951 


59 


1 0.38181 


'4 A 


1 0.17(15 


) 


0.020 0.020 


0.020 


-1 





0.32221 


3 


1 1.21l 


M 


) 0.32201 


> 5 < 


1 0.06 


(6) 




0.020 0.020 


0.020 


-1 


1 


0.32271 


3 


1 1.21( 


46 


1 0.32251 


'5 s 


1 0.06 


6) 




0.020 0.020 


0.020 


-1 


-1 


0.32221 


3 


1 1.20i 


} M 


1 0.32201 


5 A 


1 0.06 


(6) 




0.020 0.020 


0.020 








0.32191 


'3; 


i 1.21( 


'46 N 


1 0.32171 


} 5" 


1 0.06 


(6) 




0.020 0.020 


0.020 


-2 





0.32301 


3 


1 1.20l 


'46 


) 0.32281 


5 


1 0.06 


6) 




0.020 0.020 


0.020 


-2 


1 


0.32371 


3 


1 1.20i 


} 46 


1 0.32351 


5 A 


1 0.06 


(6) 




0.020 0.020 


0.020 


-2 


-1 


0.32281 


'3; 


1 1.19( 


;4e; 


1 0.32261 


} 5" 


1 0.06 


6) 




0.020 0.020 


0.020 


2 


-2 


0.32501 


3 


1 1.20l 


M 


1 0.32481 


5 


1 0.06 


(6) 




0.020 0.020 


0.020 


2 


2 


0.32321 


3 


1 1.18( 


; 46 ; 


1 0.32291 


; s ; 


1 0.06 


\i 




0.030 0.030 


0.030 


-1 





0.38911 


4 


1 1.27( 


; 47 : 


1 0.38901 


4 


1 0.06 






0.030 0.030 


0.030 


-1 


1 


0.38961 




1 1.281 




i 0.38951 




1 0.06 






0.030 0.030 


0.030 


-1 


-1 


0.389H 


$ 


1 1.271 




! 0.38901 


& 


1 0.06 






0.030 0.030 


0.030 








0.38881 




1 1.271 




l 0.38871 




1 0.06 


\i 




0.030 0.030 


0.030 


-2 





0.38991 


& 


i 1.281 




1 0.38981 


i 


1 0.06 






0.030 0.030 


0.030 


-2 


1 


0.39071 


4 


1 1.3H 


'48; 


1 0.39061 


'A" 


1 0.06 


(5) 




0.030 0.030 


0.030 


-2 


-1 


0.38971 


% 


1 1.27l 




1 0.38961 




1 0.06 


5) 




0.030 0.030 


0.030 


2 


-2 


0.39201 


% 


1 1.341 


'48, 


1 0.39191 


>4 


1 0.05 


(5) 




0.030 0.030 


0.030 


2 


2 


0.39011 


*4 


1 1.271 


47< 


l 0.39001 


> 4 < 


1 0.06 


(5) 
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TABLE XV: Summary of pseudo-scalar meson mass-squared 



splittings (xlO 3 ) obtained from the masses in Tab. XIII Lat- 
tice size 16 3 . 

m sea m 1 m 3 qi{\) qz{\) AAf 2 (xl0 3 )(cov) AM 2 (xl0 3 )(uncov) 

0.1210(14) 0.1207(16) 

0.1515 13) 0.1510(14 

0.1820 14) 0.1812(14 

0.1532(13) 0.1529(13) 

0.1844 13) 0.1838(12 

0.1869(13) 0.1865(12 

0.1222(15) 0.1222(16 

0.1236(17) 0.1237(17 

0.1551 13) 0.1549(13 

0.1596(27) 0.1583(31 

0.1800 25) 0.1788(28 

0.2013(27) 0.2001(29 

0.2005 24) 0.1997(26) 

0.2223(25) 0.2216(26 

0.2449(26) 0.2444(26) 

0.3244(43) 0.3245(47 



0.010 0.010 


0.010 





1 


0.010 0.010 


0.020 





1 


0.010 0.010 


0.030 





1 


0.010 0.020 


0.020 





1 


0.010 0.020 


0.030 





1 


0.010 0.030 


0.030 





1 


0.010 0.010 


0.020 


1 





0.010 0.010 


0.030 


1 





0.010 0.020 


0.030 


1 





0.010 0.010 


0.010 


1 


1 


0.010 0.010 


0.020 


1 


1 


0.010 0.010 


0.030 


1 


1 


0.010 0.020 


0.020 


1 


1 


0.010 0.020 


0.030 


1 


1 


0.010 0.030 


0.030 


1 


1 


0.010 0.010 


0.010 


1 


-1 


0.010 0.010 


0.020 


1 


-1 


0.010 0.010 


0.030 


1 


-1 


0.010 0.020 


0.020 


1 


-1 


0.010 0.020 


0.030 


1 


-1 


0.010 0.030 


0.030 


1 


-1 


0.010 0.010 


0.020 


-2 


-1 


0.010 0.010 


0.020 


1 


2 


0.010 0.010 


0.030 


-2 


-1 


0.010 0.010 


0.030 


1 


2 


0.010 0.020 


0.020 


-2 


-1 


0.010 0.020 


0.030 


-2 


-1 


0.010 0.020 


0.030 


1 


2 


0.010 0.030 


0.030 


-2 


-1 


0.010 0.010 


0.010 


2 





0.010 0.010 


0.020 


2 





0.010 0.010 


0.030 


2 





0.010 0.020 


0.020 


2 





0.010 0.020 


0.030 


2 





0.010 0.030 


0.030 


2 





0.010 0.010 


0.010 


2 


1 


0.010 0.010 


0.010 


2 


-1 


0.010 0.010 


0.020 


2 


-1 


0.010 0.010 


0.020 


-1 


2 


0.010 0.010 


0.030 


2 


-1 


0.010 0.010 


0.030 


-1 


2 


0.010 0.020 


0.020 


2 


-1 


0.010 0.020 


0.030 


2 


-1 


0.010 0.020 


0.030 


-1 


2 


0.010 0.030 


0.030 


2 


-1 


0.010 0.010 


0.020 





-2 


0.010 0.010 


0.030 





-2 


0.010 0.020 


0.030 





-2 


0.010 0.010 


0.010 


-2 


-2 


0.010 0.010 


0.020 


-2 


-2 


0.010 0.010 


0.030 


-2 


-2 


0.010 0.020 


0.020 


-2 


-2 


0.010 0.020 


0.030 


-2 


-2 



0.3674(43) 0.3675(44 
0.4100(44) 0.4095(44 



0.4123(42) 0.4120(41 
0.4566(43) 0.4558(41 



0.5027(43) 0.5015(41 
0.4582(68) 0.4562(75 



0.5448(57) 0.5416(62 

0.4730 76) 0.4716(82 

0.6460 59) 0.6421(61) 

0.5585(58) 0.5565(61) 

0.5747 63) 0.5734(65 

0.6614(59) 0.6590(59 

0.6797(62) 0.6783(61 

0.4899(59) 0.4885(65 

0.4947 61) 0.4944(65 

0.5002(68) 0.5003(70 

0.6176(53) 0.6164(52 

0.6252(55) 0.6244(53 

0.7513 54) 0.7495(49 

0.4456(66) 0.4425(75 

0.775(10) 0.775(11) 

0.833(10 0.834 10 

0.920(10) 0.919(10 

0.891(10) 0.891(10 



1.064(10) 1.061(10 
0.9828(97) 0.9820(95 



1.0441(98) 1.0426(94 
1.1309(99) 1.1283(95 



1.1964(98) 1.1934(93 

0.6106 54) 0.6086(56) 

0.7317(57) 0.7284(57 

0.7412(53) 0.7389(50) 

0.648(10) 0.642(12 

0.728(10) 0.723(11 

0.812(11) 0.807 11 

0.8092(98) 0.805(10 

0.895(10) 0.892(10) 
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TABLE XV - continued from previous page 
m sea mi m 3 qi{\) q 3 U) AM 2 (xl0 3 )(cov) AM 2 (x 10 3 )(uncov) 



0.010 0.030 


0.030 


-2 


-2 


0.010 0.010 


0.010 


-2 


2 


0.010 0.010 


0.020 


-2 


2 


0.010 0.010 


0.030 


-2 


2 


0.010 0.020 


0.020 


-2 


2 


0.010 0.020 


0.030 


-2 


2 


0.010 0.030 


0.030 


-2 


2 


0.020 0.010 


0.010 





1 


0.020 0.010 


0.020 





1 


0.020 0.010 


0.030 





1 


0.020 0.020 


0.020 





1 


0.020 0.020 


0.030 





1 


0.020 0.030 


0.030 





1 


0.020 0.010 


0.020 


1 





0.020 0.010 


0.030 


1 





0.020 0.020 


0.030 


1 





0.020 0.010 


0.010 


1 


1 


0.020 0.010 


0.020 


1 


1 


0.020 0.010 


0.030 


1 


1 


0.020 0.020 


0.020 


1 


1 


0.020 0.020 


0.030 


1 


1 


0.020 0.030 


0.030 


1 


1 


0.020 0.010 


0.010 


1 


-1 


0.020 0.010 


0.020 


1 


-1 


0.020 0.010 


0.030 


1 


-1 


0.020 0.020 


0.020 


1 


-1 


0.020 0.020 


0.030 


1 


-1 


0.020 0.030 


0.030 


1 


-1 


0.020 0.010 


0.020 


-2 


-1 


0.020 0.010 


0.020 


1 


2 


0.020 0.010 


0.030 


-2 


-1 


0.020 0.010 


0.030 


1 


2 


0.020 0.020 


0.020 


-2 


-1 


0.020 0.020 


0.030 


-2 


-1 


0.020 0.020 


0.030 


1 


2 


0.020 0.030 


0.030 


-2 


-1 


0.020 0.010 


0.010 


2 





0.020 0.010 


0.020 


2 





0.020 0.010 


0.030 


2 





0.020 0.020 


0.020 


2 





0.020 0.020 


0.030 


2 





0.020 0.030 


0.030 


2 





0.020 0.010 


0.010 


2 


1 


0.020 0.010 


0.010 


2 


-1 


0.020 0.010 


0.020 


2 


-1 


0.020 0.010 


0.020 


-1 


2 


0.020 0.010 


0.030 


2 


-1 


0.020 0.010 


0.030 


-1 


2 


0.020 0.020 


0.020 


2 


-1 


0.020 0.020 


0.030 


2 


-1 


0.020 0.020 


0.030 


-1 


2 


0.020 0.030 


0.030 


2 


-1 


0.020 0.010 


0.020 





-2 


0.020 0.010 


0.030 





-2 


0.020 0.020 


0.030 





-2 


0.020 0.010 


0.010 


-2 


-2 


0.020 0.010 


0.020 


-2 


-2 


0.020 0.010 


0.030 


-2 


-2 



0.984(10 
1.309(17 
1.480(17 
1.650(18 
1.659 17 
1.835(17 
2.019(17) 
0.1292(15) 
0.1577(13) 
0.1874(14) 
0.1588(13) 
0.1893 13 
0.1913(13 
0.1296 16 
0.1304(18 
0.1601(14 
0.1606(47 
0.1788 29 
0.1990(30 
0.1980 26 
0.2190 26 
0.2410(25 
0.3563 43 
0.3957(42 
0.4366 43 
0.4370(41 
0.4797 41 
0.5243(42 
0.4642 78 
0.5478(64 
0.4766 85 
0.6457 62 
0.5590(63 
0.5733 94 
0.6596(59 
0.6765 60 
0.5227(61 
0.5240 66 
0.5274(75 
0.6398 54 
0.6453(57 
0.7691 53 
0.4556(75 
0.847(10) 
0.898(10) 
0.9822(96) 
0.952(10) 
1.1218(99) 
1.0376(94 
1.0954(96 
1.1819 96 
1.2440(96 
0.6357(54 
0.7535 56 
0.7609(53 
0.651(12) 
0.723(11) 
0.803(12) 



0.982(10 
1.310(19 
1.481(17 
1.648(17 
1.658 16 
1.832(16 
2.014(16) 
0.1289(16) 
0.1575(13 
0.1875(14 
0.1586(13) 
0.1893(13 
0.1915(13 
0.1290(17 
0.1299(19 
0.1601(14 
0.1595(34 
0.1779(31 
0.1983(31 
0.1975(27 
0.2188(26 
0.2408(25 
0.3558(45 
0.3950(43 
0.4364(44 
0.4368(41 
0.4800(42 
0.5250(42 
0.4615(84 
0.5461(67 
0.4740(90 
0.6449(64 
0.5578(65 
0.5725(67 
0.6593(60 
0.6763(60 
0.5212(64 
0.5217(70 
0.5252(78 
0.6392(55 
0.6451(58 
0.7697(53) 
0.4531(82) 
0.846(10) 
0.896(10) 
0.9811(99) 
0.950(10) 
1.122(10) 
1.0372(96 
1.0959(97 
1.1828(97 
1.2456(97 
0.6349(56 
0.7538(57 
0.7612(53 
0.647(13) 
0.719(12) 
0.800(12) 
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TABLE XV - continued from previous page 
m sea mi m 3 <?i(±) q 3 (±) AM 2 (xl0 3 )(cov) AM 2 (x 10 3 )(uncov) 

0.799(10 
0.882(10 
0.968(19 
1.436(17 
1.593(16 
1.756(17 
1.757(16 
1.928(16 
2.105(16) 
0.1292(15 
0.1588(13 
0.1886 13 
0.1599(13 
0.1907(13 
0.1932(13 
0.1291 16 
0.1297(18 
0.1615 14 
0.1590(31 
0.1790 28 
0.1995 28 
0.2008(25 
0.2228 24 
0.2462(24 
0.3579 46 
0.3968(44 
0.4371 45 
0.4387(42 
0.4815 42 
0.5265(42 
0.4629 75 
0.5510 71 
0.4751(80 
0.6500 59 
0.5659(59 
0.5824 62 
0.6690(56 
0.6889 57 
0.5229(63 
0.5224 66 
0.5247(73' 
0.6446 53) 
0.6509 56) 
0.7767(53) 
0.4527(74) 
0.850(10) 
0.899(10 
0.987(10) 
0.950 10 
1.126(10 
1.0424(96 
1.1003(98 
1.1871(97 
1.2503 97 
0.6403(55 
0.7586 55 
0.7669(52 
0.645(12) 



0.020 0.020 


0.020 


-2 


-2 


0.020 0.020 


0.030 


-2 


-2 


0.020 0.030 


0.030 


-2 


-2 


0.020 0.010 


0.010 


-2 


2 


0.020 0.010 


0.020 


-2 


2 


0.020 0.010 


0.030 


-2 


2 


0.020 0.020 


0.020 


-2 


2 


0.020 0.020 


0.030 


-2 


2 


0.020 0.030 


0.030 


-2 


2 


0.030 0.010 


0.010 





1 


0.030 0.010 


0.020 





1 


0.030 0.010 


0.030 





1 


0.030 0.020 


0.020 





1 


0.030 0.020 


0.030 





1 


0.030 0.030 


0.030 





1 


0.030 0.010 


0.020 


1 





0.030 0.010 


0.030 


1 





0.030 0.020 


0.030 


1 





0.030 0.010 


0.010 


1 


1 


0.030 0.010 


0.020 


1 


1 


0.030 0.010 


0.030 


1 


1 


0.030 0.020 


0.020 


1 


1 


0.030 0.020 


0.030 


1 


1 


0.030 0.030 


0.030 


1 


1 


0.030 0.010 


0.010 


1 


-1 


0.030 0.010 


0.020 


1 


-1 


0.030 0.010 


0.030 


1 


-1 


0.030 0.020 


0.020 


1 


-1 


0.030 0.020 


0.030 


1 


-1 


0.030 0.030 


0.030 


1 


-1 


0.030 0.010 


0.020 


-2 


-1 


0.030 0.010 


0.020 


1 


2 


0.030 0.010 


0.030 


-2 


-1 


0.030 0.010 


0.030 


1 


2 


0.030 0.020 


0.020 


-2 


-1 


0.030 0.020 


0.030 


-2 


-1 


0.030 0.020 


0.030 


1 


2 


0.030 0.030 


0.030 


-2 


-1 


0.030 0.010 


0.010 


2 





0.030 0.010 


0.020 


2 





0.030 0.010 


0.030 


2 





0.030 0.020 


0.020 


2 





0.030 0.020 


0.030 


2 





0.030 0.030 


0.030 


2 





0.030 0.010 


0.010 


2 


1 


0.030 0.010 


0.010 


2 


-1 


0.030 0.010 


0.020 


2 


-1 


0.030 0.010 


0.020 


-1 


2 


0.030 0.010 


0.030 


2 


-1 


0.030 0.010 


0.030 


-1 


2 


0.030 0.020 


0.020 


2 


-1 


0.030 0.020 


0.030 


2 


-1 


0.030 0.020 


0.030 


-1 


2 


0.030 0.030 


0.030 


2 


-1 


0.030 0.010 


0.020 





-2 


0.030 0.010 


0.030 





-2 


0.030 0.020 


0.030 





-2 


0.030 0.010 


0.010 


-2 


-2 



0.797( 




0.880( 


ll\ 


0.968( 


l°i 


1.434( 




1.591( 


17) 


1.755( 


17) 


1.757 


16) 


1.929( 


16) 


2.108(17) 


0.1289 


16 


0.1584 


14 


0.1882 


(14) 


0.1594 


14 


0.1901 


(13) 


0.1924 


13 


0.1289 


17 


0.1295 


19 


0.1610 


(14) 


0.1587 


33 


0.1790 


(30 


0.1998 


(30 


0.2010 


(27) 


0.2232 


(26) 


0.2465 


25 


0.3569 


49 


0.3955 


47 


0.4354 




0.4367 




0.4788 




0.5232 


(18 


0.4627 


(80 


0.5504 


(65 


0.4759 


(86) 


0.6499 




0.5659 


in 


0.5828 




0.6690 


83) 


0.6889 




0.5215 




0.5213 




0.5238 




0.6427 


(56 


0.6489 


(58 


0.7738 


(55 


0.4517(80) 


0.848( 


11) 


0.896( 


11) 


0.984( 


10) 


0.947( 


11 


1.121( 


10) 


1.037( 


10) 


1.094( 


10) 


1.181( 


10) 


1.243( 


10) 


0.6386(57) 


0.7567 




0.7643 


n 


0.644(13) 
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TABLE XV - continued from previous page 
m sea m 1 m 3 q 3 (±) AM 2 (xlQ 3 )(cov) AM 2 (x 10 3 )(uncov) 



0.030 0.010 0.020 


-2 


-2 


0.724i 




0.7241 


12) 


0.030 0.010 0.030 


-2 


-2 


0.805i 




0.806( 


12 


0.030 0.020 0.020 


-2 


-2 


O.8IO1 


10 


0.8111 




0.030 0.020 0.030 


-2 


-2 


0.8978(98) 


0.8991 




0.030 0.030 0.030 


-2 


-2 


0.9899(96) 


0.99K 




0.030 0.010 0.010 


-2 


2 


1.442(18) 


1.4381 




0.030 0.010 0.020 


-2 


2 


1.598i 




1.5921 




0.030 0.010 0.030 


-2 


2 


1.758i 


IS 


1.7511 




0.030 0.020 0.020 


-2 


2 


1.765i 


16 


1.7561 




0.030 0.020 0.030 


-2 


2 


1.935i 


17) 


1.9241 




0.030 0.030 0.030 


-2 


2 


2.114i 


,17) 


2.101(17) 
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TABLE XVI: Summary of pseudo-scalar meson mass-squared 



splittings (xlO 3 ) obtained from the masses in Tab. XIV Lat- 
tice size 24 3 . 



0.005 0.005 


0.020 


-1 





0.1148(11) 


0.1147(11) 


0.005 0.005 


0.030 


-1 





0.1161 15 


0.1166(15 


0.005 0.010 


0.010 


-1 





0.12945 


(77 


) 0.12939 




0.005 0.010 


0.020 


-1 





0.13090 


(92 


) 0.13096 


m 


0.005 0.010 


0.030 


-1 





0.1327(10) 


0.1331(10) 


0.005 0.020 


0.020 


-1 





0.16304(85 


) 0.16335 


(76) 


0.005 0.020 


0.030 


-1 





0.16615 


(87 


) 0.16638 


80 


0.005 0.030 


0.030 


-1 





0.19980 


(86 


) 0.19990 


80 


0.005 0.001 


0.001 


-1 





0.10228(85 


) 0.10235 


96 


0.005 0.005 


0.005 


-1 


1 


0.3233 


22) 


0.3226(24) 


0.005 0.005 


0.010 


-1 


1 


0.3447 


23) 


0.3440(25 


0.005 0.005 


0.020 


-1 


1 


0.3876 


27j 


0.3878(27) 


0.005 0.005 


0.030 


-1 


1 


0.4310 


31 


0.4325(30 


0.005 0.010 


0.010 


-1 


1 


0.3665 


24 


0.3662(24 


0.005 0.010 


0.020 


-1 


1 


0.4109 


27) 


0.4114(26) 


0.005 0.010 


0.030 


-1 


1 


0.4558 


30) 


0.4572(28 


0.005 0.020 


0.020 


-1 


1 


0.4583 


29 


0.4592(26) 


0.005 0.020 


0.030 


-1 


1 


0.5062 


30) 


0.5074(28 


0.005 0.030 


0.030 


-1 


1 


0.5566 


'31 


0.5578(29 


0.005 0.001 


0.001 


-1 


1 


0.2880 




0.2889( 


27) 


0.005 0.005 


0.005 


-1 


-1 


0.1332 


n 


0.1326(16) 


0.005 0.005 


0.010 


-1 


-1 


0.1416 


i 


0.1414(15 


0.005 0.005 


0.020 


-1 


-1 


0.1603 




0.1607(17) 


0.005 0.005 


0.030 


-1 


-1 


0.1818 


20 


0.1821( 


20) 


0.005 0.010 


0.010 


-1 


-1 


0.1511 


14 


0.1512( 


14 


0.005 0.010 


0.020 


-1 


-1 


0.1713 


15) 


0.1718(14) 


0.005 0.010 


0.030 


-1 


-1 


0.1936 




0.1938(15 


0.005 0.020 


0.020 


-1 


-1 


0.1937 


11] 


0.1940(12) 


0.005 0.020 


0.030 


-1 


-1 


0.2176 


i 


0.2173(13 


0.005 0.030 


0.030 


-1 


-1 


0.2425 




0.2417(13) 


0.005 0.001 


0.001 


-1 


-1 


0.1208 


18) 


0.1202(20) 


0.005 0.005 


0.005 





1 


0.11418(75 


) 0.11384(82) 


0.005 0.005 


0.010 





1 


0.12905 


(73 


) 0.12898 




0.005 0.005 


0.020 





1 


0.15917 


(92 


) 0.15954 


85 


0.005 0.005 


0.030 





1 


0.1903 


ID 


0.1907(10) 


0.005 0.010 


0.020 





1 


0.16032(86) 0.16069(78) 


0.005 0.010 


0.030 





1 


0.1920(10) 


0.19246 87 


0.005 0.020 


0.030 





1 


0.19585(90 


) 0.19605 79 


0.005 0.005 


0.010 


1 





0.11413 


(85 


) 0.11380 


(91) 


0.005 0.005 


0.010 


-2 





0.4624 


34) 


0.4611( 


36) 


0.005 0.005 


0.020 


-2 





0.4653 


! 46) 


0.4650( 


46 


0.005 0.005 


0.030 


-2 





0.4705 


59) 


0.4723( 


i 


0.005 0.010 


0.010 


-2 





0.5232 


31 


0.5231( 




0.005 0.001 


0.001 


-2 





0.4156 


34) 


0.4159( 


38) 


0.005 0.010 


0.020 


-2 





0.5290 


*3 6 j 


0.5293( 


35j 


0.005 0.010 


0.030 


-2 





0.5365 


42 


0.5379( 


40 


0.005 0.020 


0.020 


-2 





0.6568 


34 


0.6580( 


30 


0.005 0.020 


0.030 


-2 





0.6691 


34) 


0.6701( 


32) 


0.005 0.030 


0.030 


-2 





0.8028 


ii 


0.8032( 


32) 


0.005 0.005 


0.005 


-2 


1 


0.7671 




0.7654( 


56 


0.005 0.005 


0.010 


-2 


1 


0.7947 


54) 


0.7929( 


i 


0.005 0.005 


0.010 


1 


-2 


0.8391 


53 


0.8381( 




0.005 0.005 


0.020 


-2 


1 


0.8519(66) 


0.8519(67) 
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TABLE XVI - continued from previous page 
m sea mi m 3 <?i(|) <? 3 (g) AM 2 (xl0 3 )(cov) AM 2 (x 10 3 )(uncov) 

0.9837(62) 0.9850(62) 

0.9102(80) 0.9137(78) 

1.1305(73) 1.1342(71 

0.8683(55) 0.8676(57 

0.6851 55) 0.6869(62 

0.9291(63) 0.9299(61 

1.0165 63) 1.0183(60 

0.9910(70) 0.9941(66 

1.1670(70) 1.1705(66) 

1.0848(66) 1.0868(61 

1.1539(68) 1.1565(64) 

1.2421 70) 1.2447(65 



0.005 0.005 


0.020 


1 


-2 


0.005 0.005 


0.030 


-2 


1 


0.005 0.005 


0.030 


1 


-2 


0.005 0.010 


0.010 


-2 


1 


0.005 0.001 


0.001 


-2 


1 


0.005 0.010 


0.020 


-2 


1 


0.005 0.010 


0.020 


1 


-2 


0.005 0.010 


0.030 


-2 


1 


0.005 0.010 


0.030 


1 


-2 


0.005 0.020 


0.020 


-2 


1 


0.005 0.020 


0.030 


-2 


1 


0.005 0.020 


0.030 


1 


-2 


0.005 0.030 


0.030 


-2 


1 


0.005 0.005 


0.005 


-2 


-1 


0.005 0.005 


0.010 


-2 


-1 


0.005 0.005 


0.010 


1 


2 


0.005 0.005 


0.020 


-2 


-1 


0.005 0.005 


0.020 


1 


2 


0.005 0.005 


0.030 


-2 


-1 


0.005 0.005 


0.030 


1 


2 


0.005 0.010 


0.010 


-2 


-1 


0.005 0.001 


0.001 


-2 


-1 


0.005 0.010 


0.020 


-2 


-1 


0.005 0.010 


0.020 


1 


2 


0.005 0.010 


0.030 


-2 


-1 


0.005 0.010 


0.030 


1 


2 


0.005 0.020 


0.020 


-2 


-1 


0.005 0.020 


0.030 


-2 


-1 


0.005 0.020 


0.030 


1 


2 


0.005 0.030 


0.030 


-2 


-1 


0.005 0.005 


0.005 


2 





0.005 0.005 


0.010 





-2 


0.005 0.005 


0.020 





-2 


0.005 0.005 


0.030 





-2 


0.005 0.010 


0.020 





-2 


0.005 0.010 


0.030 





-2 


0.005 0.020 


0.030 





-2 


0.005 0.005 


0.005 


2 


-2 


0.005 0.005 


0.010 


2 


-2 


0.005 0.005 


0.020 


2 


-2 


0.005 0.005 


0.030 


2 


-2 


0.005 0.010 


0.010 


2 


-2 


0.005 0.010 


0.020 


2 


-2 


0.005 0.010 


0.030 


2 


-2 


0.005 0.020 


0.020 


2 


-2 


0.005 0.020 


0.030 


2 


-2 


0.005 0.030 


0.030 


2 


-2 


0.005 0.001 


0.001 


2 


-2 


0.005 0.005 


0.005 


2 


2 


0.005 0.005 


0.010 


2 


2 


0.005 0.005 


0.020 


2 


2 


0.005 0.005 


0.030 


2 


2 


0.005 0.010 


0.010 


2 


2 


0.005 0.010 


0.020 


2 


2 


0.005 0.010 


0.030 


2 


2 


0.005 0.020 


0.020 


2 


2 


0.005 0.020 


0.030 


2 


2 


0.005 0.030 


0.030 


2 


2 



1.3172(70) 1.3196(67 
0.3861(35) 0.3846(38 



0.3878(38) 0.3870(41 
0.4322 33) 0.4321(34 



0.3967(48) 0.3970(49 
0.5282(38) 0.5299(34 



0.4112(61) 0.4122(63 

0.6312 46) 0.6322(39 

0.4367 34) 0.4370(34) 

0.3500(42) 0.3489(48) 

0.4492(39) 0.4499(36 

0.5365(36) 0.5381(31 

0.4658 45) 0.4664(42 

0.6417(40) 0.6425(33 

0.5546 35) 0.5556(29 

0.5758(37) 0.5755(31 

0.6639(36) 0.6635(29 

0.6878(36) 0.6864(30 

0.4626 30) 0.4614(32) 

0.5217 29) 0.5214(30 

0.6413(37) 0.6428(34) 

0.7649(46) 0.7667(41 



0.6458(34) 0.6474(31 
0.7718(40) 0.7735(35 



0.7870(36) 0.7878(32 
1.3062 90) 1.3034(99) 



1.3910(94) 1.388(10) 
1.561(10) 1.563 11 



1.734(12) 1.741(12 

1.4780(98) 1.476 10) 

1.654(11) 1.657 10) 

1.833(12 1.839(11) 

1.844(11 1.847 10 

2.034(12 2.039(11 

2.235(12) 2.240 12 

1.1658(95) 1.169(10 

0.5427 59) 0.5406(64 

0.5757(59) 0.5753(62 

0.6498 68) 0.6516(67 

0.7349(80) 0.7361(79 

0.6135 57) 0.6142(57) 

0.6934 61) 0.6953(55 

0.7816(68) 0.7823(60) 

0.7820 59) 0.7834(50 

0.8764(62) 0.8754(52) 

0.9747(62) 0.9718(52) 
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TABLE XVI - continued from previous page 
m sea m 1 m 3 gi(l) ggg) AM 2 (xlQ 3 )(cov) AM 2 (x 10 3 )(uncov) 



0.005 0.001 


0.001 


2 


2 


0.4943(70) 


i 0.4920(80) 


0.010 0.010 


0.010 





1 


0.13426(85 0.13325 75 


0.010 0.001 


0.001 


-1 





0.1044(10) 


i 0.1068(10) 


0.010 0.010 


0.020 





1 


0.16591(89 0.16441 80 


0.010 0.001 


0.001 


-1 


1 


0.3012(32^ 


i 0.3008(35) 


0.010 0.010 


0.030 





1 


0.19850(99) 0.19621(91) 


0.010 0.001 


0.001 


-1 


-1 


0.1162(17) 


i 0.1262(20) 


0.010 0.020 


0.020 





1 


0.16935(82 


) 0.16667(80) 


0.010 0.020 


0.030 





1 


0.20233 


(88 


) 0.19932 


(88) 


0.010 0.030 


0.030 





1 


0.20593 


(92 


) 0.20282 


90 


0.010 0.010 


0.020 


1 





0.13724 


(92 


) 0.13465 


(81) 


0.010 0.010 


0.030 


1 





0.1403(10) 


i 0.13637 


95 


0.010 0.020 


0.030 


1 





0.17274(88 


) 0.16931 


86 


0.010 0.010 


0.010 


1 


1 


0.158K 




i 0.15621 




0.010 0.010 


0.020 


1 


1 


0.18061 


15 


i 0.1762( 


15 


0.010 0.010 


0.030 


1 


1 


0.20321 


17 


i 0.19751 


16 


0.010 0.020 


0.020 


1 


1 


0.2022( 


15 


i 0.19751 


14 


0.010 0.020 


0.030 


1 


1 


0.2247( 


17 


i 0.220K 


; is ) 


0.010 0.030 


0.030 


1 


1 


0.24741 


19 


i 0.24391 


15 


0.010 0.010 


0.010 


1 


-1 


0.37871 


'25) 


i 0.37671 


25) 


0.010 0.010 


0.020 


1 


-1 


0.42571 


'27 


i 0.4218( 


'27 


0.010 0.010 


0.030 


1 


-1 


0.47461 


'30, 


i 0.46751 


'29) 


0.010 0.020 


0.020 


1 


-1 


0.47511 


\2T 


i 0.469K 


; 28 ) 


0.010 0.020 


0.030 


1 


-1 


0.52531 


M 


i 0.5170( 


)30j 


0.010 0.030 


0.030 


1 


-1 


0.57631 


;so; 


i 0.5672( 




0.010 0.010 


0.020 


-2 


-1 


0.47461 


)4o: 


i 0.4624( 


38 


0.010 0.010 


0.020 


1 


2 


0.55991 


;37; 


i 0.55081 


32 


0.010 0.010 


0.030 


-2 


-1 


0.49361 


45 


i 0.4767( 


44 


0.010 0.010 


0.030 


1 


2 


0.66631 


;3s; 


i 0.6543( 


35 


0.010 0.020 


0.020 


-2 


-1 


0.57791 


)3e: 


i 0.56581 


33 


0.010 0.020 


0.030 


-2 


-1 


0.59681 


'40) 


i 0.5837( 


36) 


0.010 0.020 


0.030 


1 


2 


0.68461 


'38 


i 0.6727( 


'34 


0.010 0.001 


0.001 


-2 





0.42391 


'40, 


i 0.43381 


'41) 


0.010 0.030 


0.030 


-2 


-1 


0.70391 


>44 


i 0.69381 


! 36 


0.010 0.010 


0.010 


2 





0.54261 




i 0.5386( 


30 


0.010 0.010 


0.020 


2 





0.55461 




i 0.544K 


32 


0.010 0.010 


0.030 


2 





0.56711 




i 0.55101 


37 


0.010 0.020 


0.020 


2 





0.682K 




i 0.67141 


32 


0.010 0.020 


0.030 


2 





0.69561 




i 0.6819( 


34 


0.010 0.030 


0.030 


2 





0.82751 




i 0.81501 


'36) 


0.010 0.010 


0.010 


2 


1 


0.45561 


$ 


i 0.45081 


>) 


0.010 0.001 


0.001 


-2 


1 


0.71251 


'73) 


i 0.71501 


80) 


0.010 0.010 


0.010 


2 


-1 


0.8982( 


'58 


i 0.89251 


'58 


0.010 0.010 


0.020 


2 


-1 


0.96651 


'62 


i 0.95431 


'60) 


0.010 0.010 


0.020 


-1 


2 


1.05171 




i 1.04291 


'63) 


0.010 0.010 


0.030 


2 


-1 


1.03801 




i 1.0173( 


'66 


0.010 0.010 


0.030 


-1 


2 


1.21091 




i 1.1953( 


'70) 


0.010 0.020 


0.020 


2 


-1 


1.12511 




i 1.1100( 


'64 


0.010 0.020 


0.030 


2 


-1 


1.199K 




i 1.17851 


'68) 


0.010 0.020 


0.030 


-1 


2 


1.28711 




i 1.26771 


'70 


0.010 0.030 


0.030 


2 


-1 


1.36311 




i 1.34161 


'71) 


0.010 0.001 


0.001 


-2 


-1 


0.34311 


;! 


i 0.36531 


)47i 


0.010 0.010 


0.020 





-2 


0.66851 


'35) 


i 0.6624( 


32) 


0.010 0.010 


0.030 





-2 


0.79801 


'39 


i 0.7886( 


'36 


0.010 0.020 


0.030 





-2 


0.813K 


'35 


i 0.80101 


'35 


0.010 0.001 


0.001 


2 


-2 


1.214(12) 


1.216(14) 


0.010 0.001 


0.001 


2 


2 


0.4768(71) 


i 0.51561 


'81) 


0.010 0.010 


0.010 


-2 


-2 


0.6414(61) 


i 0.6339( 


!62) 
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TABLE XVI - continued from previous page 
m sea mi m 3 qiU) q 3 U) AM 2 (xl0 3 )(cov) AM 2 (x 10 3 )(uncov) 



0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.010 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 



0.010 
0.010 
0.020 
0.020 
0.030 
0.010 
0.010 
0.010 
0.020 
0.020 
0.030 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 



0.020 
0.030 
0.020 
0.030 
0.030 
0.010 
0.020 
0.030 
0.020 
0.030 
0.030 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.020 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 
0.030 



-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-2 
-1 
-1 
-1 
-2 
-2 
-2 
2 
2 
-1 
-1 
-1 
-2 
-2 
-2 
2 
2 



-2 
-2 
-2 
-2 
-2 
2 
2 
2 
2 
2 
2 

1 
-1 

1 
-1 
-2 
2 

1 

-1 

1 
-1 

-2 
2 



0.7300(63 
0.8197(68 
0.8159 61 
0.9051(68 
0.9947(75) 
1.529(10) 
1.716 10 
1.911(12 
1.913(11) 
2.112(11 
2.315(12) 
0.1752(12 
0.5002(36 
0.2007 22 
0.7058(48 
1.1807 83 
0.5811(53 
2.011(14) 
0.8104(88) 
0.21666(88) 
0.6226(34' 
0.2442(20 
0.8709 35 
1.4670(77 
0.7087(44) 
2.502(13) 
0.9830(81) 



0.7129(61 
0.7972(66 
0.7970(57 
0.8864(60 
0.9806(63) 
1.519(10) 
1.698(10) 
1.880(12) 
1.887(11) 
2.078 12) 
2.278(12) 
0.17318(67) 
0.4947(23) 
0.1978(12 
0.6976(26) 
1.1680(53 
0.5732(28 
1.9905(95 
0.7987(48 
0.21507(66) 
0.6140(24' 
0.2462(11 
0.8643(26 
1.4477(54 
0.7108(25 
2.4663(97 
0.9901(44) 



lattice size 


mf m p % 2 /dof m n y^/doi 


16 3 


0.010 0.7125(57) 0.70(85) 0.7122(57) 0.70(85) 


16 3 


0.020 0.7986(40) 1.8(1.3) 0.7982(40) 1.8(1.3) 


16 3 


0.030 0.8747(36) 2.2(1.5) 0.8742(36) 2.2(1.5) 


24 3 


0.005 0.6477(53) 0.85(94) 0.6474(53) 0.85(94) 


24 3 


0.010 0.7121(31) 0.20(46) 0.7118(32) 0.21(46) 


24 3 


0.020 0.8065(25) 0.82(92) 0.8060(25) 0.84(93) 


24 3 


0.030 0.8871(23) 0.53(72) 0.8865(23) 0.53(72) 



TABLE XVII. Proton and neutron masses on unitary points for QCD+QED configurations. The 
masses are from size 16 3 box source, point sink correlation functions for both 16 3 and 24 3 lattices. 
The fit range for 16 3 lattices is 5-10. The fit range for 24 3 lattices is 6-11. The % 2 /dof is from a 
covariant fit. 
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lattice size 




2/1 £ 

X /dof 


m n 


X 2 /dof 




0.010 0.010 0.020 0.7416(49) 


0.59(78) 


0.7562(43) 


0.83(92) 


16 3 


0.010 0.010 0.030 0.7684(45) 


0.77(88) 


0.7981(36) 


1 i 


'i 1 "\ 
i.i j 


16 3 


0.010 0.020 0.030 0.8086(36) 


1.7(1.2) 


0.8238(33) 


9 n 
z.u 


'1 

/ 


16 3 


0.020 0.010 0.020 0.7553(51) 


0.88(95) 


0.7698(46) 


1.0 


l.Z ) 




0.020 0.010 0.030 0.7825(46) 


1.0(1.0) 


0.8120(40) 


9 /I 
Z.4 


'1 

l.OJ 


16 3 


0.020 0.020 0.030 0.8230(38) 


2.2(1.4) 


0.8380(36) 


9 Q 
z.y 


'1 7"l 


16 3 


0.030 0.010 0.020 0.7721(60) 


1.6(1.2) 


0.7839(51) 


1 ^ 


'1 1 
i.ii 


16^ 


0.030 0.010 0.030 0.7988(55) 


1.6(1.3) 


0.8241(42) 


1 9 
i.z 


'1 1 
i.ii 


16^ 


0.030 0.020 0.030 0.8361(41) 


1.7(1.3) 


0.8496(38) 


1 7 


l.OJ 




0.005 0.005 0.010 0.6676(85) 


1.3(1.1) 


0.6747(73) 


1 n 

l.U 


l.UJ 


24 d 


0.005 0.005 0.020 0.6992(68) 


1.5(1.2) 


0.7225(51) 


i n 

l.U 


l.UJ 


24 3 


0.005 0.005 0.030 0.7279(59) 


1.5(1.2) 


0.7680(41) 


1 9 
l.Z 


J LI 


24 d 


0.005 0.010 0.020 0.7225(51) 


1.6(1.2) 


t-\ m— ann / A A \ 

0.7383(44) 


1.4 


1.2) 


24 d 


0.005 0.010 0.030 0.7502(44) 


1.7(1.3) 


0.7824(37) 


1.7 


1.3) 


24 3 


0.005 0.020 0.030 0.7928(33) 


2.4(1.5) 


0.8090(32) 


2.5 


1.6) 


24 3 


0.010 0.010 0.020 0.7304(77) 


0.70(96) 


0.7461(61) 


0.9 


1.1) 


24 3 


0.010 0.010 0.030 0.7575(71) 


0.63(94) 


0.7895(48) 


1.3 


1.2) 


24 3 


0.010 0.020 0.030 0.7980(46) 


1.0(1.2) 


0.8146(40) 


1.8 


1.5) 



TABLE XVIII. Proton and neutron masses on QCD configurations with non-degenerate u, d quark 
masses. The p-n masses are from wall source, point sink correlation functions for 16 3 and 24 3 
lattices. The fit range for 16 3 lattices is 5-10. The fit range for 24 3 lattices is 7-12. The x 2 /dof is 
from a covariant fit. 



82 



Appendix B: Partially quenched chiral perturbation theory framework 

The aim of this appendix is twofold: (1) to obtain all possible terms in the chiral La- 
grangian relevant to the kaon mass-squared at order 0(e 2 ) and 0(e 2 p 2 ) for the partially 
quenched SU(2) + kaon system, and (2) to derive the expression for the EM correction 
to the kaon mass-squared to order 0(e 2 p 2 ). The appendix is compact, mainly summa- 
rizing results and defining notation. Much of the machinery, of course, has been worked 
out before, and we refer the interested reader to the literature. Here we follow closely the 
works in Refs. [TUHT3] . The new contributions in this work are 0(e 2 )-terms and electro- 
magnetic one-loop chiral logarithmic correction to the kaon mass-squared. We also list the 
0(e 2 p 2 )-operators relevant to the kaon-mass-squared, which serves as a check of the possible 
dependence of 0(e 2 p 2 ) corrections on charges and masses. 

We begin by reminding the reader of the important details and notation, then construct 
the Lagrangian density, and finally compute the corrections to the kaon mass-squared to the 
order of our interest. 



1. SU(2) pion sector 

In the partially quenched system composed of Ny valence quarks, N$ sea quarks and 
Ny ghost quarks, the field U(x) representing the Nambu-Goldstone multiplet is the local 
coordinate of the coset space G/H (G = SXJ(N S + N V \N V ) L x SU(iV s + N V \N V ) R , H = 
SXJ(N S + N V \N V ) V ) at each x 

u[U{x)} = exp (i UM.) . (Bl) 



V2F 

With the normalization of F such that F ~ 92 MeV, the leading-order (LO) chiral La- 
grangian reads 

£qcd, 2 = ^ + , (B2) 

where (, ) denotes the supertrace in the partially quenched light quark sector, whose flavors 
can be indexed as 

1 = 1, ■ ■ ■ , Ny '■ light valence quark flavors , 

I = N v + 1, ■ ■ ■ , Ny + N s : light sea quark flavors , 

J = N v + N s + 1, ■ ■ ■ , 2N V + N s : light ghost quark flavors . (B3) 
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The variables appearing in Eq. (B2), 



Uft = i {u 1 (d^u - iR^u) - u (d^u 1 - iL^u ] ) } , 

X± = v) X u ] ± u X ] u , x = 2B M , (B4) 

are given in terms of the spurion field Ai in place of the ordinary quark mass matrix, and the 
external fields R^, L^, which transform under the local chiral rotation (g^x), g R (x)) G G 
as 

Rfi^- R'n = 9rRp9r + WR d^9R , 
L » ^ L 'n = 9LLf,g[ + ig L d M # z , 

M^M' = g R Mg{ . (B5) 

For 

u[U] ^ u[W] = g R u[U] h((g L , g R ); n)t = h((g L , g R ); IT) u[U] g\ , (B6) 

with h((gi, g R )'i n) G H, it turns out that and x± transform covariantly with respect to 
h 

A^A' = hAh ] , (B7) 

and that £qcd,2 is invariant under the local chiral transformation. 

The high frequency modes of photons coupled to quarks also generate local interactions 
in the low-energy effective Lagrangian of QCD. The coupling of quarks to photons preserves 
chiralities; 

A^q^QqL + qR^QqR) , (B8) 
where Q represents the charge matrix and takes the form 

Q = ediag(g ul/ , <?dv, q u s, QdS, QuV, Qdv) , (B9) 

for two-light flavors. The systematic dependence on these quark charges can hence be traced 
back once Q is promoted to a set of hermitian spurion fields, Q Ri l, that transform under 
chiral rotations as 

Ql^ Q' l =9lQl91, Qr^Q' r = 9rQr9 r . (BIO) 
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On the other hand, to construct the effective Lagrangian, it is convenient to define quantities 



that transform covariantly by h((g L , g R ); II), i.e., as in Eq. (B7) 



Q L = uQ L u\ Q R = u^Q R u. (Bll) 

Since we will set Ql, Qr to the diagonal EM charge matrix Q after constructing the effective 
Lagrangian, we impose the chiral-invariant condition 

(Qr) = (Ql) = (Q) , (B12) 

which reduces just to the charge matrix in the end. The leading order (0(p 2 ) ~ 0(e 2 )) 
Lagrangian involving Nambu-Goldstone bosons is thus given by 

L*,2 = Y ( u » ull + X+) + C (QrQl) , (B13) 

The QED corrections from the low frequency photons can also be included by coupling the 
Nambu-Goldstone bosons to the [/(l)-gauge potential A^(x) and by setting the external 
fields L M , i? M along the direction of Q in the end; 

L fl =QA lt = R lt . (B14) 



2. The kaon sector 



In SU(2) chiral perturbation theory, the strange quark is treated as being heavy, and 
hence the kaons are no longer treated as Nambu-Goldstone bosons. Since the EM charge 
of the sea strange quark, ss, differs from that of the valence strange quark, sy, in our 
simulation, these together with the ghost strange quark, s, are regarded as constituting the 
partially-quenched strange sector. Nevertheless, for the purpose of the analysis of our lattice 
data, it suffices to write down the effective Lagrangian with respect to the kaon multiplet 
including the valence anti-strange quark sy, keeping track explicitly of the dependence on 
the electric charges Q s ,v, Qs,s (including e) of sy and ss, respectively, with the low energy 
constants having implicit dependence on the sea strange quark mass. The relativistic form 
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of the kinetic and mass terms of the kaon multiplet K (U, D denote constituent quarks) 

( [Uvs v ] \ 

[D V Sv] 
[UsSv] 
[DsM 

[U G S V ] 

V [DgM J 

which is subject to the chiral rotation 



K 



(B15) 



K^h[U{x), (g L , g R )]K 



(B16) 



is given by 



C 



K, kin 



V»K - M 2 K ] K 



(B17) 



where M is the LO mass of the kaon and the covariant derivative V^K is with respect to 
the Maurer-Cartan form r„ 



r, 



2i 



{u ] (d^u - iR^u) + u (d^u ] - iL^u ] ) } . 



(B18) 



As is well-known [23], K is not suitable for the chiral order counting since that variable also 
carries the high frequency modes. The fluctuation is decomposed into the high frequency 
modes originating from M and the low frequency modes represented by k = k v 



K{x) 



i iMv - x k{x) 



(B19) 



where v is a light-like four-vector. In terms of k, Eq. (B17) becomes 



(B20) 



The field k carries the momentum of the order p < 4nF n , M, and the above Lagrangian is 
0(p). In the succeeding sections the effective Lagrangian is constructed in terms of k and 
is converted to the relativistic form described by K. 



86 



3. 0(e ) and 0(e p ) Lagrangian for the kaon sector 

Having established the partially quenched framework and notation, we construct the 
electromagnetic part of the chiral Lagrangian by writing down all possible 0(e 2 )-terms 
possessing the symmetries of massless (QCD + QED) in the non-relativistic theory and 
their relativistic counterparts, and 0(e 2 p 2 )-terms that can induce the tree-level contribution 
to the kaon mass-squared. 

building block definition order P C 



X ± Eq. \m\ 0(p 2 ) ±x±(x) ix±) 

T 



9 2 ) ± X ±(x) (x±f 
(Q) Q± (Q) (Qr ± Ql) 0(e 2 ) ± (Q) Q ± {x) ± (Q) (q± 
S(±) (Qr) 2 ± (Ql) 2 O(e') ±Q 2 ±) (x) ±(Q 2 ±) ) r 
Qrl,± QrQl±QlQr 0(e 2 ) ±Qrl,±(x) (q rl ,± 



TABLE XIX. Parity (P) and charge conjugation (C) transformation properties for operators at 
chiral order 0(p 2 ), 0(e 2 ) that do not contain kaon fields and transform as A \— > h A h) . Under P, 
x = (x°, x) transforms to x = — x). 



definition order P C 



ktf O(l) fcfct(5) (Jfefct) 



k±,it Eq. ( 
kfrv] Eq. ( 
k±,fiu Eq. ( 


B21 
B22 
B23 


) 0( P ) fc£(x) ±(fc±, M ) T 
)0( P 2 ) fcC^](S) ±(fc (H ) T 

)0(p 2 ) fcf(x) ±(£^) T 


kf®+ ^ ( 
kf Q - Eq.( 


B24 
B24 


io(e 2 ) fcj^ +(i) ± ^Q + y 

)0(e 2 ) -k^-(x)^{kf^) T 



TABLE XX. Parity (P) and charge conjugation (C) transformation properties of operators at 
chiral order 0(p 2 ), 0(e 2 ) that are bilinear in kaon fields and transform as A \— > h A h) . W is either 
one of Q s ,v or Q St s- Under P, x = (x , x) transforms to x = (x , — x). 



To this end, Table XIX and XX list the building blocks of chiral order 0(p 2 ) and 0(e 2 ) 
that transform as Eq. (B7). The definition of various variables appearing in Table XX are 
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as follows; 



k ± ^ = i{{V^ ±k{V^ , (B21) 

= - (V^kVX + V u kV^) 
k\ptv] = V[^kV v ]k^ 

= r (V^V.^ - V.W^) , (B22) 

fc±,^ = (V^A;)A;t ± A ; (v^A ; ) t , V M „ = V^V, + V„V M , (B23) 

fc w,a± = ^ ^ f Q ± ± Q±A;A; t ) . (B24) 



In Table XIX, u^, for instance, is omitted, as it will be not be used hereafter. Table XIX 



and XX include 0(e )-terms but not 0(e) and O(ep), because EM charges are left in pairs 



in the low-energy effective theory after the high frequency photon modes are integrated out. 



From Table XIX and XX, we find that there are 13 0(e )-operators bilinear in kaon fields 



that are invariant under chiral, P and C transformations 



(jfeifeU) , W(kk^B), WiW 2 (ktf) 



(B25) 



where W is Q s ,v or Q St s, and 



Ae |Q ( 2 +) , Q RL , + , (Q) Q+, (Q 2 ) , (Qrl,+) , («2» 

(W h W 2 ) G {(Q fl , y, Q s ,y) , (Q s ,v, Q SiS ) , (Q S)5 , Qs, S )} ■ 



(B26) 



The relativistic forms of the individual operators are read from the relation (B19), and 



8N 



0(e 2 )-Lagrangian density in the kaon sector is hence given by 



C 



K, e 2 



-A 



(1,2) 
K 



Qr) 2 +(Ql) 2 ) k^k 



-A%' 1] (q r Q l + q l q r ) K - Ag' 2) ( Q R Q L + Q L Q R )K*K 



-Af (Q) tft (q h + Q L ) K - A<? (Q) 2 K^K 

-4* 1} v K*K ~ 4' 11 2) 02, 5 ^ - 4 '* 3) Qs, v Qs, s K^K 

-4 3) (Q)gs,v^-4' 3,2) (2) Qs, S K*K 

-A% 3 ' 3) Q s , s kUq r + Q l )k. (B27) 



There are 0(e 2 p)-terms bilinear in kaon fields that are allowed from the symmetries. All 



possible terms are obtained from Eq. (B26 ) by the replacement kk) — » v ok-- These operators 



generate no 0(e 2 )-contribution and chiral-logarithmic corrections to the kaon mass-squared. 
They, however, induce 0(e 2 p 2 )-contribution to the kaon mass-squared after the renormal- 



ization of the kaon field. We shall come back to this point in Section B 4 



definition 



C 



Eq. 
Eq. 



k W, Q± E 

(H(l)>±(2) H ' 

h W,Q+ 

±(1)>±(2).P 
±(1),±(2),/^ 

fc ± (1) ,±( 2) ^ ^q- 

7 W,V„Q_ -p 



B28) C' ,H «> ® (CIUJ 



B29 



B29 



B30 



B30 



B31 



k 



W,Q+,nv 

±(1),±(2) 



(2) 



i(l),±(2) 
*±(1),±(2) 



kY:^:"(x) 



;± (1) i) (± (2) i) (^ + ±(2) , 
(± (2) i) 
i) (± (a) i) 



(i),±(2:, /"' 
Q 

).±(2)>P 



Z (l) 



±(2),M 



A- 



W,V v Q-,n 



(2) 



(± (1) l) (± (2) 1) (O^.M 



TABLE XXI. Parity (P) and charge conjugation (C) transformation properties of operators at 
chiral order 0(e 2 p 2 ) that are bilinear in kaon fields and transform as A i— > h A tv . 



Next we turn to listing up 0(e 2 p 2 )-terms that induce the corrections to the kaon mass- 
squared at the tree level. The building blocks are those in Table IXIXl IXXI and IXXIl The 
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definition of the quantities in Table XXI is as follows; 



k W,Q± 



<rt£ = W (*<HG±±fi±*(H 



k±Z±%, = W {^VuQ* ±(2) (v„Q ± ) fc ±(lH 



The 0(e 2 p 2 )-terms with no derivatives are 



(kk^{ X+ , C + }) , (kk^Q +x+ Q + ) , (ktf Q-X+Q- 
(ktf\x-,C-]) , (ktf{ X -, Qrl,-}) , 

k^k ( X+ C + ) , (kh*C+) ( X+ ) , <A^x + ) , fctfc ( X+ ) ( C+ ) 
' ktf Q + ) ( X+ Q+) , 



where 



lyfctfc , W (fcfct x+ ) (Q) , W (kk^Q + ) ( X+ ) , 

Wk^k( X+ )(Q) , 
W 1 W 2 (kk^ X+ ) , WW^A; (X+) , 



C + e{(Q)Q + , Q { 2 +) , Qrl,+| , C_g{(Q)Q_, Q ( 2 _ } } . 



(B28) 
(B29) 
(B30) 
(B31) 



(B32) 



(B33) 



The 0(e 2 p 2 )-terms with two derivatives that contribute to the kaon mass-squared at the 
tree level are (i] fJ-l/ is the metric in Minkowski space) 

W {rT V) B) , W (rT Q-) > W (*T S> > 

WW 2 (?T^)) , (v^k+^v) , (B34) 



where ^4, B are the same as in Eq. (B26) 
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4. EM correction to kaon mass-squared 



In this subsection, the explicit expression for the 0(e 2 ) and the 0(e 2 p 2 ) chiral-logarithmic 
correction to the kaon mass-squared is obtained by setting the charge matrices Q in light 
flavor partially quenched system as in Eq. (B9), Q s ,v = eg s y and Q s ,s = e<lss- The EM 



contribution to the K + mass-squared at order e 2 is 



(M 



2e\(A^ + A^'hlv + e\2Ar> + 2A™ + + q 2 dS ) 

+e 2 A%\q uS + q d s)quv 

-\-e a k q sV -+- e n K q sS -+- e n. K q s vqsS 



l(2,l)\„2 



(1.2) 



I (2,2) 



,2 A (s,3) 



+e 2 A^ 6, (q uS + g dS )g s y + e 2 A { ^' 2 \q uS + g d5 )g sS 
A K q uV q sS . 



In quenched QED this becomes 



(M^ + ) 2 = 2e\A^ + A^)^v + ^A^'qiv + 2e 2 A^ q uV q sV . 



l( 2 ,l)\„2 



2 /I MO 



,2 4 (». 2 ). 



(B35) 



(B36) 



The 0(e 2 )-correction to the neutral kaon mass-squared, (Mj?? ) 2 , is given by substituting q^y 



for q u y in Eq. (B35). 



We next consider the one-loop contribution to kaon mass squared. The scalar QED La- 



grangian density (B17) gives the correction from the diagrams, in each of which a photon 
propagates explicitly, but these contributions are absorbed by the redefinition of the coeffi- 
cients in Eq. (B27) and those of 0(e 2 p)-operators in the infinite volume. The leading EM 



chiral-logarithmic correction comes only from the tadpole diagrams induced by Eq. (|B27|) 



e 2 A 



K 



167T 2 


F 2 < 

u n 


e 2 


,(2,1) 


167T 2 


P2 £ 

^0 n 


e 2 


2A$*> 



( q iV ~ fns) Xin In 



Xin 

~7 



16tt 2 



Xin 

~7 



53 (<ZnS — Q'ms) 2 Xmn hi 



e 2 A®iV s g + ^^^ + ^ 



n, m : sea. n=£m 
(«,2) 



i (-.3, 3) 



167T 2 
X 



F 2 



5^ (?tv - ?ns) Xm In f j 

-» ■ son ^ ' 



(B37) 
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where i = u or d, fi is the renormalization scale, and 

Xm Xn 



Xr, 



X, 



2B m n 



Q 



1 

iVc 



(B38) 



In our simulation, all sea quarks are neutral and the two light sea quarks are degenerate in 
mass mts)- Hence (M 1 ^) 2 reduces to 



(M 



log \2 
K.i) 



— 2 



16vr 2 F 2 



Mil ^ 



(2,1) 
if 



QiVQsV A K 



Xi(s) In 



Xi(5) 



(B39) 



where Xi(S) = (mi + m {S )). 

There are two types of finite volume corrections induced at the one-loop level. The first 
type is given by the scalar QED diagrams 



A(M j 



{L) = {q K fe 2 



I EM, photonic 



K 1 1 
4^ + ^)2 



L 2 



(4tt)' 



T-L{rriKL) 



(B40) 



where « and various functions are defined in Eqs. (14), (16) and (17). Another type is the 

2 



finite volume correction to the terms |b3t|, A «,) (t), whose expression is obtained 

xithm in 

M{mL) 



by making the following substitution to each logarithm in Eq. (B37) 

,2 



m In 



m 



L 2 



(B41) 



with ftA(x) in Eq. (13). The finite size scaling effect on the 0(e ) wave function renormal- 
ization could induce 0(e 2 p 2 )-correction to kaon mass squared after the renormalization of 
the kaon field. The explicit calculation, however, shows that such effects do not exist. 



There are as many LEC's as 0(e p )-operators in Eqs. (B32) and (B34) participating in 
the 0(e 2 p 2 )-contribution to kaon mass-squared, while our lattice study here can determine 
at best the linear combinations of LEC's of terms with the same charge and light quark mass 
dependence of order e 2 m, from the response of the data to the variation of these parameters 
in the (QCD + QED) action. The dependence on those parameters can be read off from 



Eqs. (B32) and (B34). In effect, Eq. (B32) alone leads to the following form of the charge 
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and mass dependence of the 0(e 2 m)-correction (i = u, d) in quenched QED, as anticipated, 
M £ P i) = e 2 m iv (x ( 3 K) (ftv + Qsv) 2 + x{ K) (ftv - Qsvf + 4*° {Viv ~ llv) 



+e 2 m (S) (xf > + + ^ (?*v - dsv? + X™ - <&)) 



(B42) 



We note that rriiv an d nvsrj are denoted by m\ and m 4 = m 5 , respectively, in Eq. (23). 
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